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I.  INTRODUCTION 


In  order  to  apply  workability  theory  to  metalworking  processes,  we 
must  first  establish  flow  and  fracture  theories.  Then,  by  combining  the  de- 
tailed information  on  mechanics  and  the  ductile  fracture  criterion,  workabil- 
ity of  materials  in  metalworking  processes  can  be  predicted.  Based  on  this 
prediction,  workability  control  for  preventing  fracturing  can  be  achieved 
by  selecting  the  proper  set  of  process  parameters. 

In  the  previous  report  [l]*  theories  on  flow  and  fracture  in  metal- 
working processes  were  developed  with  an  emphasis  on  applying  the  workabil- 
ity theory  to  metalworking  processes  where  the  occurrence  of  internal 
fracturing  is  a limiting  factor.  The  present  investigation  is  concerned 
with  the  application  of  these  theories  to  the  prediction  of  workability 
of  materials  in  axisymmetric  bar  extrusion  and  drawing,  with  special  refer- 
ence to  center  bursting. 

Workability  of  materials  is  the  extent  to  which  materials  can  deform 
without  forming  cracks  during  a mechanical  working  process.  Workability, 
therefore,  depends  on  the  conditions  imposed  by  the  working  process. 

Critical  stress  and  strain  conditions  involved  in  the  mechanical,  working 
processes  must  be  known  and  should  be  specified  in  te.:"is  of  process  vari- 
ables, such  as  height  reduction,  friction  at  the  interface,  and  workpiece 
geometries  and  dimensions.  Part  I of  the  present  investigation  deals  with 
the  deteimination  of  deformation  mechanics  in  extrusion  and  drawing. 

The  role  of  the  detailed  mechanics  in  the  worl^ability  study  is  to 

provide  the  stress  and  strain  paths  at  a critical  site  of  a deforming 
•Numbers  in  brackets  refer  to  the  references  ett  the  end  of  this  report. 


®®terisl • Therefore,  'the  method  of  analysis  should  he  capable  of  accurately 
determining  not  only  overall  quantities  involved  in  metalworking  processes, 
such  as  forming  loads,  but  also  stress  and  strain  distributions  during 
deformation.  Furthermore,  it  is  required  to  determine  the  stress  and  strain 
distributions  under  various  process  conditions.  Therefore,  to  Justify  the 
approach,  the  computation  involved  in  the  method  must  be  efficient.  The 
matrix  method  developed  by  Lee  and  Kobayashi  [2]  comes  close  to  fulfilling 
these  reqiiirements.  In  the  previous  report,  this  matrix  method  was  refined 
and  it  was  demonstrated  that  the  method  is  effective  for  the  analysis  of 
steady-state  processes  as  well  as  non-steady-state  processes.  In  the  pre- 
sent investigation,  the  matrix  method,  with  further  improvements,  is  vised 
to  determine  steady-state  deformation  characteristics  as  functions  of  mater- 
ial property,  die-workpiece  interface  friction,  die  angle,  and  reduction 
in  bar  extrusion  and  drawing. 


II.  METHOD  OF  ANALYSIS 


The  basic  concepts  of  the  matrix  method  for  rigid-plastic  deformation 
problems  are  the  use  of  the  Lagrange  multiplier  in  a variational  formula- 
tion and  linesirization  of  nonlinear  stiffness  equations.  In  the  present 
formulation  the  dynamic  effects,  i.e.,  the  inertia  effects  on  the  forces 
and  the  strain-rate  effects  on  the  material  properties,  are  neglected. 

For  rigid-plastic  materials  the  condition  of  incompressibility  is 
imposed  on  the  admissible  velocity  fields.  This  constraint  can  be  removed 
by  introducing  a Lagrange  mviltiplier.  Consider  a body  V whose  surface  S, 
consists  of  Sy  and  S^.  The  body  is  composed  of  a rigid-plastic  material 
that  obeys  the  von  Mises  yield  criterion  and  its  associated  flow  rule, 
under  the  boundary  conditions,  such  that  the  entire  body  is  deforming 
plastically.  Body  forces  are  asstuned  to  be  absent  in  the  region  V.  It 
can  be  shown  [3]  that  for  the  actual  solution,  the  functional  (l)  becomes 
stationary  with  respect  to  the  multiplier  X and  the  velocity  fields  that 
satisfy  the  velocity  boundary  conditions  on  Sy  but  not  necessarily  the 
incompressibility  condition  (kinematically  complete): 


ae  dV  + I 

XC^’e  dV  - 

t'^U  dS, 

(1) 

iv  J 

V ~ J 

where  e is  the  effective  strain-rate;  5,  the  effective  stress;  T,  the 

traction  vector  specified  on  the  boundary  S^;  U,  the  velocity  vector; 

T» 

C is  the  proper  vector  notation  of  the  Kronecker  delta  such  that  C g = 0 
implies  the  incompressibility  condition. 


U 


The  formulation  of  the  discrete  variational  problem  follows  the  same 
procedure  as  that  used  in  the  finite-element  method  [U],  [5].  The  body  V 
is  divided  into  M elements  interconnected  at  N nodal  points.  The  approxima- 
tion of  the  functional  ♦ by  a function  $ is  performed  on  the  elemental  level 
by  replacing  U with  a kinematically  complete  distribution,  given  by 

y - 2iJ»  (2) 

where  G is  the  interpolation  function  euid  u is  the  vector  whose  components 
are  velocities  at  nodal  points  associated  with  the  element.  The  stredn- 
rate  vector  is  then  derivable  in  the  form 

e = (3) 

Assembling  the  function  at  the  elemental  level  into  an  approximate  finite- 
element  model  over  all  the  elements  and  applying  the  stationary  condition 
to  the  function  we  obtain  the  stiffness  relations  consisting  of  a large 
system  of  nonlinear  equations.  In  order  to  solve  the  stiffness  eq\iations, 
we  adopt  the  following  procedure.  The  nonlinear  stiffness  equations  were 
lineeurized  by  considering  a small  perturbation  Au  in  the  velocity  vector  u, 
such  that 


»(n)  ' »(n-l)  * ''»(»)• 

for  the  n-th  iteration  process."^  We  then  obtain  the  global  perturbation 
equation  for  the  n-th  iteration  as 


In  the  actual  calculations,  U/  » = u/  - v + oAu,  x is  used,  where  a Is 

~vn;  ~(n-lj  ~(n) 

the  deceleration  coefficient. 


^n-1) 


(n) 


(5) 


. ^ - 


“ R / T \ • 

•“(n-l) 


Specific  form\ilations  for  the  matrix  and  the  vector  R,  \ising  a quadri- 
latereil  element  with  a bilinear  velocity  distribution,  are  presented  else- 
idiere  [l],  [6].  It  must  be  noted  that,  since  the  method  is  an  iterative 
process,  it  is  necessary  to  provide  an  initial  guess  for  the  velocity  field 
Uq  (but  no  need  for  X).  The  solution  of  linear  equation  (5)  yields 
close  to  zero  eind  a proper  veilue  of  the  mean  stresses  ~(n-l) 

close  to  the  actual  solution. 

A.  Convergence 

In  the  previous  studies  the  convergence  of  the  solution  has  been 
measured  by  the  quantity  |Ay|/|u|,  where  the  Euclidean  vector  norm  is  de- 
fined by 

U = /f  Ib.I^  (6) 

i=l 

where  N = total  nxmiber  of  nodal  points. 

The  convergence  criterion  requires  that  the  error  norm  at  the  n-th 

iteration  (|Au/  \l/lu/  -xl)  be  less  than  that  at  the  previous  iteration 
~(Dj  n-l ) 

This  convergence  criterion,  with  suitable  selections  of  the  deceleration 
coefficient  a,  has  worked  out  well  for  the  solutions  of  various  metal- 
working problems.  At  the  same  time,  however,  it  was  realized  that  solution 
divergence  has  been  indicated  according  to  this  criterion,  although  the 
solution  was  converging  in  its  true  sense.  This  can  appau*ently  be  seen, 
particularly  for  a function  which  is  not  well  behaved.  In  this  case,  the 


6 


criterion  resulted  In  the  use  of  unnecessarily  small  values  of  the  decelera- 
tion coefficient,  and  thus  In  Increased  computing  time.  In  the  present 
program.  Instead  of  the  error  norm,  the  quantity  f defined  by 


(m)  1 2\ 


where 


Is  utilized  for  solution  convergence.  In  Eq.  (7),  the  superscript  (m) 
denotes  the  values  of  the  m-th  element  and  summation  Is  made  over  all  the 
elements.  The  criterion  for  convergence  now  Is  that  the  magnitude  of  f at 
the  n-th  Iteration  be  less  than  the  magnitude  at  the  previous  Iteration. 

In  this  convergence  scheme,  the  proper  value  of  the  deceleration  coefficient 
at  each  Iteration  can  be  selected  efficiently  from  the  previous  Information 
on  the  function  behavior  and  the  convergence  requirement  at  the  current 
Iteration. 


B.  Rigid-body  treatment 

The  matrix  method  described  In  this  section  applies  only  If  the  entire 
body  Is  deforming  plastically  and  no  rigid  zone  or  unloading  exists  during 
the  deformation  process.  In  practical  problems,  however,  situations  do  arise 
where  the  rigid  zone  as  well  as  rigid  unloading  are  Involved.  If  these 
regions  of  no  deformation  are  contained  within  the  control  volume  V,  the 
extremum  principles  do  not  apply  to  the  problem  of  obtaining  Internal 
distributions.  There  are  seviral  approaches  to  handling  this  difficulty. 
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but  a most  effective  technique  is  one  which  involves  the  approximate  deter- 
mination of  the  boundary  of  a nearly  rigid  zone. 

A nearly  rigid  zone  can  be  characterized  by  its  very  low  value  of 
effective  strain  rate  in  comparison  to  the  deforming  body.  Ihiring  the  iter- 
ation process  for  an  incremental  solution  over  the  entire  body,  the  effec- 
tive strain-rate  in  the  possible  rigid  region  approaches  zero  as  the  solu- 
tion converges.  Since  the  effective  strain- rate  appears  in  the  denominator 
of  the  stiffness  matrix  S,  the  component  of  the  normalized  stiffness  matrix 
will  tend  to  become  infinity  if  the  nodal  point  associated  with  this  compo- 
nent is  contained  in  the  rigid  zone.  Therefore,  the  elements  for  which  the 
effective  strain-rate  is  smaller  than  a certain  value  (say,  0.0001)  are 
considered  to  be  in  the  rigid  zone.  The  effective  strain- rates  of  the  ele- 
ments lying  inside  the  rigid  zone  are  then  kept  at  this  value  in  the  per- 
turbation relationship,  and  the  iteration  is  continued  for  the  solution  in 
the  plastically  deforming  region  until  a desired  convergence  is  achieved. 

It  should  be  noted,  of  course,  that  the  converged  solution  gives  the  stress 
distribution  only  in  the  plastically  deforming  region. 

The  complete  program  of  this  improved  version  of  the  matrix  method 
is  listed  in  the  appendix  of  this  report. 


Several  investigators  have  analyzed  the  extrusion  process  by  using 
the  finite-element  method,  mostly  by  elastic-plastic  analysis  [7],  [8],  [9]* 
However,  almost  all  the  analyses  have  been  performed  for  the  case  of  loading 
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of  the  workpiece  that  fits  the  die  and  container  and  of  extruding  it  a 
small  amount,  instead  of  extruding  the  workpiece  until  a steady  state  is 
reached.  The  exceptions  are  the  work  by  Lee,  Mallett,  and  Yang  [lO]  for 
the  plane-strain  extrusion  with  frictionless  curved  dies  using  the  elastic- 
plastic  finite-element  method  and  the  rigid-plastic  analysis  of  axisymmetric 
extrusion  through  frictionless  conical  dies  by  Shah  and  Kobayashi  [6].  In 
the  latter  work,  the  authors  investigated  a possibility  of  applying  the 
matrix  method  to  steady-state  metalworking  processes  and  demonstrated  that 
the  analysis  of  a steady-state  extrusion  can  be  made  by  the  matrix  method. 

The  analysis  of  extrusion  in  this  section  is  an  extension  of  this  work. 

A.  Computational  conditions  and  procedures 

Boundary  conditions  and  mesh  system  The  boundary  conditions  smd  the 
mesh  system  used  for  the  analysis  of  extrusion  through  conical  dies  are  shown 
in  Fig.  1.  The  material  in  the  container  moves  axially  with  the  uniform 
velocity  of  unit  magnitude.  The  container  is  assumed  to  be  frictionless, 
and  along  the  conical  die  surfaces,  the  tangential  traction,  which  is  equal 
to  the  frictional  stress  at  the  die-workpiece  interface,  is  prescribed. 

The  extruded  material  moves  axially  with  the  uniform  velocity  of  the  magni- 
tude determined  from  the  area  reduction  and  the  incompressibility  relation- 
ship. Also,  no  traction  acts  along  the  surfaces  of  the  extruded  part. 


Along  the  axis  of  symmetry,  the  conditions  are  that  the  shear  traction  and 
the  radial  velocity  must  vanish. 

It  must  he  noted  here  that  the  die  corners  were  slightly  modified  by 
connecting  the  two  material  nodal  points  located  closest  to  the  die  cornei*s. 
This  modification  w€is  made  in  order  to  avoid  high  singularity  of  the  velocity 
components  near  the  die  corners . 

B.  Work-hardening  materials 

In  the  analysis  of  a non-steady-state  process , the  effect  of  work- 
hardening can  be  readily  incorporated  into  the  analysis  by  computing  the 
incremental  strains  and  modifying  the  flow  stress  at  each  deformation  step  ac- 
cording to  the  increase  in  the  total  effective  strain.  In  the  analysis  of 
steady-state  processes,  however,  the  flow  stress  distribution  must  be  con- 
sistent with  the  final  effective  strain  distribution  according  to  the  mater- 
ial’s work-hardening  characteristics.  This  requirement  can  be  achieved  by 
using  the  following  computational  procedure.  During  the  iteration  process 
for  a converging  solution,  the  flow  lines  corresponding  to  the  latest 
velocity  field  are  constructed  after  each  iteration.  The  network  of  grid 
distortions  and  the  effective  strain  distributions  are  determined  from 
these  flow  lines.  The  effective  strains  for  all  elements  are  then  inter- 
polated from  these  values,  and  using  a given  stress-strain  relationship, 
corresponding  flow  stress  distributions  for  elements  are  determined.  Using 
this  new  flow  stress  distribution,  the  next  iteration  for  the  velocity  field 
is  carried  out  and  the  same  procedure  is  repeated  until  the  converged 
solution  for  the  velocity  field  is  obtained.  Since  the  solution  depends 
not  only  on  ^ut  also  on  the  flow  stress  distribution,  when  the 
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velocity  solution  has  converged,  the  flow  stress  distribution  and  effec- 


tive strain  distribution  also  match  each  other  according  to  the  stress- 
strain  behavior  of  the  materied.  Thus  the  convergence  of  velocity  solu- 
tions gives  the  correct  solution  for  the  work-hardening  material. 

C.  Extrusion  conditions 

The  extrusion  process  conditions  under  which  the  con5)utation  was 
carried  out  are  summarized  in  Table  1.  Since  the  final  solution  for  the 
extrusion  with  nonhardening  material  was  used  as  the  initial  guess  for 
work-h€urdening  materials,  the  results  for  non-work-hardening  materials 
are  8^.so  presented. 


Table  l(a):  Extrusion  process  conditions  for  computation, 
a (die  semi-cone  angle)  = 30® 

^ = 2.366,  area  reduction  1 - (5—)  = 0.82 

^0  ^i 
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Fig.  2 Stress-strain  curves  for  SAE  1112  steel 
(Yq  = T9,800  psi;  550.21  MN/m^)  and  for 

aluminum  alloy  202lt-T351  (Yq  * 1*7,778  psi 
329.1*2  MN/m^). 

SAE  1112  steel  The  stress  distributions  computed  differs  from  the 
actual  distrlbvttion  by  a hydrostatic  component.  This  hydrostatic  component 
was  determined  by  setting  the  net  axial  force  along  the  exit  boundary  eqvtal 


to  zero,  and  the  magnitude  of  stresses  were  corrected  accordingly.  The 
average  extrusion  pressures  were  then  determined  from  the  stresses  along 
the  entrance  boundary.  They  are  summarized  in  Table  2. 


Table  2:  Average  extrusion  pressure, 
a = 30® 


The  general  trend  of  the  velocity  distributions  is  the  same  for  all 
friction  conditions.  A typical  example  is  given  in  Fig.  3.  The  two  velocity 
components  are  plotted  as  functions  of  the  radial  coordinate  at  various 
locations  in  the  axial  direction  within  the  deformation  zone.  These 
distributions  are  in  general  agreement  with  the  measurements  in  the  visio- 
plasticity  study  of  axlsymmetrlc  extrxision  of  lead  [ll]. 

Detailed  differences  of  deformation  characteristics,  due  to  material 
properties  and  friction  at  the  die-workpiece  interface,  are  more  clearly 
indicated,  for  exanqile,  in  grid  distortions.  The  steady-state  grid  distor- 


tion patterns  are  compared  for  non-work-hardening  and  work-hardening  cases 
for  the  two  friction  conditions  in  Fig.  U.  With  reference  to  Fig.  U(a), 


Fig.  U Grid  distortion  patterns  for  non-workhardening  and  workhardening 
materials  with  (a)  die  frictional  stress  f » 0,  and  (b)  f = O.UY-. 


it  is  clearly  seen  that  the  grid  lines  that  are  originally  perpendicular 
to  the  axis  of  the  workpiece  distort  and  show  the  double  peak  in  the  ex- 
ti*uded  part  for  non-work-hardening  material.  However,  for  a work-hardening 
material,  this  tendency  lessens.  With  increasing  friction  at  the  die- 
workpiece  interface,  as  shown  in  Fig.  U(b),  the  double-peak  distortion 
almost  disappears  for  both  work-hardening  and  non-work-hardening  materials. 
Another  important  feature  of  grid  distortions  can  be  seen  from  the  vertical 
grid  distortion  at  the  exit,  as  plotted  in  Fig.  5-  The  distortion  is  greater 
for  a work-hardening  material  compared  to  that  for  a non-work-hardening 
material,  for  small  friction  at  the  interface,  but  the  trend  is  reversed 
for  large  interface  friction. 

The  effective  strain-rate  distributions  are  shown  in  Fig.  6.  The 
figure  also  shows  the  boundaries  between  the  rigid  and  plastically  deforming 
regions  by  a near-zero  strain-rate  contour.  The  rigid  plastic  boundaries 
at  the  exit  are  almost  identical  in  shape  and  location  for  all  cases.  How- 
ever, the  rigid-plastic  boundaries  at  the  entrance  moves  backwards,  indi- 
cating larger  deforming  zones,  with  increasing  interface  friction.  This 
observation  applies  to  both  work-heu'dening  and  non-work-hardening  materials. 
Also,  the  results  show  that  the  deformation  zone  size  becomes  larger  for 
work-hardening.  The  difference  is  more  pronounced  with  increasing  friction. 
Neglecting  small  details,  the  effective  strain-rate  distribution  is  identicsil 
for  both  materials . The  strain-rate  increases  gradually  from  the  entrance 
tow€u:d  the  exit  and  near  the  exit  there  is  a sharp  drop.  As  might  be  ex- 
pected, there  is  some  degree  of  strain-rate  concentration  near  the  die  corner. 

The  effective  strain-rates  are  integrated  along  the  flow  lines  to  yield 
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the  effective  strain  distributions.  A typical  effective  strain  distribu- 
tion is  given  in  Fig.  7.  As  for  vertical  sections, the  strain  is  the  high- 
est near  the  die  and  the  lowest  near  the  extrusion  axis  throughout  the 
deformation  zone.  The  effective  strain  distribution  at  the  exit  section 
is  also  shown  in  the  figiu*e,  which  indicates  that  the  strain  is  the  lowest 
near  the  extrusion  axis  and  increases  toward  the  periphery.  The  strain 
distributions  in  the  final  product  of  the  two  materials  are  plotted  for 
several  friction  values  in  Fig.  8.  For  both  materials  nonuniformity  of 
deformation  increases  with  increasing  die-workpiece  interface  friction. 

The  degree  of  nonuniformity  is  greater  for  small  friction  and  less  for  large 
friction  in  work-hardening  materials,  while  the  opposite  holds  true  for  non- 
work-hardening materials.  These  results  reflect  exactly  the  findings  in 
the  vertical  grid  line  distortion  discussed  with  reference  to  Fig.  5. 

The  distributions  of  the  hydrostatic  pressure  shown  for 

the  two  friction  conditions  in  Fig.  9.  The  distribution  patterns  are  the 
same  for  both  friction  conditions  and  for  both  materials  with  and  without 
work-heu”dening.  Note  that  the  hydrostatic  pressure  increases  its  magnitude 
throughout  the  deformation  zone  with  increasing  friction  at  the  die- 
workpiece  interface.  It  can  be  observed  also  that  along  the  die-workpiece 
interface,  the  hydrostatic  pressure  is  largest  at  the  entrance  and  shows  a 
minimum  at  some  distance  from  the  exit.  This  particular  feature  is  seen 
in  the  normal  pressure  distributions  along  the  die  interface  shown  in 
Fig.  10.  The  pressure  is  the  highest  at  the  entrance  and  it  decreases 
toward  the  exit,  then  increasing  again  near  the  exit.  With  Increasing 
friction,  the  pressure  Increases  but  the  distribution  pattern  remains  the 
same.  The  trend  of  the  curves  for  nonhardening  and  work-hardening  materials 
is  also  the  same  with  difference  in  magnitude  of  a constant  amount. 
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In  other  results  obtained  by  the  finite-element  analysis,  it  is  not 


only  interesting  to  find  the  directions  of  principal  stresses  in  the  deforma- 
tion zone,  but  surprising  to  find  them  to  be  almost  identical,  reg€u:dless 
of  interface  friction,  except  near  the  die,  in  spite  of  definite  deviations 
of  some  quantities.  A typical  result  is  given  in  Fig.  11  for  SAE  1112 
steel  with  f = O.UYq.  The  lines  indicate  the  direction  of  the  larger 
principal  stress  in  the  meridian  plane. 

The  results  shown  in  the  foregoing  were  selected  from  among  informa- 

f 

tion  obtained  from  the  computation.  Emphasis  in  presentation  was  placed 
mainly  on  the  effect  of  die-workpiece  interface  friction  on  the  detailed  J 

deformation  mechanics  and,  to  some  extent,  on  the  effect  of  materials 
properties . 

Aliiminum  alloy  202U-T3$1  The  results  for  this  material  are  also 
utilized  for  the  workability  study  in  Part  II  of  this  report  with  regard 
to  center  bursting.  The  major  variables  are  the  reduction  in  area  and 
the  friction  condition  at  the  die-workpiece  interface.  The  average  extru- 
sion presstires  for  various  reductions  and  for  the  two  friction  conditions 
are  summarized  in  Table  3- 


Table  3:  Average  extrusion  pressure, 
a = 45° 


Friction 

1.25 

1.6 

1.8 

2.0 

f = 0 

2.102 

2.913 

3.289 

3.567 

f = 0.4y- 

2.310 

3.4o8 

3.826 

4.202 

i 

When  the  average  extrusion  pressure  was  plotted  as  a function  of  ln(— ), 
the  relationships  were  linear  for  both  friction  conditions,  as  shown  in 
Fig.  12.  The  figure  shows  that  the  slope  of  the  straight  line  increases 
with  increasing  friction  and  that  the  inclusion  of  work-hardening  char- 
acteristics of  a material  increases  the  slope. 

The  velocity  distribution  for  al  alloy  202l*-T351  do  not  differ  from 
those  for  SAE  1112  steel  in  general.  However,  direct  comparison  of  the 
computed  velocity  fields  with  those  measured  in  the  visioplasticity  study 
[ll]  is  now  possible,  and  such  a comparison  is  shown  in  Fig.  13.  The  com- 
puted results  are  given  for  the  case  of  a frictionless  die,  since  the 
general  picture  of  the  velocity  distribution  was  altered  little  by  friction 
at  the  die-workpiece  interface.  The  comparison  of  the  computed  and  measured 
velocity  distributions  in  Fig.  13  reveal  that  the  theoretical  resialts  in 
both  velocity  components  are  in  agreement  with  the  experimental  results  in 
every  detail  of  distribution  characteristics.  This  is  one  of  the  convincing 
evidences  of  the  accuracy  of  the  rigid-plastic  finite-element  analysis. 
Althoxjgh  in  comparing  Fig.  13  with  Fig.  3 for  SAE  1112  steel,  some  minor 
differences  in  the  velocity  distribution  may  be  noted,  the  differences 
are  mainly  due  to  the  deformation  zone  characteristics  of  the  two  materials. 

Fig.  lU  shows  the  total  distortion  of  the  grid  line  perpendicular  to 
the  extrusion  ajcis  as  influenced  by  the  material  property  and  friction 
for  several  reductions.  Inhomogeneous  deformation  increases  with  increasing 
reduction  in  area  and  with  larger  friction  at  the  interface.  An  interesting 
feature  is  that  a double-peak  distortion  is  more  pronounced  with  increasing 
reduction  for  a non-work-hardening  material,  while  the  contrary  is  true  for 
a work-hardening  material. 
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The  effective  strain-rate  distributions  for  various  area  reductions 


with  two  friction  conditions,  are  shown  in  Fig.  15.  With  regard  to  the 
effect  of  friction  on  the  deformation  zone,  the  exit  boundary  remains  about 
the  same  but  the  deformation  zone  expands  as  friction  at  the  die-workpiece 
interface  increases.  This  concliision  is  the  same  as  that  drawn  for  SAE 
1112  steel. 

As  the  reduction  in  area  increases,  naturally  the  deformation  zone 
becomes  larger  and  the  magnitude  of  strain- rate  increases.  Common  features 
of  the  strain- rate  distributions  for  all  the  reductions  investigated  are: 

(l)  strain-rate  concentrations  occur  near  the  corners  at  the  entrance  and 
at  the  exit,  and  (2)  along  the  axis  of  symmetry  the  strain-rate  peak  appears 
in  the  middle  of  the  deformation  zone.  The  effective  strain  distributions 
across  the  extruded  bar  are  plotted  in  Fig.  l6.  The  strain  is  largest  at 
the  STirface  and  smallest  at  the  center.  An  interesting  result  is  that 
the  degree  of  nonuniformity  in  strain  (difference  between  the  largest  and 
smallest  strains)  is  greatest  for  the  smallest  area  reduction.  Also,  the 
effect  of  friction  on  the  final  strain  distribution  is  practically  none, 
except  that  the  magnitude  of  the  surface  strain  increases  slightly  with 
increasing  interface  friction. 

The  distribution  of  the  mean  stress  is  similar  to  the  pattern  obtained 
for  SAE  1112  steel.  As  the  reduction  decreases,  the  mean  stress  increases 
and  becomes  tensile  in  the  zone  near  the  center,  as  shown  in  Fig.  17. 

Another  finding  is  seen,  with  respect  to  the  die  pressure,  in  Fig.  l8. 

The  die  pressure  is  highest  for  the  smallest  reduction.  This  implies  that 
the  pressure  distributions  for  smaller  reductions  are  critical  for  die 
design.  The  distribution  pattern  is  the  same  as  that  for  SAE  1112  steel 
showing  a minimum  at  some  distance  from  the  die  exit. 
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Finally,  the  directions  of  the  largest  principal  stresses  are  plotted 
for  two  reductions  in  Fig.  19.  Again,  the  effect  of  die-workpiece  inter- 


face friction  is  negligible  and  the  pattern  appears  to  be  determined 
solely  by  the  geometriced  constraints. 


The  drawing  process,  which  is  characterized  by  smaller  die  angles  and 


smaller  area  reductions,  when  compared  with  those  for  extrusion,  has  not 
yet  been  treated  by  the  finite-element  method.  Although  it  appears  that 
the  differences  between  extrusion  and  drawing  are  merely  geometriced.  and 
in  the  magnitude  of  hydrostatic  components  (extrusion  being  in  comp-ession 
and  drawing  in  tension  predominantly),  it  is  not  known  if  the  results  ob- 
tained in  extrusion  can  be  extrapolated  to  those  in  drawing  according  to 
the  geometrical  conditions  and  a simple  concept  of  pushing  and  pulling. 
Therefore,  we  discuss  the  analytical  results  in  drawing,  emphasizing  a 
comparison  with  t hose  obtained  for  extrusion  in  the  previous  section. 

A.  Computational  conditions  and  procedures 

Pig.  20  shows  the  mesh  system  and  the  boiindary  conditions  used  for 
the  analysis  of  drawing.  The  cylindrical  surface  of  the  entering  bar  is 
traction- free.  The  die  land  is  provided  at  the  exit  of  the  die,  along 
which  the  frictional  stress  is  prescribed.  Other  botindary  conditions  are 
the  same  as  those  for  extrusion.  For  the  analysis  of  extrusion  the  shape 
of  the  die  comers  were  modified  in  order  to  avoid  singularity.  However, 
this  was  not  necessary  for  the  analysis  of  drawing  because  of  the  small 
die  angle. 

The  materieil  selected  for  the  computation  was  SAE  11^^  cold-drawn 
steel  whose  stress-strain  relationship  is  expressed  by 
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with  Yq  = 8U,800  psi 

= 58I+.68  MN/m^. 


This  material  was  selected  mainly  for  the  workability  study  in  drawing 
given  in  Part  2 of  this  report. 

The  drawing  process  conditions  for  which  the  computation  was  performed 
are  summarized  in  Table  U. 


Table  U:  Drawing  process  conditions. 
Material:  SAE  lll+lt  cold-drawn  steel 
Semi-die  angle  ct  = 6°  and  8° 


R^/Rq  =1.25  was  used  as  an  initial  guess  for  drawing  with  a ~ 6°  and 
R^/Rq  = 1.2  by  modifying  the  solution  according  to  geometrical  proportions. 
Thirty- four  (3^)  iterations  were  required  to  achieve  the  converged  solution 
with 


I This  solution  was  then  used  as  an  initial  guess  for  the  computation  of  all 

the  other  cases.  The  convergence  was  excellent  and  only  6 to  10  Iterations 
were  necessary  to  obtain  the  solutions  in  most  cases. 


I 
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B.  Results  and  discussion 

The  results  for  a = 6°  are  presented  and  discussed  here.  The  results 
for  a = 8®  are  presented  £uid  discussed  in  terms  of  workability  in  drawing 
in  Part  2 of  this  report. 

An  example  of  the  velocity  distribution  is  shown  for  = 1.2  and 

f = 0.25Yq  in  Fig.  21.  The  distribution  pattern  looks  somewhat  different 
from  those  for  extrusion.  This  is  because  of  small  reduction  and  small 
die  angle  and  beca\ise  the  toteil  distortion  is  small.  Fig.  22  shows  vertical 
grid  line  distortion  in  drawn  bars.  Greater  friction  produces  more  distortion. 
The  difference  in  the  total  distortion  due  to  die  friction  increases  with 
increasing  reductions. 

The  general  pattern  of  the  effective  strain-rate  distribution  is  the 
same  as  that  for  extrusion,  as  seen  in  Fig.  23.  Concentrations  of  strain 
rate  occur  at  the  die  corners.  Along  the  axis  of  symmetry,  strain-rate 
distribution  shows  a peak  at  the  midpoint  of  the  deformation  zone.  The 
deformation  zone  expeinds  with  larger  die-workpiece  interface  friction  as 
it  did  for  extrusion. 

The  variation  of  effective  strain  distribution  in  drawn  bars  is 
comparatively  small  due  to  a small  die  angle  and  small  reductions  (Fig.  2k). 
Fig.  2k  shows  also  that  the  effect  of  the  die  friction  on  the  strain  distribu- 
tion is  negligible.  By  comparing  Fig.  2k  with  the  results  of  extrusion 
in  Fig.  l6,  it  can  be  seen  that  the  die  angle  is  a most  important  variable 
in  controlling  nonuniformity  of  the  extruded  or  drawn  bar  properties. 

The  hydrostatic  pressure  distribution  is  shown  in  Fig.  25.  It  shows 
the  combined  characteristics  of  those  for  small  reductions  similar  to 


Fig.  17  and  those  for  relatively  extended  deformation  zone  similar  to 
Fig.  9*  The  magnitude  of  the  hydrostatic  component  increases  with  larger 
die  friction  in  drawing.  This  is  contrsu^r  to  the  results  with  extrusion 
shown  in  Fig.  9*  The  die  pressure  distribution  is  plotted  for  variotis 
reductions  and  for  two  die  friction  conditions  in  Fig.  26.  Higher  die 
pressure  is  obtained  for  smaller  reduction  as  it  was  for  extrusion.  How- 
ever, contrary  to  the  case  of  extrusion,  the  die  pressure  is  higher  with 
frictionless  dies  than  those  with  friction.  These  findings  on  die  pressure 
are  in  agreement  with  experimental  results  [12],  although  empirical  values 
are  average  die  pressures . 


Fig.  22  Total  distortion  of  a line  originally  perpendi 
cular  to  the  axis  of  drawing  for  variovis  area 
reductions. 


V.  SUMMARY  AND  CONCLUSIONS 

In  recent  years  the  matrix  method  has  been  applied  to  the  anedysis 
of  various  metalworking  processes , including  the  steady-state  extrusion 
process.  The  present  matrix  method  program  incorporates  the  treatment  of 
friction  at  the  die-workpiece  interface  and  an  improvement  of  computational 
efficiency  by  adopting  a new  convergence  criterion.  Using  this  program, 
the  detailed  mechanics  of  steady-state  extrusion  and  drawing  were  obtained 
under  various  process  conditions.  The  stress-strain  properties  of  SAE  1112 
steel  and  aluminum  alloy  202U-T351  were  -.ised  for  extrusion  and  the  stress- 
strain  curve  of  SAE  llUl*  cold-drawn  steel  was  used  for  drawing.  The 
variables  included  for  the  computation  were  the  die  angle,  die-workpiece 
interface  friction,  and  the  area  reduction.  The  computed  resiilts  were 
presented  in  terms  of  velocity  distribution,  grid  distortion,  strain-rate 
distribution,  and  stress  and  strain  distributions. 

Althov^i materials  properties  apparently  influence  metal  flow,  as  seen 
in  grid  distortion  patterns,  their  effects  on  the  overall  deformation  char- 
acteristics appear  to  be  insignificant  in  extrusion  and  drawing  processes. 
Among  other  variables,  friction  at  the  die-workpiece  interface  plays  an 
important  role  in  determining  the  detailed  mechanics  in  these  processes. 
With  increasing  friction  the  degree  of  grid  distortion  becomes  leirger  and 
the  deformation  zone  size  expands  in  both  processes.  But  the  effect  of 
friction  on  nonuniformity  of  product  properties  is  less  significant,  while 
the  die  angle  is  an  in5)ortant  variable  in  controlling  this  nonuniformity. 
The  contrast  between  the  extrusion  process  and  the  drawing  process  cam  be 
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found  in  the  die  pressure  distribution  and  the  distribution  of  hydrostatic 
stress  components.  In  both  processes , the  die  pressure  decreases  with 
increasing  reduction  in  area.  However,  the  die  pressure  is  greater  for 
larger  interface  friction  in  extrusion,  while  the  reverse  is  true  in 
drawing.  Similarly,  the  magnitude  of  the  hydrostatic  stress  component  is 
less  with  larger  interface  friction  in  extrusion,  but  increases  with  in- 
creasing friction  in  drawing. 

It  is  concluded  that  the  matrix  method  is  indeed  an  efficient  numerical 
method  which  provides  usefial  information  on  the  detailed  deformatio;v  char- 
acteristics for  various  process  variables.  In  order  to  ascertain  the 
accuracy  of  the  solutions  presented  here,  however,  an  extensive  theoretical 
and  experimental  investigation  is  still  needed. 


I . INTRODUCTION 


For  occtirrence  of  cracks  at  the  free  surface,  such  as  in  edge-cracking 
in  sheet  rolling  and  surface  cracks  in  upsetting,  the  fracture  criterion 
can  be  constructed  experimentally.  However,  for  predicting  internal 
fracturing,  formulations  of  fract^lre  criteria  under  general  deformation 
are  required.  In  this  part  of  the  present  investigation,  the  validity  of 
the  theory  on  ductile  fracture  developed  in  the  previous  report  [l]  was 
examined  by  the  experimental  data  found  in  the  literature.  Then,  combining 
the  formulation  of  fracture  criterion  with  the  deformation  mechanics  fo\ind 
in  Part  1,  the  workability  of  materials  in  extrusion  and  drawing  was  deter- 
mined. 

In  an  attempt  to  develop  a general  fracture  criterion,  a model  for 
void  growth  was  proposed  to  examine  the  stress  and  strain  fields  around 
voids.  The  solutions  were  obtained  by  the  finite-element  method.  The 
results  show  that  deformation  is  concentrated  along  the  narrow  band  in 
the  maximum  shear  stress  direction.  The  fracture  criterion  assumed  that 
voids  of  orders  of  meignitude  smaller  than  the  primary  voids  are  distributed 
throughout  and  that  frsu:turing  occurred  when  the  small  voids  grow  and 
touch  each  other  along  the  band  connecting  the  large  voids.  The  McClintock 
analysis  was  used  for  predicting  the  growth  and  coalescence  of  small  voids. 
Based  on  this  concept,  the  fracture  strain  formulation  by  McClintock  was 
modified.  Then  the  modified  form,  along  with  the  formulation  by  Cockroft 
and  Latham,  was  tested  by  the  experiments  on  surface  cracks  and  by  the 
tension  data  by  Bridgman.  Its  validity  was  further  examined  by  applying  the 
criterion  to  the  determination  of  workability  in  bar  extrusion  and  drawing. 


II . FRACTURE  CRITERIA 


1 


McClintock  and  his  coworkers  [13],  [1^+]  developed  solutions  for  void 
growth  under  the  transverse  stress  state.  Their  model  consists  of  a single 
elliptic  cylindrical  void  extending  in  one  direction  and  imbedded  in  a 
rigid  plastic  media.  The  major  and  minor  axes  of  the  void  coincide  with 
the  principal  stress  directions.  They  derived  an  approximate  solution 
for  the  rigid  work-hardening  materials,  a = Ke*^,  as 


where  R and  Rq  are  the  current  and  initial  mean  radii  of  the  hole,  respec- 
tively; and  a^,  principal  stress  components  along  the  major  and  minor 
axes,  respectively;  5 is  the  effective  stress  and  i is  the  effective  strain; 
m is  the  eccentricity  of  the  ellipse  defined  in  terms  of  the  semi-major 
and  semi -minor  axes  a and  b,  as 


m 


a - b 
a + b ’ 


(11) 


and  m^  is  the  initial  eccentricity.  They  assumed  that  the  cylindrical  hole 
is  contained  in  the  cell  and  that  the  material  is  composed  of  these  cells. 
Fracturing  was  assumed  to  occur  at  the  point  where  a growing  void  touches 
the  cell  boundfury  whose  deformation  is  assimied  to  be  the  s£une  as  that  of 
the  boundary  at  infinity.  The  fracture  strain,  neglecting  the  void  inter- 
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action,  is  then  expressed  (dropping  a transient  term)  by 


K ‘j 

f -1 
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¥ i 


log. 


e 2a, 


e,.  = 


{2(1  - n)  2 


’/3(1  - n) 


a + a, 

a b 


]]  * ^(^)} 


(12) 


where  and  a are  the  initial  values  of  hole  spacing  I and  a,  respectively. 
d>  U 8> 

The  above  derivation  was  given  for  generalized  plane  strain  as  an  approxi- 
mation for  three-dimensional  deformations.  In  three-dimensional  configura- 
tions, there  arejsix  modes  of  fractiore,  two  in  each  of  the  three  perpendi- 
cular pleines.  Whichever  one  of  these  six  modes  gives  the  smallest  fracture 
strain  is  the  actual  mode. 

The  model  considers  only  the  case  where  holes  form  at  zero  strain  by 
complete  separation  of  the  particle-matrix  interface.  The  review  of  the 
investigations  on  ductile  fracture  revealed  that  void  nucleation  is  a 
complex  process  and  that  the  nucleation  pattern  depends  on  particle  size, 
particle  composition,  and  possibly,  particle  distribution.  The  general 
observation  is  that  void  nucleation  occiurs  at  the  large  particles  first. 

If  a number  of  voids  form  and  grow,  then  these  voids  may  act  as  a stress 
raiser  in  the  matrix  and  cause  further  void  nucleation  at  other  particle 
sites.  Also,  there  are  many  indications  that,  during  deformation,  shear 
bands  between  large  voids  develop  and  that  voids  form  at  small  particles 
in  these  bands,  while  the  particles  outside  the  bands  remain  relatively 
inactive.  Thus,  the  formation  of  shear  bands  between  large  voids  and 
growth  and  coalescence  of  small  voids  in  the  shear  band  may  be  a mechanism 
of  eventual  fracture. 


J 
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Assuming  this  mechanism  for  fracture,  a model  for  the  growth  of  large 


voids  has  been  proposed.  Fig.  27  shows  an  assumed  distribution  of  large 


voids,  extending  in  arrays  to  infinity.  The  plane-strain  condition  is 


assumed  in  the  direction  perpendicular  to  the  (x^,x^)  plane.  The  minimum 


and  maximum  principal  stresses  act  along  the  bo\jndary  at  infinity  in  the 


directions  of  coordinate  axes  (o  ,0.  ).  The  solutions  of  this  void  growth 


model  were  obtained  by  the  elastic-plastic  finite-element  model  [l].  A 


result  of  the  computation  is  shown  in  Fig.  28.  As  Fig.  28  indicates. 


deformation  is  concentrated  along  the  narrow  oand  in  the  maximum  shear 


stress  direction,  and  the  concentration  of  strain  along  the  band  is  about 


2 /V/  3 times  the  overall  strain  for  most  of  the  band  and  about  h >>'  3 greater 


near  a large  void.  This  imp.lies  that  the  band  is  the  most  favorable  place 


for  small  voids  to  nucleate  and  grow. 


The  modified  form  of  fracture  criterion  was  obtained  by  assuming  that 


voids  of  orders  of  magnitude  smaller  than  the  primary  voids  are  distributed 


throughout  and  that  fracturing  occvirred  when  the  smaJ.1  voids  grow  eind  touch 


each  other  along  the  band  connecting  the  large  voids.  When  the  McClintock 


analysis  is  applied  locally  under  the  condition  prevailing  within  the 


shear  band,  we  obtain 


_a_ 

*^(1  - n)  ^ ^ 


o + a,  •>  a - a 


^(1  - n)  % ^ ^ !b 


1*  a 

where  K = — log  {- — ) and  f is  a factor  indicating  the  strain  concentration 
it  e 


along  the  shear  band. 
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When  the  stress  ratios  vary  during  deformation,  McClintock's  model 
leads  to  the  condition  given  by 


r r n)  , ^1 

0 ^»^(l  - n)  ^ o •'  a -* 


de  = K, 


(lU) 


instead  of  Eq.  (13),  where  K is  the  same  factor  as  that  in  Eq.  (13). 

Another  fracture  criterion  which  will  be  considered  in  the  present 
investigation  is  the  one  proposed  by  Cockroft  and  Latham  [15],  [l6].  They 
proposed  that  fracture  occurs  when 


-f 

re 


a(-*)  de  = C 


(15) 


for  a given  temperature  and  strain-rate,  where  o*  * highest  tensile  stress 
and  C = material  constant.  The  proposed  criterion  is  phenomenological 
and  implies  that  ductile  fracture  is  dependent  both  on  shear  strain  and 
on  tensile  stresses.  From  the  viewpoint  of  void  growth,  however,  it 
appears  more  reasonable  to  include  stress  ratios  in  the  formulation  of 
fracture  criterion.  Thus,  the  Cockroft-Latham  criterion  is  modified  to 


r (5*)  de  = C. 

'0  a 


(16) 


It  may  be  of  interest  to  note  that  Eq.  (lU)  reduces  to  Eq.  (l6)  with 
K = 2C,  if  the  argument  of  sinh  is  small  so  that  the  hyperbolic  sine 
function  is  approximated  by  a linear  function.  The  magnitude  of  the 
argtunent  of  sinh  depends  on  the  stress  state  as  well  eis  on  the  work-hardening 
coefficient.  When  n = 1,  Eqs.  (lU)  and  (l6)  are  identical.  For  other  vedues 

5U 


of  the  work-hardening  coefficient,  the  difference  between  the  two  criteria 


Oa  + 

is  small  if  < 1,  but  increases  with  increasing  values  of  , 

0 5 

as  shown  in  Fig.  29. 

We  now  examine  in  more  detail  the  two  criteria  expressed  by  Eq.  (lU) 
and  by  Eq . (16). 


A.  Surface  fracture 

Since  the  fractxire  condition  applicable  to  the  free-surface  cracks 
can  be  observed  experimentally,  several  investigators  [l7]-[20]  have  used 
upsetting  of  cylindrical  specimens  for  the  study  of  ductile  fractures.  It 
was  concluded  from  these  investigations  that  the  fracture  criterion  is  ex- 
pressed by 

= a.  - ^ (it) 

where  a is  a materieil  constant  and  and  are  the  circumferential  strain 
and  the  axial  strain  at  the  equatorial  free  sxurface,  respectively.  The  two 
principal  stresses  in  the  surface,  and  (o^  > Cg),  during  plastic 
deformation,  are  expressed  by 
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<j  2 + a 
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(18) 


2 

where  ot  = and  is  a hydrostatic  stress  acting  on  the  surface, 
incremental  effective  strain  di  is  given  by 

de  = — A.  + a + de.  , 


The 


(19) 


Under  these  conditions,  the  fracture  criteria  given  by  Eqs.  (lU)  and  (l6) 
are  both  path-dependent,  except  when  Eq.  (l6)  becomes  path-independent 
for  = 0.  Taking  linear  strain  paths  under  o^/o  = const.,  and  noting 
that  and  are  the  minimum  and  maximum  t;..resses,  respectively,  Eq.  (lU) 
gives,  for  non-work-hardening  materials  (n  = 0),  the  fracture  strain  as 


^ ^ f K ^ 
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for  a £ -2»  and 
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-t  sinh(^ 


for  “2  £ “ £ criterion  given  by  Eq.  (16)  becomes 
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from  which 


- — (-^)  = 

2A  + or7^  |(2  + a)  + —(A  + a + a^{^) 


Eqs.  (20),  (21),  and  (22),  with  0^/5  = 0 are  plotted  for  K = 0.5  and  C = 
0.25  in  Fig.  30. 

It  can  be  seen  in  Fig.  30  that  the  fracture  line  given  by  Eq.  (22) 
has  a slope  of  for  0^/5  = 0 (free-surface  cracks),  and  Eqs.  (20)  and 


Fig.  31  Comparison  of  two  fracture  criteria  for  various  values 
of  a /a. 
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(2l)  result  in  almost  the  same  fracture  strains  as  those  according  to 
Eq.  (22)  if  the  relationship  K = 2C  is  maintained.  Similar  plots  of  Eqs. 
(20),  (21),  and  (22)  for  various  values  of  oJa  are  given  in  Fig.  31. 

B.  Fracture  in  uniaxial  tension 

Bridgman  [21 ] investigated  plastic  flow  and  fractiure  with  special 
emphasis  on  the  effects  of  hydrostatic  pressure,  using  uniaxial  tension. 

The  two  criteria  given  by  Eqs.  (lU)  and  (I6)  are  examined  for  internal 
fracturing,  using  the  Bridgman  experimental  results.  In  uniaxial  tension 
of  a bar,  the  critical  site  for  fract\ire  is  at  the  axis  of  symmetry  in  the 
neck  section  where  e,  = £ and  a = Therefore,  0=0=  o and 

= 0^  = 0 + 0^.  Then,  Eq.  (l^t)  is  expressed  by 

-f 

re 

(F  + 1)  de  = K,  (23) 

where  F^  = h when  h > -1, 

F^  = -1  when  h ^ -1 

and 

h . § . 2 !n)l 

/3(1  - n)  a ^ 

assuming  that  the  negative  damage  rate  is  not  permitted. 

The  criterion  given  by  Eq.  (16)  becomes 


where 


F-  « — for  — > -1 

a 5 

and 

F-  = -1  for  — <-l. 

0 

Eqs.  (23)  and  (2U)  are  con5>ared  with  experiments  in  Figs.  32  and  33, 
respectively.  The  experimental  results  by  Bridgman  used  for  comparison 
are  summarized  in  Table  5. 

In  deducing  experimental  points  in  the  figures  from  the  data  in 

Table  5,  several assvmptions  and  approximations  eure  introduced:  (l)  the 

effective  stress  5 was  assumed  to  be  independent  of  the  hydrostatic  pressure, 

(2)  the  stress  0 at  the  neck  was  estimated  according  to  the  Bridgman  analysis, 
r 

and  (3)  the  effect  of  work-hardening  coefficients  was  neglected. 

The  comparison  shown  in  Figs.  32  and  33  reveals  that  the  fractiore 
criteria  are  both  reasonably  good  for  predicting  the  fracture  strains. 


Table  5 Bridgman  experimental  data 
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Table  5 (continued) 
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III.  WORKABILITY  IN  EXTRUSION 


For  workability  in  extrusion  with  reference  to  center  bursting, 
experimental  observations  are  limited  and  information  about  the  stresses 
developed  in  metal  being  deformed  by  extrusion  is  insufficient.  Thus, 
Latham  and  Cockroft  [15]  speculated  the  conditions  in  axisymmetrlc  extru- 
sion by  assuming  that  the  slip-line  theory  can  be  applied  even  though  this 
is  appropriate  only  to  plane-strain  conditions.  Also,  it  was  assumed  that 
lubrication  during  cold  extrusion  is  sufficiently  good  for  the  calculations 
to  be  based  on  zero  friction.  Based  on  the  tensile  plastic  work  density, 
they  determined  a damage  factor  for  a range  of  extrusion  conditions  for 
several  metals.  Avitzia*  [22],  [23]  Investigated  theoretically  center  burst 
ing  in  extrusion,  using  the  upper-bound  approach.  Because  the  fracture 
condition  was  not  incori>orated  in  the  approach,  some  objections  may  be 
redsed  regarding  the  validity  of  the  results.  However,  the  study  led  to 
defining  the  range  of  successful  extrusions  in  terms  of  die  angle,  friction 
conditions,  and  the  reduction. 

Roffmanner  [2U]  also  performed  slip-line  analysis  for  extrusion  and 
compared  with  vlsloplasticlty  results  for  axisymmetrlc  extrusion.  He 
found  that  the  calculated  meucimum  center  line  principal  stresses  were  in 
good  agreement  with  vlsioplasticity  results  when  the  pleuie  strain  and 
axisymmetrlc  results  were  compared  on  the  basis  of  equal  cross-sections. 
Combining  the  slip-line  analysis  for  deformation  and  the  fracture  criterion 
proposed  by  Cockroft  and  Latham  [l6],  HoffSanner  examined  workability  in 
extrusion.  In  Fig.  34,  the  three  parabolic  curves  represent  the  process- 
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Fig.  3^  Workability  criteria  for  center  burst  based 
on  a maximum  tensile  stress-strain  energy 
criterion  [23]. 
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required  strain  energy  as  a function  of  the  reduction  ratio  for  three  die 
angles.  Whenever  the  process-required  energy  exceeds  the  energy  available 
for  successful  deformation  of  the  material,  center  burst  should  be  observed. 
The  extrusion  observations  are  superimposed  on  these  ciurves  and  demonstrate 
the  near-perfect  agreement  between  practice  and  prediction. 

A.  Experimental,  data 

Pepe  [25]  conducted  a detailed  investigation  of  the  phenomenological 
development  of  center  burst  defects  dxaring  hydrostatic  extrusion.  Central 
bursts  were  developed  in  IO8O  pearlitic  steel  during  a 25-percent  reduction 
in  area  pass  throiigh  conical  dies  having  an  included  angle  of  U5°.  The 
material  had  been  previously  cold  reduced  by  hydrostatic  extrusion.  Central 
bursts  were  formed  diiring  the  second  path.  The  study  revealed  that  central 
bursts  develop  from  regions  of  structural  damage  and  grow  to  completion  by 
a two-stage  fracture  process. 

Zimmerman  et  al.  [26]  tested  the  occurrence  of  center  bursting  using 
I02U  steel.  Hot-rolled  102U  steel  bars  were  drawn  and  then  extruded  in 
three  steps.  The  final  extrusion  step  was  utilized  to  examine,  through 
a variety  of  die  semi-angles,  whether  central  bursting  was  produced  or  not. 
They  found  that  among  UOOO  (lOOO  shafts  from  each  of  the  four  heats)  shafts, 
center  bursting  was  detected  in  about  shafts  from  one  of  the  heats. 

The  above  observations  of  center  biirsting  in  extrusion,  as  some  other 
experimental  studies  show,  were  all  made  in  material  which  has  been  pre- 
deformed. Although  the  present  theoretical  method  is  applicable  to  the 
determination  of  stress  and  strain  distributions  in  extrusion  for  the 
materiel,  which  heis  distributed  strength  and  is  capable  of  following  the 


i 

I previous  deformation  history,  the  experimental  observation  of  center 

btirsting  during  a single  pass  of  deformation  has  been  sought  because  the 
process  condition  in  the  computation  can  then  be  specified  more  faithfully 
according  to  the  experimental  procedure.  In  the  present  investigation, 
experiments  reported  by  Hofflnanner  [2U]  were  used  for  examining  the 
validity  of  the  workability  theory. 

Hoffmanner  observed  center  bursting  in  a single-step  extrusion  for 
aluminum  alloy  202U-T351  at  room  ten5)erature.  The  experimenteuL  conditions 
and  the  results  were  as  follows: 

1.  Material : The  specimens  were  obtained  from  a 3-inch  diameter  of 
aluminum  alloy  202^*  in  as-received  T351  condition. 

2.  Tensile  test : Testing  was  performed  at  room  temperature  on  an 
Instron  testing  machine  at  a nominal  strain  rate  of  0.1  in/in/min. 
Measurements  of  extension  were  performed  continuously  until  fracture 
with  an  extensometer  exhibiting  a strain  sensitivity  of  at  least 

I 0.0005  in/in.  Measurements  of  the  minimum  cross-sectional  radius 

j 

(a)  and  neck  radius  of  curvature  (R)  were  performed  in  two  directions 

at  90®  to  each  other  either  continuously  by  photographing  the 

I 

I specimen  or  discontinuously  by  removing  the  specimen  from  the  test 

j fixt\xres  and  performing  these  measurements  on  an  optical  comparator. 

3.  Compression  test:  Con^jression  testing  was  investigated  by  using 
simple  upsetting  of  cylinders.  For  strain  measurements,  an  ortho- 

) gonal  array  of  Vickers  diamond  pyramid  impressions  was  accurately 

I 

I 

j placed  at  0.050-inch  separations  about  the  exact  center  of  the 

^ specimens  along  subsequent  directions  of  principsil  normal  stress. 

These  directions  corresponded  to  the  direction  of  loading  and 


the  direction  normal  to  it  through  the  center  of  the  specimen. 
Eight  impressions  were  placed  in  each  of  the  two  directions. 
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U.  Extrvision;  Extrusion  was  performed  at  a ram  speed  of  12  in/min 
with  dies  of  (60°),  90°,  and  (120°)  Included  angles  with  anhydro\is 
lanolin  as  a lubricant.  Both  the  billets  sind  the  container  were 
lubricated  prior  to  extrusion  and  load-time  curves  were  recorded 
on  an  oscillograph  during  deformation.  The  billets  were  machined 
with  one  end  contoured  to  match  the  die  and  approximately  0.005 
inch  under  the  container  diameter.  A surface  finish  in  the  20  to 
30  microinch  rms  range  was  specified.  Initial  billet  diameters  of 
1.780,  I.U12,  and  0.812  were  used.  The  two  larger  diameters  were 
used  in  extrusion  through  90°  dies  to  study  single-  and  multiple- 
pass  extrusion  and  the  0.8l2  diameter  billets  were  used  in  single- 
pass reductions  through  all  three  die  angles. 

i 

S 

Resiilts  ^ 

j 

1.  Tension  test  re  splits  j 


Table  6 


Specimen 

number 

Gauge 

section 

diameter 

(inch) 

Initial 

a/R 

ratio 

Final 

a/R 

ratio 

Fracture 

strain 

Remarks 

TNLHl 

0.250 

0 

Longitudinal  specimen 

TNLH2 

0.250 

0 

0.342 

0.354 

Longitudinal  specimen 

TNTHl 

0.250 

0 

0.ll*6 

Transverse  specimen 

TNTH2 

0.250 

0 

0 

0.156 

Trsuisverse  specimen 
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2.  Compression  test  results 


Table  7 


Initial  Final 

Initial 

height/diameter 

ratio 

Fracture 

strain 

0.50  1.61 

I.1I+5 

0.265 

0.50 

2.2k 

- 

0.50  1.58 

I.II+5 

0.281 

Remarks 


Teflon  lub. , 
sheared  over 


3.  Extrusion  results 


Table  8 


Specimen  no. 


HN-IOIA 
HN-IOIB 
HN-IOIC 
1+51 
U52 
HN-990 
HN-999 
HN-9OB 


Included 
die  angle 

Remarks  ; 

1.2lt5 

90° 

Unlubricated,  billet  stuck 

1.326 

90° 

Centerburst 

1.2ii9 

90° 

Centerburst 

1.157 

90° 

Centerburst 

2.328 

90° 

i 

1.988 

90° 

Centerburst  j 

1.089 

90° 

1 

2.132 

90° 

Comnutation.  results,  and  discussion 


Using  the  matrix  method,  the  extrusion  process  was  analyzed  for  the 
material,  al  alloy  2021+-T351,  and  the  detailed  mechanics  were  presented 
in  Part  1 of  this  report.  The  process  conditions  were:  semi-cone  angle 
a = 1*5°,  area  reduction  ^'^5,  1.6,  1.8,  and  2.0  for  the  two  friction 

conditions,  f = 0 and  0.1+Yq.  The  critical  site  for  occurrence  of  center 
bursting  is  along  the  axis  of  extrusion.  The  stress  and  strain  states 
along  the  axis  of  extrusion  resulted  from  the  computation  described  in 
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to  that  in  uniaxial  tension,  the  queintities  necessary  for  workability  are 
[ (F^  + l)  di  from  Eq.  (23)  for  one  criterion 

and 

[ (Fp  + l)  di  from  Eq.  (.2k)  for  another. 

•'0 

These  quantities  are  plotted  as  functions  of  area  reduction  in  Figs.  35  and 
36.  From  the  figures,  some  conclusions  immediately  follow.  The  two  quanti- 
ties are  almost  identical,  with  a factor  of  2.  The  magnitudes  decreased 
greatly  when  work-hardening  of  materials  was  included.  The  magnitudes 
decreased  further  with  friction  at  the  die-workpiece  interface.  Note  that 
the  curve  for  non-work-hardening  materials  in  Fig.  36  is  in  good  agreement 
with  the  Hoffmanner  result  given  in  Fig.  3^  for  90°  included  angle.  This 
is  somewhat  surprising  because  the  curve  by  Hoffmanner  was  obtained  by 
approximations  based  on  the  slip-line  analysis. 

The  critical  values  of  these  quantities  (K  in  Eq.  (23)  and  C in  Eq. 
(2U))  above  which  fracturing  is  predicted  can  be  obtained  from  the  fracture 
data  in  tension  or  compression  tests.  Hoffmanner  estimated  C = 0.3  (K  = 2C) 
from  the  tension  test  results.  On  the  other  hand,  another  estimate  of 
C = 0.08  was  obtained  from  the  coj.pression  test  results.  If  we  apply  the 
critical  value  of  C = 0.3  to  the  results  in  Fig.  36,  fracturing  should 
not  occiir  in  all  the  reductions.  If  we  assme  C = 0.08  to  be  a critical 
value,  then  center  bursting  would  be  observed  over  the  range  1.08  < 

< 1.86  with  frictionless  die  and  l.l6  < Rj^/Rq  ^ 1.60  with  friction  f = O.I+Yq 
Although  the  prediction  of  workability  with  a critical  value,  C = 0.08, 
appears  to  be  in  good  agreement  with  the  experimental  res\ilts  in  Table  8, 
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Fig.  35  Variations  of  (F^  + l)  de  as  a function  of  area  reduc 

•'0 

tion  in  extrusion.  Materials:  non-work-hardening  and 
al  eLl.loy  202Jt-T351. 


Fig.  36  Variations  of  I (Fp  + l)  di  as  a fiinction  of  area  reduc- 

■'0 

tion  in, extrxision.  Materials:  non-work-hardening  and 
al  alloy  202U-T351. 
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the  proof  of  the  validity  of  the  fracture  criteria  is  not  conclusive. 

First,  it  is  not  certain  which  of  the  critical  values,  C (or  K),  is  proper 
to  use  for  workability  predictions.  Evidently,  the  fracture  behavior  is 
not  homogeneous . This  is  due  to  the  complex  microstructure  of  the  material 
and  the  role  of  microstructure  in  determining  the  fracture  strains  is  not 
yet  understood.  Second,  the  observations  of  center  bursting  in  Table  8 
revealed,  as  Hoffmanner  reported,  that  it  always  occurred  by  complete 
separation.  In  spite  of  the  fact  that  the  fracture  surfaces  always  pos- 
sessed the  centrally  located  conical  cavity  typical  of  the  center  burst, 
whether  fracturing  initiated  at  the  center  is  in  doubt. 
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IV.  WORKABILITY  IN  DRAWING 


To  determine  workability  in  drawing,  the  sajne  approach  as  that  for 
extrusion  was  applied.  The  process  mechanics  were  analyzed  by  the  matrix 
method  and  the  stress  and  strain  distributions  along  the  ajcis  of  drawing 
were  computed  for  various  reductions  of  area.  The  tension  and  compression 
tests  were  performed  to  obtain  the  stress  and  strain  property  of  the 
material  and  to  determine  the  criticeil  value  for  fracturing.  Combining 
these  data,  the  workability  in  drawing  was  predicted.  The  drawing  experi- 
ments were  conducted  at  room  temperature  to  examine  the  predictions.  The 
material  was  SAE  llUl*  cold-drawn  steel.  The  stress-strain  property  of 
the  materieil  is  given  in  Part  1. 


A.  Process  mfechanics 

The  process  conditions  were:  semi-die  angle  a = 8°;  area  reduction 
R^/Rq  = 1.05,  1.10,  1.2,  1.3;  end  friction  f = 0,  0.25yQ.  The  stress  and 
strain  distributions  were  obtained  and  the  significant  quaintities  for 
fracturing  were  determined  for  each  area  reduction. 

Because  the  two  fracture  criteria  discussed  in  the  previous  sections 
resulted  in  almost  identical  predictions  of  fracture  conditions,  only  the 
quEintities  (normalized  energy  density)  defined  by  (Fp  + l)  di  was  shown 

Jo 


as  a function  of  the  area  reduction  in  Fig.  37.  In  calculating  the  energy 
density,  the  stress  components  were  determined  by  correcting  the  Hydro- 
static part  of  the  stress  components  from  the  total  force  balEuice.  It  is 
seen  in  Fig.  37  that  the  energy  densities  for  the  critical  path  along  the 
axis  of  drawing  increased  monotoniceilly  with  increasing  area  reduction. 
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1.05 


1.2 


1.3 


Rj/Ro 
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Pig.  37  Variations  of  (F_  + l)  de  as  a function 

■'o 

of  area  reduction  in  drawing.  Material: 
SAE  llUU  cold-drawn  steel. 


They  were  larger  for  larger  dle-vorkpiece  interface  friction.  This  is 
contrary  to  the  extrusion  results  shown  in  Figs.  35  and  36.  To  determine 
workability,  this  energy  density  level  must  be  compared  with  the  critical 
value  C for  fracturing,  estimated  from  the  fracture  data  in  tension  and 
compression  tests. 


B.  Tension  and  compression  tests 


Tension  specimens  were  machined  longitudinally  from  a 0.5-in.  diameter 
bar.  Testing  was  performed  on  an  Instron  machine  with  a plate  speed  of 
0.02  in/min.  The  load-displacement  curve  was  recorded  and  the  instantaneous 
neck  radii,  as  well  as  neck  diameter,  was  measured  until  fracture.  The 
critical  energy  density  was  estimated  from  the  tension  test  as  C(t)  = 0.271. 
The  Bridgman  analysis  was  used  for  an  approximate  stress  analysis  during 
neck  formation. 


Compression  specimens  were  machined  from  the  same  bar  of  0.5-in. 

diameter.  The  height  and  the  diameter  of  the  specimens  were  0.50  in.  eind 

0.1+9  in.,  respectively.  The  grid  lines  were  printed  on  the  cylindrical 

surface  of  the  specimens,  and  the  strains  at  fracture  were  determined  from 

the  grid  distortion  of  the  deformed  specimens.  The  critical  value  C was 

estimated  from  the  fracture  strains  according  to  Eq.  (22)  and  was  given  by 

C/  V = 0.093.  It  is  to  be  noted  that  the  fracture  behavior  is  different 
(c) 

in  tension  and  in  compression. 


C.  Workability 

The  energy  density  in  drawing  shown  in  Fig.  37  is  compared  with  the 
critical  value  for  fracturing  in  order  to  predict  workability  of  SAE  lll+l* 
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steel  in  drawing.  It  is  seen  again  that  the  prediction  is  drastically 
different  depending  upon  which  critical  value  of  tension  and  compression 
is  applied.  With  the  critical  value  of  0.09  which  was  estimated  from  the 
compression  tests,  center  bursting  would  occur  in  all  drawing  with  the 
area  reduction  larger  than  = 1-03,  with  very  little  effect  of  die- 

workpiece  interface.  On  the  other  hand,  if  is  taken  as  the  critical 

value,  the  prediction  is  that  for  frictionless  dies  no  center  curst  would 
occur  up  to  = 1-36,  and  that  this  limiting  value  decreases  with 

increasing  die  friction;  for  example,  for  friction  f = 0.25Yq, center  burst 

would  not  be  observed  for  R./R  less  than  l.l6. 

1 o 

In  order  to  examine  the  workability  predictions , bar  drawing  experi- 
ments were  performed.  The  bar  specimen  configuration  is  shown  in  Fig.  38. 
A threaded  grip  was  used  for  drawing.  This  scheme,  however,  sets  the  draw 
ing  limit  according  to  the  maximum  drawing  force  which  can  be  sustained 
by  the  threaded  part  of  the  specimen.  Commercially  available  carbide 
drawing  dies  with  8°  semi-included  die  angle  were  used.  The  steady-state 
drawing  force  was  measured,  sind  the  occurrence  of  center  bursting  was 
checked  with  drawn  bars.  The  experimental,  conditions  and  results  are 
summarized  in  Table  9. 

Table  9 Experimental  conditions  and  results. 


Center 

bursting 


Speci- 

men 

no. 

(Outlet 
diam. ) 

2Ro 

(Inlet 
diam. ) 
2R, 

(Reduc . 

ratio) 

Ri/Eo 

Lubrication 

Ave . drawing 
pressure 

(L/UE^Y„) 

k 

0.388" 

0.1*85" 

1.25 

White  lead 

0.907 

in  oil 

8 

0.388" 

0.1*65" 

1.198 

rr 

0.713 

5 

O.U195" 

0.1*85" 

1.156 

ft 

0.691 

l6 

0.1*195" 

0.1*65" 

1.108 

If 

0.361 

1 

0.1*51" 

0.1*85" 

1.075 

II 

0.1*35 

9 

0.1*51" 

O.U65" 

1.031 

II 

0.1*28 

21 

0.1*51" 

0.5" 

1.109 

Dry 

0.661* 

7 

0.1*65" 

O.U65" 

1.198 

Dry 

0.997 

In  Fig.  38  the  experimental  and  theoretical  average  drawing  stresses  were 


plotted  as  fxmctions  of  the  reduction  ratio  R^/Rq-  Theoretical  curves 
were  for  frictionless  dies  (f  = O)  and  for  the  frictional  stress  f = 0.25Yq. 
The  experimental  values  with  and  without  lubrication  fell  between  these 
two  theoretical  curves.  The  drawing  limit  due  to  the  present  experimental 
scheme  is  also  shown  in  Fig.  38.  Under  the  lubricated  condition,  drawing 
over  the  reduction  range  of  R^/Rq  = 1.031  1.25  produced  no  center  burst- 

ing. This  suggests  that  the  vise  of  the  critical  value  C,  estimated 
from  the  compression  fracture  data  for  workability  prediction,  was  not 
appropriate.  With  regard  to  the  criticeLL  value  of  = 0.271,  it  was 

fovind  that  the  drawing  limit  set  by  the  present  experimental  scheme  was 
very  close  to  the  critical  value  as  shown  in  Fig.  37.  In  fact, 

drawing  was  not  possible  for  three  cases  of  drawing  under  the  dry  condition. 
As  a result,  it  was  not  possible  to  determine  the  applicability  of  as 

a fracture  criterion  to  the  occvirrence  of  center  bursting  in  bar  drawing. 


V.  SUMMARY  AND  DISCUSSION 

For  workability  in  metalworking  processes,  formulation  of  a ductile 
fracture  criterion  must  be  simple  (or  in  the  form  of  easy  applicability) 
and  yet  reflect  physical  mechanisms  of  fracture  so  that  the  predicted 
fractxire  behavior  is  close  to  the  one  of  real  materials  within  reasonable 
accuracy.  With  a two-size  voids  model  a modified  interpretation  to  the 
McClintock  formulation  of  fracture  strains  was  given.  It  was  shown  then 
that  this  formulation  and  the  one  by  Cockroft  and  Latham  resulted  in 
almost  identical  fracture  strains  for  occurrence  of  free-surface  cracks 
and  for  fracturing  in  uniaxial  tension.  It  was  demonstrated  further  that 
these  formulations  predicted  reasonably  well  the  experimentally  observed 
fracture  strains.  In  examining  the  validity  of  the  prediction  by  experi- 
ments, the  stress-strain  histories  at  a critical  site  were  determined  by 
xising  the  continuous  experimentad  observations  on  distortion  in  the  case 
of  surface  cracks  and  by  the  Bridgman  analysis  combined  with  experimented 
meeisurements  of  load  and  neck  geometry  in  the  case  of  uniaxial  tension. 

In  applying  the  criteria  to  occurrence  of  center  bursting  in  axisymmetric 
extrusion  and  drawing,  the  determination  of  stress  and  strain  paths  at  a 
critical  site  was  provided  theoretically  by  the  matrix  method  of  analysis. 
The  critical  value  of  the  material's  capability  against  fracturing  vas 
obtained  by  tension  and  compression  experiments.  Workability  in  extrusion 
was  examined  for  aluminum  alloy  202U-T351  using  the  data  found  in  the 
literature  and  the  experiments  of  workability  in  drawing  was  attempted 
for  SAE  llUU  cold-drawn  steel  to  test  the  predictions.  The  results  of 
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validation  were  inconclusive.  This  is  attributed  to  the  fact  that  the 
fracture  behavior  has  directionality  as  evidenced  from  the  tension  and 
compression  fracture  data.  Furthermore,  most  of  the  previous  investiga- 
tions on  center  bursting  in  extrusion  and  drawing  have  reported  observa- 
tions of  center  bursting  in  pre-deformed  materisds . Although  al  alloy 
2024-T351  was  reported  to  have  produced  center  bursting  in  a single 
extrusion  path  at  room  temperature,  it  occurred  always  by  complete  separa- 
tion. With  regard  to  bar  drawing  of  SAE  llltU  steel,  the  drawing  limit 
set  by  the  present  experimental  scheme  was  very  close  to  the  limiting 
energy  density  level,  thereby  making  it  extremely  difficult  to  perform 
drawing  near  the  critical  boundary  of  the  sound  and  defect  zones. 

In  spite  of  inconclusive  results  of  workability,  the  present  investi- 
gation revealed  several  significant  findings.  In  extrusion,  the  energy 
density  level  critical  for  center  bursting  was  less  for  work-hardening 
materials  when  compared  with  that  for  a non-work-hardening  material. 

This  energy  density  level  is  reduced  with  increasing  friction  at  the  die- 
workpiece  interface.  It  increases  first  with  increasing  area  reduction 
and  then  decreases  indicating  a maximvun  as  a fxmction  of  eirea  reduction. 

In  drawing,  the  energy  density  level  monotonically  increases  with  the 
area  reduction,  and  increases,  as  opposed  to  that  in  extrusion,  with 
increasing  interface  friction. 

Although  conclusive  validity  of  the  present  workability  theory  in 
extriision  and  drawing  awaits  more  extensive  and  systematic  experimental 
investigations,  as  well  as  theoretical  calciilations , the  method  of  computa- 
tion is  available  and  the  approach  has  been  cleared  toward  complete  under- 
standing of  ductile  fractvtre  in  metalworking  processes. 
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The  four  surrotmding  nodal  points  (I,  J,  K,  L)  of  an  element  are  arranged 
in  a clockwise  direction  with  the  lowest  index  assigned  to  I (e,g.,  for 
the  first  element  in  the  figure  (I,  J,  K,  L)  = (l,  9,  10,  2)). 

The  output  values  of  stress  and  die  force  are  all  normalized  by  the 
initial  yield  stress  and  are  positive  in  the  sense  of  positive  r or  z 
direction. 

In  addition  to  preparing  the  input  data  cards,  the  following  variables 


THETA 

NIPTS 

NJPTS 


m\ 


NTIMES: 

NMAX: 

RENTER: 


YEXIT: 


Semi-included  angle  of  conical  die 

Number  of  lines  in  mesh  system  parallel  to  the  axis  of  symmetry 

Number  of  lines  in  mesh  system  perpendicular  to  the  axis  of 
symmetry 

Number  of  flow  lines  to  be  constructed 

Total  number  of  increments  for  strain  integration 

Radius  of  the  billet  at  entrance  (it  is  convenient  to  choose 
RENTER  = 1,0) 

z-coordinate  of  the  billet  at  the  die  exit  (YEXIT  = O.O) 


A.  Sequence  of  input  data  card  preparation  for  the  main  program 

1.  Heading  card  (12A6) 

Columns  1-72:  Output  title  to  be  printed  with  output,  up  to  72 
characters 

2.  Columns  1-10:  The  initieilly  assigned  value  of  the  deceleration 

coefficient  for  the  first  iteration.  Note  that  in 
a subsequent  iteration,  this  value  will  be  changed 
* automatically  according  to  the  behavior  of  functions 

in  order  to  ensure  fast  convergence  (usually  a = 0.5) 

I 3.  Material  property  and  program  control  cards  (Ul5,  FIO.O) 

J 

i Colxanns  1-5:  Value  of  ITER,  total  number  of  iterations  assigned 

j Columns  6-10:  Value  of  ITCONT;  1,  for  non-workhardening  materials 

I 0,  for  workhardening  materials 


Columns  11-15:  Value  of  NPUNCH;  0,  if  only  velocity  distributions 

to  be  punched 

1,  punch  all  informations;  the  out- 
put punched  cards  are  arranged 
as  follows: 

a.  Nodal  point  velocity  distributions  (u  ,u  ) (8F10.6) 

X*  z 

b.  Total  effective  strain  (e)  distribution  (8F10.6) 

c.  Strain  rate  (e^,e^,eg ,e^^,e)  distributions  (8F10.6) 

d.  Stress  (a  ,a  ,0^,0  ,5,a  ) distributions  (8F10.6) 

r’  z’  6 rz’  ’ m 

e.  Boundary  nodal  point  forces  and  velocities 
(F^,F^,u^,u^)  (8F10.6) 

Columns  15-20:  Value  of  NPRINT;  1,  if  nodal  point  and  element  data 

are  to  be  printed 
0,  otherwise 

Columns  21-30:  Value  of  FLIMIT,  assigned  value  of  accuracy  desired 

(lAul/lul),  program  will  stop  if  this 
value  is  achieved  (recommend  value  is 
0.00008) 


8U 


' 


I 


U.  Geometry  and  traction  boundary  control  card  (U15) 

Columns  1-5:  Value  of  NUMNP,  total  number  of  nodal  points 

Columns  6-10:  Value  of  NUMEL,  total  number  of  elements 

Columns  11-15:  Value  of  NUMPC,  number  of  traction  boundary  condition 

cards  to  be  read  (see  Sequence  6,  below) 

Columns  16-20:  Value  of  NBF,  total  number  of  nodal  points  along  the 

die  surface  at  which  force  calculation 
is  required  (see  Sequence  6,  below) 

5.  Nodal  point  data  cards  (15.  F5.0,  2F10.0) 

The  numbering  system  of  nodal  points  for  the  boundary  conditions 

Columns  1-5:  Nodal  point  nmber 

Columns  6-10:  Code  for  this  nodal  point 

0.0  if  R-  and  z-forces  are  specified 

1.0  if  R- velocity  and  z-force  are  specified 

2.0  if  R- force  and  z- velocity  are  specified 

3.0  if  R-  and  z- velocity  are  specified 

5.0  for  special  boundary  condition  (nodal  points  on 
die  surface) 

Columns  11-20:  R-coordinate  of  the  nodal  point 
Columns  21-30:  z-coordinate  of  the  nodal  point 

6.  Force  calculation  nodal  points  cards  (1615) 

These  cards  provide  the  nodal  point  numbers  at  which  the  force 
calculations  are  desired.  Each  card  provides  I6  nodal  points  punched 
in  every  five  coliimns.  The  total  number  of  these  nodal  points  should 
be  the  same  as  NBF  specified  in  Sequence  1+. 


Traction  boundary  condition  cards  (215,  *+F10.0) 

These  cards  provide  prescribed  traction  stresses  along  boundaries. 

If  the  prescribed  stresses  are  all  zero  or  no  traction  stresses  are 
prescribed  (e.g.,  for  the  case  of  extrusion  (or  drawing)  with  friction- 
less dies),  these  cards  can  be  omitted  and  put  NUMPC  = 0 of  Columns 
11-15  of  Sequence  T^/T^ 


|r^  S(J)_ 


Nodal 

I 


Nodal 

J 


For  the  case  with  constant  frictional  stress  (in  terms  of  initial 
yield  stress),  these  cards  are  needed.  The  total  number  of  cards 
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shoiild  be  the  same  as  the  value  of  NUMPC.  The  shear  stress  euid  the 
normal  stress  are  positive  in  the  direction  shown  in  the  figure, 
which  is  opposite  to  the  output  positive  stress  and  die  force  directions. 

Columns  1-5:  Nodal  ijoint  number  for  point  I 

Columns  6-10:  Nodal  point  number  for  point  J 

Columns  11-20:  Prescribed  normal  stress  in  terms  of  initial  yield 
stress  at  point  I 

Columns  21-30:  Prescribed  nonn£Ll  stress  in  terms  of  initial  yield 
stress  at  point  J (if  normal  stresses  are  not 
prescribed,  put  0.0  into  these  two  values) 

Columns  31-UO:  Prescribed  shear  stress  in  terms  of  yield  stress  at 
point  I 

Columns  ^*1-50:  Prescribed  shear  stress  in  terms  of  yield  stress  at 

point  J (e.g.,  for  the  case  of  constant  friction  stress 
with  f = O.UYq,  put  ”-0.U"  into  these  two  values) 

8.  Flow  line  data  cards  (8F10.0) 

Specify  the  R-coordinates  of  every  flow  line  to  be  constructed  before 
entering  the  die.  The  total  number  of  these  points  should  be  the  same 
as  NTIMES  specified  in  the  program  PORT.  Each  card  is  to  be  punched 
8 entering  points  each  in  10  columns. 

9.  Element  data  cards  (515) 

For  each  element  these  cards  provide  the  four  surrounding  nodal  point 
numbers  (l,J,K,L)  in  a clockwise  direction 

Columns  1-5:  Element  number 

Columns  6-10:  Nodal  point  ntamber  of  point  I 

Columns  11-15:  Nodal  point  number  of  point  J 

Columns  16-20:  Nodal  point  number  of  point  K 

Colxanns  21-25:  Nodal  point  number  of  point  L 

10.  Input  velocity  field  cards  (8F10.0) 

These  cards  can  be  an  initial  guess  of  velocity  field  or  a new 
velocity  field  obtained  from  the  previous  iteration.  These  input 
cards  are  arranged  in  the  following  order  (from  the  nodal  point  to 
the  last  nodal  point) 

Columns  1-10:  R-velocity  of  the  1st  nodal  point 

Columns  11-20:  z-velocity  of  the  1st  nodal  point 

Columns  21-30:  R-velocity  of  the  2nd  nodeil  point 

Columns  31-^0:  z-velocity  of  the  2nd  nodal  point 

Each  card  provides  four  noded  point  velocities 


B.  Sequence  of  Input  data  cards  for  program  EXDRSUM 

The  output  stress  distributions  from  the  matrix  iteration  method  are 
different  from  a uniform  hydrostatic  stress  from  the  actual  solution  over 
€ill  elements  in  extrusion  or  drawing. 

This  program  calculates  the  average  extrusion  (or  drawing)  stress  by 
either  the  energy  baleuice  method  or  the  method  of  zero  toted  aucial  force 
at  the  entrance  boundary  (for  drawing)  or  at  the  exit  boundary  (for  extru- 
sion). After  the  converged  solution  is  obtained  from  the  main  program,  the 
output  punched  ceirds  are  then  used  in  this  progreun  to  cedculate  the  correct 
stresses  according  to  the  following  sequence  of  input  data  cards: 

1.  Problem  control  card  (15) 

Columns  1-5:  Value  of  lEXDR;  1,  for  the  extrusion  problem 

0,  for  the  drawing  problem 

2.  Geometry  control  card  (Ul5,  FIO.O) 

Columns  1-5:  Value  of  NUMNP,  same  as  in  the  main  progreun 
Columns  6-10:  Value  of  NUMEL,  seime  as  in  the  main  problem 
Columns  11-15:  Value  of  NIPTS,  same  as  in  the  main  problem 
Columns  16-20:  Value  of  NBF,  same  as  in  the  main  problem 
Columns  21-30:  Value  of  THETA,  semi-included  cone  die  angle 

3.  Initial  yield  stress  card  (FIO.O) 

Columns  1-10:  Vedue  of  YIELD,  initial  yield  stress  of  the 
material  considered 

1+.  Nodal  point  data  cards  (I5»  F5.0,  2F10.0) 

Same  as  Sequence  5 of  the  main  program 

5.  Element  data  cards  (515 

Same  as  Sequence  9 of  the  main  program 

6.  Input  strain-rate  distribution  cards 

From  output  of  the  main  program  (Sequence  3-c) 

7.  Input  stress  distribvttion  cards 

From  output  of  the  main  program  (Sequence  3-d) 

8.  Input  boundary  nodal  point  force  and  velocity  cards 
From  output  of  the  main  program  (Sequence  3-e) 
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vooGHAw  exosMti  r*muT«oc'rPUT«ouNCMi 


EXOPSUM 

2 

exoPsuM 

3 

exor sum 

4 

EXDR SUM 

8 

exorsum 

6 

EXORSUM 

T 

EXORSUM 

0 

EXOR SUM 

9 

EKDRSUM 

10 

EXDRSUM 

11 

E KOR  SUM 

12 

EXOPSUM 

13 

EXOPSUM 

14 

EXOPSUM 

to 

EXORSUM 

16 

EXOP  SUM 

17 

EXOPSUM 

10 

EXORSUM 

19 

EXOPSUM 

20 

EXOPSUM 

21 

EXORSUM 

22 

EXORSUM 

23 

EXORSUM 

24 

EXOPSUM 

28 

EXOPSUM 

26 

EXOR  SUM 

27 

EXOPSUM 

20 

EXORSUM 

29 

EXORSUM 

30 

EXORSUM 

31 

EXORSUM 

32 

EXDRSUM 

33 

EXOPSUM 

34 

EXOPSUM 

3S 

EXORSUM 

36 

EXOPSUM 

17 

EXOPSUM 

38 

EXORSUM 

39 

EXDRSUM 

40 

EXDRSUM 

41 

EXOPSUM 

42 

EXDRSUM 

43 

EXORSUM 

44 

EXORSUM 

48 

EXDRSUM 

46 

EXORSUM 

47 

EXDRSUM 

40 

EXDRSUM 

A9 

EXOPSUM 

80 

EXOPSUM 

SI 

EXOPSUM 

S2 

EXOPSUM 

S3 

EXORSUM 

84 

EXORSUM 

S5 

EXOPSUM 

S6 

EXDRSUM 

57 

EXDRSUM 

80 

EXORSUM 

09 

EXORSUM 

60 

EXDRSUM 

61 

EXORSUM 

42 

EXORSUM 

63 

EXDRSUM 

64 

EXDR  SUM 

65 

EKDRSUM 

66 

EXDRSUM 

67 

EXORSUM 

60 

EXORSUM 

69 

EXDRSUM 

70 

EXDRSUM 

7| 

EXDRSUM 

72 

EXDRSUM 

73 

EXDRSUM 

74 

EXOPSUM 

75 

EXOPSUM 

76 

EXOPSUM 

77 

EXDRSUM 

70 

EXDRSUM 

79 

S XORSUM 

SO 

EXDRSUM 

91 

EXOPSUM 

02 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

r, 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SINCe  ^09  eXTPUSION  op  DPAPING*  THR  OUTPUT  STPRSP  OlSTPiSUTfON  PPOP 
MATPtX  ITCPATfOM  PCTNOD  APR  fN  Of'PfOFfwCC  OP  A WJPOPM  HVOPOSTATtC 
PPCSSUPR  PPOM  THE  ACTUPE  SOLUTION  OVPP  ALL  ELCMfNTS*  THIS  PPOOPAM 
CALCULATES  THE  AVEPAGE  EXTPUSION  OP  OPAXtNG  PPE5SUPE  AND  THE  COPPEC* 
TION  PACTOP  OP  THE  NEAN  STRESS*  THE  OUTPUT  STPPSS  AND  DIP  PPESSUPE 
OISTPfPUTfONS  OP  this  program  ARP  THE  ACTUAL  SOLUTIONS* 

•PPHOTE  THAT  ALL  OUONTITIPS  INVOLVING  VOLUMN  OP  SUPPACE  INTEGRAL  APE 
measured  pep  PAOIANCt*e*  OUANTtTTES  DIVIDED  PV  |?PPA1I|*PPA 
APPNOTE  ALSO  ALL  OUTPUT  OP  STRESSES  AND  DIE  PRESSURE  APE  IN  TERMS 
OP  INITIAL  VIELO  STRESS  PPAPAPAPPPPPPPPP 
****«•««***«*«««*«*****•*•***«******•«*••«>•••••«******«*•••#««**••*•** 

Ot MCnSION  RIEAEI t?fPAPI*UR(20t«UZI20l,P(20>*EPSTRN( 2 tO  >• lELf  210*41 

1 »P7r  20l«PRI20l*STSf 4*21 01 *eP$Tpsf 2101 »ICALC20) »STSCRf4*2lO»« 

2 U01C(20ltNPf20l*COOEC2101 

*4**444P**44*44444«4*4«4*  44*4*«*********pp«ip4pp«**«««p4*p«ttP4ttM*ppp*** 

READ  PROGRAM  AND  GEOMITRV  CONTROL  CAPOS 
tEXDRvI •••AXI^SVMMETRIC  EXTRUSION  PPOBLEM 
0***AXT*SVPMETPf C DRAMING  RROPLEM 
NUMNPsTOTAL  nodal  POINTS 

numel^total  elements 

NtPTS«TGTAL  LINES  IN  MESH  SYSTEM  PARALLEL  TO  AXIS 

NOPvTOTAL  input  CAPOS  OP  SUPPACE  NODAL  POINT  PDPCC  PROM  OUTPUT  OP  MAIN 
PROGRAM* 

THCTA^SEMt-INCLUDEO  CONIC  DIE  ANGLE 
V|ELO«tN|TIAL  VIELO  STRESS 

***P4P*4«44P44««444444**«4**P44««*4ft*PP444«p*«p**p««PP*ppA***pp«4*P*»p 

READ  2300, lEXDR 

READ  2311 ,NUMNP*NUMEL*NtPTS*NPP •THETA 
READ  3312,  VIELO 
PAI ■A*4ATANI1 *1 
ANGaPA  1 4THCTA  / 1 80« 

RCOSaCOSlANGl 

RStNaSlNIANGI 


READ  NODAL  POINT  DATA  AND  ELEMENT  DATA  CAPOS 


DO  too  t•l,NUMNP 
N«| 

100  READ  1000,N,CODEfNl «P<NltZfNI 
REX  ITaRINUMNPI 
RENTEPaRCNlPTSl 

VCXlTa-f REN TER/REXITlPf RENTER /REXXTI 
N«0 

130  read  1001,M,|  IELI  M.  1 1 ,X«|  ,41 
140  NaN4l 

lNfM.Nl  170,170,150 
150  DO  160  J»l,4 
160  lELIN, J1«1EL<N-1 , J141 
|70  tP(M-NI  100,150,140 
mo  TPCNUMEL^NI  190,140,130 
190  CONTINUE 


READ  input  DATA  CAPOS 

ALL  input  DAT as  BELOw  APE  PROM  THE  OUTPUTS  OP  THE  MAIN  OR06RAM 

EPSTRNaMPPECTlVE  STRAIN  PATE 

STSf l,Nl«R-STRCSS  OP  N-TM  ELEMENT 

STSf 2,N|aZ-STReSS  CP  N-TH  ELEMENT 

STSI 3,Nl«TH-STReSS  OP  N-^M  ELEMENT 

STSC4«N|aRZ-STRCSS  OP  N-TH  ELEMENT 

EPSTRS-EPPECTlve  STRESS 

AX,  AY  ,AZ,AM  STAND  POP  INPUT  OATAS  MHICH  ARE  NOT  USED  IN  DOING  THIS 
PROGRAM, 

UR,UZ,PR,PZ  ARE  THE  VELOCITY  AND  POPCES  AT  NODAL  POINT  IN  P AND  Z-OIP* 


READ  1 003,1 IAX,AY,AZ,AM,SPSTRNrNl I ,N«1 • NUMEL I 

READ  100  3, 4 f (STS I I«NI,  I« I , 4 1, EPSTRS (N 1 , AX I •N* | • NUMEL 1 

00  220  laltNEP 

220  READ  1D02,NP( I 1,UPI ! 1,UZ( I l,PR( n .PZI t I 
1P( IC«0R,E0«1 t PRINT  3303 
tP( TBXDR*eOeOl  PRINT  330A 
PRINT  2201 
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1 

tXfW  SlIM 

P3 

PP1  NT  2202«(IC«COOEC  tf  1 •RlKltXCKl.Ksi  «NUMNP1 

ft 

5KOPSUM 

S4 

POINT  9204 

1 

•KOP^UM 

PS 

00  250  TPt •NRP 

1 

rx^'psup 

PE 

rei  llvMpf  n/MlFLO 

• PICWSIIM 

ST 

Tri  ll•■^m/vI■l.o 

f 

t Vf>K  SUM 

ss 

JsNPf I 1 

ex OP SUM 

EP 

2S0  PPINT  K203tNPII  I.CnDFIJI.UBf  n ,1<TM  I.PP(T  I .P7<  1 1 

e xoPsuM 

90 

C 

MXDPSUM 

PI 

c 

FXOPSUM 

P? 

c 

rpTOTLMTOTAL  DIE  ECRCE  TN  r-OIPfCTION  SESOPr  COPPEC^ION 

SxoppyM 

9X 

c 

CXOO SUM 

P4 

c 

E xnpsuM 

PS 

EPTCTL«0, 

EXORSUM 

96 

00  300  f«|«Nes 

rxopsuM 

9T 

J«NS( I 1 

EXDR SUM 

9S 

fMlCOOEC Jl«EO«5«0*CP«CCnE(Jl.EO*lS*OI  60  TC  301 

EXOP SUM 

OP 

60  TO  300 

EXOPSUM 

t 00 

3©l  PRTCTLPPBTOTC^PTI  n 

f 

EXDP SUM 

101 

300  CONTINUE 

r 

CXDPSUM 

loa 

c 

E XOP  SUM 

103 

c 

ft  ) exop.^uM 

104 

c 

METMOO«l...tM  ENEP6V  8ALANCE  Mchoo  TS  USED 

I FxnnsuM 

1 08 

c 

2»*«tS  2EP0  total  7*E0PCE  AT  PtOlC  ROUNDAPV  msthCO  tS  USED 

EXDP SUM 

106 

c 

EXOPSUM 

tor 

c 

. 

E XOP  SUM 

1 os 

OC  *00  M«TM0O«l  ,2 

EXDP SUM 

lOP 

tPf  METHOO«EO«n  60  Tr>  901 

E XOP  SUM 

no 

IMf  MeTHCO.EO.2  ) 6C  TO  P02 

EXOPSUM 

in 

c 

EXOPSUM 

11? 

c 

E XDP  SUM 

113 

r 

44PENEP6V  8ALANCE  METHCC*« 4**M*9«p*« 

EXOPSUM 

t 14 

c 

CALCULATE  dissipation  ENEPGV  PATE  DIE  SOOPACE  PUE  TO  EPTCTION 

EXOPSUM 

MS 

c 

EXOPSUM 

lit 

r 

EXOPSUM 

llT 

901  continue 

EXOPSUM 

1 1 9 

fpictnmo# 

EXOPSUM 

1 IP 

00  103  I»1*NPP 

EXOPSUM 

120 

J«NM|  I 1 

EXOPSUM 

121 

tricoDSi ji*eo«s«o«rp«coDEf ji«Fo*}E*oi  go  to  102 

' 

EXOPSUM 

122 

GO  TO  103 

E XDP  SUM 

123 

1 0?  MPlCTN»Mp|CTN>|MPCn4PSIN*r2f  X fPRCOS  1*(UPC  1 imps  IN+Urc  I I4PC0SI 

EXOPSUM 

12* 

103  CONTINUE 

EXOPSUM 

125 

c 

EXOP  S<iM 

126 

c 

' 

sxoosu** 

!?▼ 

c 

CALCULATE  TOTAL  STRAIN  ENEPGV  PATE  OVEO  tEM  OXFOPMEO  oCFTfON 

EXOPSUM 

126 

c 

EXOPSUM 

129 

c 

EXOPSUM 

130 

osengv«o • 

EXOPSUM 

131 

00  210  N»1 tNUMEL 

EXOPSUM 

132 

EFSSmEESTPNIN 1 

EXOPSUM 

133 

IP(EF«TRN|NI«LE.O.on  EFSS-0.0 

EXOPSUM 

134 

M«ICLCN«n 

EXOPSUM 

135 

I2»  IMLIN.2 1 

EXOPSUM 

136 

I3«J5Lf N»3) 

; 

EXOPSUM 

13*^ 

T4«  TELCN«41 

: 

EXOPSUM 

1 3S 

VOLI«  < Tf  12  1-71  n nxf  p(  !1  1*PC  T IMPf  r21»P<  I2I>X(  mspf  1?|  1 

EXOPSUM 

1 39 

VOL  ?■(  71  I4»-2II3n»fP|  14l*or  lAMPf  I31*RC  131  ♦Of  XAl«e(  13  > 1 

r 

EXOPSUM 

140 

V0L3«I  7f  13  1-71 12n*fPCT31»PfI?1«'Pft21*P(y2nP(l31MPCI211 

EXOO SUM 

1*1 

V0L4«f  7f  t 1 1-7C14I  1*(e<  t41*Pf  141  ♦oftl  >*PC  !1  MPf  11  |PRf  141 1 

» 

EXOPSUM 

142 

VOLUMN«f  VOLnVOL?*VOL36VOL4|/6. 

EXOPSUM 

1 43 

DSSNGVsOSENGV«EMSTP5fN 1«EESSPV0LUMN 

\r 

EXOPSUM 

144 

210  CONTINUE 

EXOPSUM 

148 

c 

‘ 

FXOOStfM 

146 

c 

EXOPSUM 

1*7 

c 

calculate  AVERAGE  EXTRUSION  PSRSSUPE  CP  AVEPAGE  OCAMlNG  PRESSURE  fEXT- 

‘ 

EXOPSUM 

14P 

c 

PSS1  AND  CORRECTION  MVOPOSTATXC  ©OESSUPE C ALAKOA 1 RV  ENEP6V  ELANCE- 

EXOPSUM 

14P 

c 

METWOO. 

i 

EXOPSUM 

I SO 

c 

1 

EXOPSUM 

181 

c 

EXOPSUM 

152 

PRINT  3300 

■ 

EXDP  SUM 

1*^ 

tPC  tCXOP.EO.I  1 CO  TO  bo2 

1 

EXOPSUM 

IS* 

501  €XTeSR*<0SENGV-FPTCTN|42./<-VEX|T4QExITPPCXTT> 

I 

EXOPSUM 

ISS 

ALANOA«f2*»FPTOTL-EX‘^PS5PPEKlT»PxX!T|/CPENTxPPRENTER-PEXtT*PEXlT1 

• 1 

SxnQfUM 

ISS 

GO  TO  503 

,1 

EXOPSUM 

1ST 

80?  EXTPSS*  t OSENCV-FP  TCTN IP 9, / 1 PXNTMOXPFNTEP 1 

EXOP  SUM 

189 

ALAND A«( 2.MFqTrTL-CXTPS£l/l PENTFOmoEnTE P-oEx IT TS ex  IT  1 

EXOPSUM 

159 

GO  TO  «03 

t 

E XOP SUM 

160 

c 

EXOPSUM 

161 

002  CONTI*. UE 

EXOPSUM 

1 S2 

c 

EXOO  SUM 

1 S3 

r 

EXOPSUM 

|44 

c 

♦ xp*44M»T-fOO  OF  7?P0  TOTAL  E-Ft^^E  fVEc  THE  P|GIR  •OUNCAPVPPPPPPPpPP 

1 

EXOPSUM 

165 

c 

T 
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CKOffSUM 
f KO»$UM 

e»f>»suv 

e«o*suM 
exD*suM 
evtMsuv 
e vo»si>M 
exonsuM 

exD»suM 
rve* SUM 
eXDRSUM 
exn*suM 
evDMsuM 
exovsuM 
exoMtuM 

fXPMSUM 
CKO*  SUM 
rxoPsuM 
exoPsuM 

KXOP  SUM 

rxDPsuM 
exDPsuM 
EXOPSUM 
EXOPSUM 
EKOP SUM 
EXOPSUM 
S XOP  SUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
SXOP SUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOP  SUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOP SUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 

EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPStIM 
CICOPPUM 
EXOPSUM 
EXOPSUM 
EXOPSIP* 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EKfMSUM 
EXOPSUM 
EXOPSlP* 
EXOPSUM 
E XDPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
EXOPSUM 
t XOP  SUM 


IPS 

1ST 

tss 

ISP 

ITO 

|Tl 

ITS 

ITS 

ITS 

ITS 
ITT 
ITS 
ITS 
IPO 
IPI 
IPS 
I PS 
IPS 
IPS 
I P6 
1ST 
ISP 
IBP 
IPO 
IPI 
IPS 
IPS 
IPP 
IPS 
1P6 
IPT 
I OS 
IPP 
POO 
EOl 
EOE 
EOS 
E04 
EOS 
E06 

EOT 

Eoe 

EOP 

210 

Ell 

212 

213 

214 

215 

216 
2IT 
2IP 
21P 
220 
221 
222 
EPS 
224 
EES 
224 
22T 
22P 
2tP 

250 

231 

232 

233 

234 

tss 

f3P 

23T 

tss 

23P 

240 

241 

242 

243 

244 
24f 

245 
24T 
24S 
24P 

fse 

251 
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ESS 
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iPCt2XDP.E0.il  00  re  tos 


CMSex  ECEMCMTS  NEXT  Tn  PfSIO  POUMOAPT  4T  ENTPANCf  POP  0P4M1M6 

00  TOP  l•ltNt 
JM| 

TOO  tPCPMSTPNOl.LE.O.OII  60  TO  TOT 

icALf  n«j 

60  TO  TOP 

TOT 

60  TO  TOS 
TOS  eON  TfutUE 
60  TO  POO 

CMECX  CLE^^ATS  next  to  Ptcio  POUsrAPv  AT  EXIT  POP  SXTPUStON 

• ••«••••*  ft*  ftPvXiVOXlOTPPPP***  A* 

TOS  DO  TIS  1«1 .N1 
J«NUMEL«Nl4| 

714  IPCEPSTPNC J|.LE*0.0II  60  TO  TIT 
tCAUttlMj 
60  TO  TIP 
TIT  Jaj-Nl 

60  TO  T14 
TIS  CONTINUE 
POO  CONTINUE 

ftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftaftftftftftftftftftftftVftftftftftftsftftpftftftftftftftftftftOftftft 

CALCULATE  total  E^PQPCP f PEPOPE  COPPECTIONI  POP  ELEMENTS  NEXT  TO 
PIGtO  bounoapv 

***aftftftftaftftftftftftftftftftftftft*ftftftp«ft**ftftft**ftft*ftpftaftft*ftft**ftftftftftftftftftft*ftftftftftftftftft 

TPCPCEaO. 

iMl 

INOCxal 

SOI  XXMlCALCtl 

tl*TELfKX,ll 

T2«rELfKK,2l 

ISalELCKX.S) 

|4«|ELCKX«4t 

ZlMQ. 

22af2f Ill*2C12lt/2. 

PEafPf  tSIftPf 141 1/2. 
tPCtNOSX.EO.il  60  TO  SOS 
APEftl«CP2PP2-'Pt4Pl  1/2. 

APB42«P|PZ2 
60  TO  SOP 

SOS  PtMCPCf t >4Prt2>>/2. 

APEA|«CP2PP2-»P1PP1  I/p. 

APE42«0. 

SOP  TM0PCE«TP0PeEAAPEAlPSTSI2.XXI4APEAP*STSC4«XK| 

TPCt.EO.NlI  60  TO  BIO 

P|»P2 

XABOVSaXXftI 

tPCK4S0VE-lCALC  !♦!  1 1 TES.TSOtTS? 

TSO  iPtftl 

tNOEXPt 
60  TO  SOI 

TSS  fPC SAP0VB.LT.1C4UC ITI 1 1 XNEXTaS AMOVE ftNl 
tPCXASOvB.ST.ICALf |T|I 1 XNE XT«X ABOVE*Nt 
tlalELCXNEXT*! t 
12*  fELCXNEXT«2l 
X3«  reLCXNEKT* Jl 
t4alELCKNexT.4l 
PP«CPC II4PC2IAPC3IAPC All/A. 

APEAlftlPfPPS-Pt  PPI 1/2. 

APEASftPlPEt 

APCA2«P2P22 

▼P  OPCCvTPO  PCBtSTSCP  •KASOVEIPC  APE  A1«>APE  AS  1ASTSC4tK  ABOVE  ISA  PEAS 
P!«Pf 

EEPCECIt l•ECI2ll/2• 
riPPE 

tPCXNEXT  .fO.  tCALCfAIII  60  TO  SOS 
XASCVEpXNEXT 
60  TO  TSS 
SOS  tafftl 

INOEXat 
60  TO  POt 


calculate  COPPEerlON  MVDPOSTATIC  PPPSSUXE  pv  PEAUAING  total  x*POPCE 
rOUAL  TO  TEPO  APTfP  COCPPCTTON  CALANOAI.  CALCULATE  AVEPACf  EXTPUStON 
MPfESUPE  OP  AVEPACE  0PAVTN6  PPEtSUPf  CfXTPSSI  PT  POPCf  BALANCING. 

ftaftftpaftftft*ftftftftftftftftftft«ft««**«ft*ft«ftftftft*ft*ft*ftftft**fta«*ftft»ftftftftftftftftftftftftftftftftftft 
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exDRSUM 

2SS 

c 

exoosuM 

2S6 

610 

PRINT  3301 

eicovsuM 

25T 

PRINT  2307,(|CALI II .1 ■! .Nil 

exDnsu»« 

2Sa 

TPIiexOR.eO.ll  GO  TO  611 

eirottauM 

282 

4U4NDA»-CTE0RCE/(REKrER9RENrERl 1*2« 

exoosii« 

E60 

exTPSl»-ALAN0A*«ReNTeR*»eNTeR-RexlT*RexIT)A2.*PRT0TL 

exonsuM 

261 

exTPss»exTPsi/«RExiT*Rex  iti 

exoRsuw 

2D2 

GO  TO  803 

cxonsu*! 

263 

61  I 

4L4N04v*f  TRORCE/IREX  JT9REXI7I  >92« 

EXD«  SU«) 

264 

exTPSS»-ALANOA*«PENTeR*ReNTeR-REX  ir*PCXlT)*2.*PRT0TL 

exDRSUM 

268 

503 

CONTINUE 

EXDASUM 

266 

PRINT  230t.exTPsS 

EXM  SUN 

2ST 

PRINT  2302«4L4NO4 

EXOASUM 

266 

c 

EXDR  SUN 

262 

C *9*1 

EXDRSUM 

2T0 

C CALCULATE  THE  ACTU2E  DIE  PRESSURE 

exoosu** 

2Tl 

C 9*91 

1 

1 

9 

1 

9 

9 

« 

EXDRSUM 

2T2 

C 

EXDRSUM 

2T3 

PRINT  3316 

EXDRSUM 

274 

00  400  1«1«N09 

EXDRSUM 

278 

JsNP< I 1 

EXDRSUM 

276 

IPICOOei Jl.eo.8.0.CR.COOE<JI.EO.I6.ei  go  to  AOO 

EXDRSUM 

277 

GO  TO  400 

EXDRSUM 

278 

409 

J2»J*NIPTS 

EXDRSUM 

272 

ji»j-nipts 

EXDRSUM 

260 

TPC  JI.LT.NIPTS)  JIaJ 

EXDRSUM 

261 

IPCJ2.GT.NUNNPI  J2>J 

EXDRSUM 

262 

R1  BfR(  JI^RC  J211/2. 

EXDRSUM 

263 

R2«(RfJlt*R<JII/2. 

EXDRSUM 

264 

H><2(JI  l-ZCJ2t  1/2. 

EXDRSUM 

268 

AREAsIR 1 7221*80021 ( R2-R1 l»l R2-R1 1 AHPHI /2. 

EXDRSUM 

266 

PB  ( CER(  1 19RCOS<>E2f  T 19eslNt/AREA  UALANOA 

E XD»  SUM 

2RT 

PRINT  2303.  J.COOelJl.RlJI.riJI.  AREA.P 

EXDRSUM 

266 

400 

continue 

EXDRSUM 

269 

c 

EXDRSUM 

290 

C 999< 

EXDRSUM 

291 

C CALCULATE  THE  ACTUPAL  STRESS  OISTRIPUTION 

EXDRSUM 

292 

c 9991 

EXDRSUM 

293 

c 

EXDRSUM 

294 

IPIMETHOO.CO.l  1 PRINT  3300 

EXDRSUM 

295 

tPf METH0D.E0.2)  POINT  3301 

EXDRSUM 

296 

PRINT  2304 

EXDRSUM 

227 

DO  *l>5  Ns]  .NUMEL 

EXDRSUM 

296 

ST  SCR  f 4 • N 1 « ST  S f 4 f N 1 

EXDRSUM 

299 

DO  406  l»t«3 

E XDR  SUM 

300 

406 

ST SCO ( t ,Nl>STSt 1 fNl^AUANDA 

EXDRSUM 

301 

AVES*<STSCR< 1 fNI^STSCR*  2tN>*STSCPf 3*N> 1/3- 

EXDRSUM 

302 

405 

PRINT  2305 «N, (STSCR( I tNl, 1« 1 • 4 1 , XVESt EPSTPSIN » » EESTRNIN  ) 

EXDRSUM 

303 

900 

CONTINUE 

EXDRSUM 

304 

C 

EXOR SUM 

305 

1000 

PORMATf  I5tP5«  0«SPt0.0| 

EXDRSUM 

306 

loot 

rORMATI 16151 

EXDRSUM 

307 

1002 

PORMATf 19<4F17«0I 

EXDRSUM 

306 

1003 

PORHATl 6610.01 

EXDRSUM 

309 

2201 

PORMATI/.8x.*NOOAL*tSX.*COOE*.SX.*R-COORO...*'5x.*Z-COORO...*l 

EXDRSUM 

310 

2202 

90QMAT(5X«lS.SX,P5.2«5X,Fl0.6f5X  ,F10.6I 

EXDRSUM 

311 

2203 

P0RMATI2X.  18.91  2X.P10.e)l 

E XDR  SUM 

312 

22  04 

PORMAT1///.8X,*  input  OAT  a*/.  2X.  •nodal*  . 6X.  *0008*  , 6 X.  *6- VELOCITY* 

EXDRSUM 

313 

1 

1 .•2X  . *2-  VELOC 1 TV*  . 2X  .*  p-  POPCE  ...*.2X.*  Z-POOCE*  1 

EXDRSUM 

3I« 

2300 

FORMAT C IS) 

EXDRSUM 

3tS 

2301 

PORHATl //,SX,*AVERAGE  EXTPUSIONI  OR  ORAVINGI  PRESSURES*. P 12.61 

exDPSUM 

316 

2302 

FORMATf //.5x,9TME  CORRECTION  HYDROSTATIC  PRESSURE** .F| 2.6 ) 

EXDRSUM 

317 

2303 

FORMAT I2X,*N0CAL  POINT b*, I A« 2x , *C00E* 9F5, 2«  2K , 4R.COOPO«*«Ft 0 .6* 

e»oi*su« 

316 

1 

1 2X,*Z-COORO.s*P|0.6.2X.*APSAs*.P|0.6.2X.*DIE  PRESSURES*. PIS. 81 

EXDRSUM 

312 

2304 

PORNATI/.IOX.*ACTURAL  STRESS  DISTRIBUTION*// 

evp^suM 

320 

\ 

1 •5x,*CI.E.N0.*.5X.*Rpp5rR6SS*,5X,*2-STPESS9,SX.9rH.STPFS6** 

exoosuM 

321 

f SX.*RZ-STRESS*,*..mEAN  stress. .EP-STRES6*.3X.*EPP-STRATN  RATe*l 

EXDRSUM 
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230S 

F0RMAT(7x«I6»7F13.6 1 

EXDR  SUM 
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EXDRSUM 

324 

231  1 

pormatiais.pio.oi 

E XDR  SUM 
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3300 

FORMAT!  1H1  .SX.FENERGV  BALANCE  METHOD*) 

exoosuM 
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3301 

FORMAT! IHltSX.FTOTAL  7ERO  Z-FOPCE  AT  RIGID  FOUNOARV  METHCD*) 

EXDRSUM 
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3303 

FORMAT! 1 HI .*••• .EXTRUSION  PROMLCM , . . «* ) 

EXDRSUM 
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FORMAT! IHl «9...«DRA«TNG  PR OSLEM. • . ) 

EXDRSUM 
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331  0 

FOPMAT! ///fSX.9THE  ACTUAL  DIE  PRESSURE  DI STR IBUT ION* ) 

EXDRSUM 

330 

3312 

FORMAT!F10.0) 

EXDRSUM 

331 

STOP 

C XORSUM 
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END 
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EXTRUDE 

EXTPtJDE 

EXTRUDE 

extrude 

EXTRUDE 


2 

3 

4 

5 

6 
7 
A 
Q 

10 
1 1 
tE 

13 

14 

15 

16 
IT 
IB 
IR 
20 
21 
22 

23 

24 

25 
76 
27 
2B 
20 

30 

31 

32 

33 

34 

35 

36 

37 
3H 
30 

40 

41 

42 
*3 

44 

45 

46 

47 
4A 

40 

BO 

Et 

52 

53 

54 

55 

56 

57 
3A 

50 
60 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 
7! 

72 

73 

74 

75 

76 

77 

78 

79 
RO 

51 
62 


PROGRAM  PURTI TAP61 • t NPUT tOUTPUT •TAPE5 •! NPUT •TAPe6«0UTPUT « PUNCHl 

C 

c ♦•♦♦♦•••••••♦••••♦•••4P44404RP644P44P44444R44444PP4P44PPP44449PP4P8PP4 

C PROGRAM  POP  STEADY  STATE  E iTRUSION  OR  DRAPING  BY  MATRIX  METHOD* 

C NOTE  THAT  THE  POLLONUNG  PROGRAM  IS  WRITING  ACCORDING  TO  THE  RIGHT 
C ORDERING  OP  NODAL  POINTS  AND  plEMENT  NUMBERS  fSEE  EXPLANATION! 


COMMON/GENCON/NUMNP,NUM£L*H£Df  |2) •VOL*NEO*NS*  t TERNO* ISTOP* 

1 YlELDtMBAN0*TEST«M01 AGt NHP, NUMPC*  f TER« NRP 2* NSC ALE tNPUNCH*NPR I NT * 

2 NCHFCK *ACFINI 

common  /wall/  theta, PTtTANTH 

COMMON  /STRPATH/  ST EP , Y STAR T , r D I F , Yg X I T , V M 1 N, P ENTER , REX  I T , VEX IT 
COMMON  /DIMEN/  M | «M 2, M3 .MA ,M5 ,M6, M7, MB « MO ,M10 * Mf f tMI 2* M| * ,M| 4 

C 

C PROGRAM  PUOT  IS  PCR  CONTROLLING  THf  DIMFNSION  OP  THE  COMPLETE  PP06RAM 
C ITS  PURPOSE  IS  TO  PREVENT  ASSIGNING  A LARGER  THAN  NECESSARY  DIMENSION 
C POP  ANY  ARRAY  THROUGH  THE  USE  DP  POLLOWING  STATEMENT*** 

C 4**4«44iW4i4t««4*W«»*»4*«»«44«4«4«*«*»PP4*4444444*P4«'W«W44»4**4«4P4P#44P9 

C 

common  Af246S01 

c 

c *W4i*«4«WP4«P4W«««PPPPP*P*P4*4*»*P*P4PP4»**4P44P*4«4«W44*44**44**4P4P«W 

C NFIELO  IS  THE  dimension  OP  ARP AV  A*  ITS  VALUE  CAN  BE  DETERMINED 
C PRECISELY  BY  RUf.NING  THF  PROGRAM  DNCt  • 

C *«4*««444**4*«««»4«A*4*i»W44*4»«W4t««44444W4W<i*««**#«4*4*«*«4*P4444W4W4P 

C 

NP IFLO«2A650 
OAI s4*«ATANC 1 *0 1 

r 

C **♦»**« 

C THfcTA 
C YSTART 
C 

c 

C STEP 
C YD  IE 
C VEX  IT 
C VMI  N 
C RENTER 
C OSXI7 
C VEXIT 
C NTIMFS 
C NIPTS 
r NJPTS 
C NMAX 
C **«44**i 

C 

THE  TA  «6  • 

TMETA*TMETAWPA| /I60. 
tanth«tan<theta I 
YEX|T*0. 

OFNTCR«1*0 
NTIMES*! 1 
N|PTS»7 
NJPTS*15 
NMAXbTSO 

c 

C READ  THE  INPUT  DATA  CCNTROL  CARDS 

r HEDbOUTOVIT  title  up  to  72  characters* 

C ACPINiBiNiTf AL  OECCELARATION  COEPPICINTS,  POR  PtoST  ITERATION  ONLY* 

r iteobtotal  number  op  iterations  assigned* 

C ITCDNTbI,  fop  NONWORKHAROEN  ing  materials 
C Of  POR  WORK-HARDENING  MATERIALS* 

C NPUNCHbO,  IP  NOTHING  OTHER  THAN  VELOCITIES  TO  BE  PUNCHED* 

C 1*  OTHERWISE. 

C NPRINTbI,  IP  THE  NODAL  POINT  OATAS  TO  BE  PRINTED* 

C 0*  OTHERWISE* 

C PLIMXTbvalUE  of  CEPROR  NORMl/f SOLUTION  NORM|  REQUIRED  FOR  PINAL  RESULT 
C NUMNPsNUMBEP  op  nodal  POINTS! total  1 • 

C NU»CLRNUMBSP  OF  ELE  ME  FT  S I TOTAL  1 

C NUMPC«NUNBER  OP  TRACTION  BOUNDARY  CONDITIONS  CARDS  TO  BE  READ* 

C NBPbNUMBER  op  nodal  points  at  which  PORCE  calculations  are  REQUIRED* 

c *W«4«««*4*«««**«***««4P«*«««***«***4*P*****««4«««*4«**4*«4*«4*4*4*4R*4 

C 

READ!5*t000l  HED 
READ! 5* inOAIACPINX 

REAO(5«1003 ! ITER* ITCONT , NPUNCM, NPR INT , PL IM IT 
READ! 5* lOOSl  NUMNP* NUMEL* NUMPC ,NBP 

93 


SEMI- INCLUDED  CONIC  DIE  ANGLE* 

/-COOPOIKATc  OF  NOTING  POINT  TC  CALCULATE  TOTAL 
EFFECTIVE  STRAIN  AN^  IS  SETTING  SfGMTLY  LESS  THAN 
THE  CENTFO  OF  FIRST  ELEMENT* 

THF  INCREMENT  SIFF  FOR  STRAIN-RATE  INTEGRATION 
7-COOPDfNATE  AT  ENTRANCE  TO  THE  DIE 
/-COORDINATE  AT  EXIT  FROM  DIE 

MfMMUV  /-COOROINATE  OP  THE  CONTROL  VOLUME 

RADIUS  OF  BILLFT  PEPQRE  FXTPUSION 

RADIUS  OP  BILLET  AFTFR  EXTRUSION 

EXIT  VELOCITY  OP  THF  PRODUCT 

NUMBER  OP  FLOW  LINES  TO  BE  CONSTRUCTED 

NUMBER  DP  LINES  IN  MESH  SYSTEM  PARALLEL  TO  AXIS 

NUMBER  OP  LINES  IN  MESH  SYSTEM  PERPENDICULAR  TO  AXIS 

TOTAL  NUMBER  OP  INCREMENTS  POP  INTEGRATION 


EXTQUOf 

03 

NtOT5»N|OT5Tl 

exTQUO* 

EA 

N692»M9292 

rxToooe 

es 

N9«««90 

sxrouDS 

04 

N82«N922 

eXTROOE 

8* 

NELsNUMEL 

ExTounf 

09 

NPCoNU^oC 

6KT»UnE 

•9 

N«Qs39NU9MO 

EVTBUDS 

00 

NOoNEO 

EXTPUOe 

91 

Nt«KlPT5 

€XT»UOE 

9? 

NJSNJOTS 

•XTPUOf 

93 

C 

e»TP  UDE 

94 

c 

CKTPU05 

9S 

c 

OETeO'iJMF  THE  location  CF  THE  STAPTIMG  POINTS  OF  DIFFERENT  ARPAVS  ON 

SXTPUOS 

94 

c 

AOOAV  A 

eXTPilOE 

9? 

c 

PXTPUOc 

90 

c 

ex  TOUCH 

99 

N1  •! 

exTouoe 

100 

N2«N| ♦NUMNP 

exTouof 

101 

N3«N24>NUMNP 

exTOUCH 

102 

N9«N34NU9NP 

sxTour^e 

103 

N5«N9TNUNNP 

exTouce 

1 09 

N6aK54NUMNP 

SXTO(iOE 

lOS 

N7«N6499NUMeL 

exToune 

106 

NOvNT^NUMEL 

HXTOUOC 

lOT 

N9«N946*NUNEL 

CXTOUOH 

106 

N|0«N9^5«NUMEL 

E XTOUOe 

109 

NI1«N104>SONUNEL 

exToime 

110 

N12«N1|»NBF 

exTouoe 

111 

Nt  3«Nt24N8F2 

exTovioe 

112 

N 1 9 «N 1 34  29  NUMPC 

EXTOUOE 

113 

NI5*N|9499NUMPC 

EXTOUOE 

119 

Nt  6«N1 54N0F 

exTOUOE 

116 

Nt7«N1E4N8F 

ExTOUOE 

1 le 

N19*N174N9F 

eXTOUOE 

117 

N19«Nld4NTtNcS 

EXTOUOE 

1 16 

c 

EXTOUDH 

119 

c 

EXTOUOE 

120 

c 

OSAO  THE  input  data 

EXTOUOE 

121 

c 

EXTOUOE 

122 

c 

EXTOUOE 

123 

CALL  POELIM  C AINl },A(N2»,ACN3I .AINAl ,A(N51 •A(N6I .AtNIl I .AINISI, 

e X TO  UO£ 

129 

1 A(NIA), A(Nt8l«NELf NPC»NT |ME St N IPTS • TSTAQT • VK IN , R£ XI T) 

EXTOUOE 

128 

VOt  «• ( PENT6R-OE  XI T 1 /TAKTH 

EXTOUOE 

126 

VEXITo-IOENTEO/OExlT  >»9  2 

EXTOUOE 

127 

ST» P«  f V ST AOT-VM I N 1 / ( NMAX 9 1 • 1 5 1 

EXTOUOE 

129 

0O|TE<6, 10071  ACElNl .ITCONT.FLIMIT 

SXTOUCE 

129 

W9ITEC6* 1006 1 srro, TSTAortTOfE#  TEX IT, TNiN.PENrsR.oex/r, vsxi t* 

EXTOUOE 

1 30 

1 NTtNES«NlPTS*NJPTS*NNAX 

EXTOUOE 

131 

c 

E XTOUOE 

132 

N20«N1o^nE0 

EXTOUOE 

133 

c 

EXTOUOE 

139 

c 

ExT»UOC 

138 

c 

EXTOUOE 

136 

c 

since  AFTEP  The  SOLUTION  OF  stiffness  equations*  the  stiffness  matrix 

EXTOUOE 

137 

c 

IS  NOT  NEOEEO  FCO  that  ITERATION*  THE  SPACE  PROV lOcD  FOP  STIFFNESS 

ExTO-jne 

136 

c 

MATO  lx  FILL  BE  USED  FC®  CCNSTRUCTION  OF  FLOW  LINES  AND  FfjP  STRAIN* 

SxToure 

130 

c 

PATE  INTEGRATION.  THE  NUMB5RS  Mi  TO  X 1 4 AFE  FPO  DETERMINING  TmE 

rxTotjoe 

190 

c 

LOCATION  OF  OROPEP  AOOAV 

EXTPUO? 

191 

c 

EXTfcrUO* 

192 

c 

EXTOUOE 

193 

M|«l 

EXTOUOE 

199 

M2«M| ♦NI9NJ 

e»T*ope 

t9S 

M3aF24Nt9NJ 

EXTOUOE 

196 

M9«N34NJ 

E XTOUOE 

197 

MSaFA 4NI 

BXTOUOS 

19A 

M6aMS4N|9NJ 

EXTOUOE 

199 

M7«M64NI9NJ 

EXTOUOE 

ISO 

M8sM74NT1 me  S9KMAX 

EXTOUOE 

151 

MOaMRANTlMESANMAX 

EXTOUTH 

152 

••1  0aM94NUMEL 

EXTOUOE 

153 

M11aW|04NUMEL 

EXTOUOE 

159 

M12aMiiaN19NJ 

EXTOur^r 

155 

X|3aM12«>NJ 

EXTOUOE 

1S6 

Ml4aM134NTlMES9NMAX 

E XTOUOE 

157 

c 

•XTOUOE 

139 

c 

EXTOUOE 

159 

HOWEVEO*!^  THE  SPACE  OF  STIFFNESS  MATRIX  IS  NOT  ENOUGH  FOR  STRAIN  OATS 

EXfOUTE 

1 *0 

c 

INTFGOATIDN*  CREATES  mqrf  SPACE  AS  REOUIPEO* 

ExToyne 

161 

c 

EXTPUOE 

162 

c 

ExTO..»Ot 

163 

NNN aM«AN09NSO 

•XTPUOE 

169 

N21 aN204NNN 

EXTOUOE 

16« 

IF(  Ml  A* 5* * NNN 1 N21aK?04M1A 

EXTOUOE 

166 

N22aN2l4N9F29I 29MFANO- 11 

•XTOUOE 

|6^ 

N2TaN?24NE0 

EXTOUOE 

169 

N29aN?34N?G 

•xTRUOE 

164 

EXTRUDE 

1 TO 

EXTRUDE 

|71 

f «truoe 

177 

extrude 

I’^.T 

EXTRUDE 

124 

EXTRUDE 

ITS 

EXTRUDE 

l?6 

EXTRUDE 

177 

EXTRUDE 

|7R 

EXTRUDE 

179 

EXTRUDE 

190 

EXTRUDE 

161 

EXTRUDE 

IB2 

exTRUO? 

193 

EXTRUDE 

ISA 

EXTRUDE 

169 

EXTRUDE 

196 

EXTRUDE 

167 

EXTRUDE 

IBB 

EXTRUDE 

169 

EXTRUDE 

190 

EXTRUDE 

191 

EXTRUDE 

192 

EXTRUDE 

193 

EXTRUDE 

194 

EXTRUDE 

199 

EXTRUDE 

196 

EXTRUDE 

197 

EXTRUDE 

196 

EXTRUDE 

199 

EXTRVOE 

200 

EXTRUDE 

201 

EXTRUDE 

202 

EXTRUDE 

203 

EXTRUDE 

204 

extrude 

209 

EXTRUDE 

206 

EXTRUDE 

207 

EXTRUDE 

706 

eXTHUOE 

exmuoe 

exTEuoe 

exTAUDC 

exTAuoe 

EXTAUOe 

exT«uoc 

exTuuoc 

exTuuoe 

exTituoe 

cxTwuoe 

exTuuoe 

exTEuoe 

extuuDe 

exT»uoe 

CXT^UOC 

CXTAUOC 

EXTAUOC 

exTftuoe 

Kxrnuoe 

CXTHUDC 

CXTEUOC 

EXTWUOe 

exYwuoe 

exTEUDC 

fXTAUOe 

CXTl«UOC 

exrivuoe 

extuuDf 

ffXTmiOC 

CHT*uec 

tXTilUOC 

tXT«UOt 

CXtvuDC 

CXT»«JOC 


2t0 

ZtX 

212 

213 

214 

219 
216 
21T 
216 
216 

220 
221 
222 

223 

224 

229 
226 
227 
226 
226 

230 

231 

232 

233 

234 

239 

236 

237 
236 
236 

240 

241 

242 

243 

244 
249 
246 
24*r 
244 


296 

261 


* ^74  IS  total  SPACC  REOUIALO*  12  N2|«L0  • THE  INfTlAt  ASSI6N6D  SPACE* 
C tS  NOT  ENOUOM,  STOP  THE  PPQGAAN. 


|PCN2A*l€«NP lELO > OC  TC  100 

VPItei 6* lOOtl  N24 

STOP 

100  CONTINUE 

•AITEIOt 10021  N24 


C THE  POLLOPtNO  SUOKOUTINE  PLAST  DOING  TMf  MAIN  PART  02  MATRIX  ITERATION 
C meThOO  02  PRESENT  PROELEm. 

C •A444PRRRAR44RA444A4R4A4A4AA4AAAAAARAAAARRRRARR4ARRA4PARR4ARAMRPRM4MAA 


CALL  PLASTI ACNI |«A(NEI ,Af N3I • A|N4 I « A INS » » A<N6 > « AIN7 I • Af N6 1* A<N61« 
lAf N10l«A(N|ll«A(Nt2l« A4N13I. AIN14|« A(N|S1*A|N|6I •ACN17I*ACN16I • 
2A(NI6t  *A|N20I •A<N2t 1 * A ( N22 > , A I N23 I vNO* NELtNPC tNBE* 

32LI«fT« 1TCONT«N7IMES«N|PTS«NJOTS«NMAXI 


1000  PORMATf 12A6I 

1001  FORMATf///*  T*-E  DIMENSION  02  THE  ARRAY  (A  I IS  TOO  SMALLP/ 

16  THE  SIZE  02  THE  ARRAY  <AI  MUST  BE  6,  17) 

1002  20RMAT4//6  THE  NECCESSARY  SIZE  02  THE  ARRAY  fAI  1S6*  171 

1003  FORMAT <419*21 0*0) 

1004  20RMAT< 6210*01 
1009  20RMAT<419) 

1007  20RHATf//6  AC0C2  • 6,26*S*10K«6  ITCONT  ■ 6* 12* 

1 10X«6  2lIM1T«6«2| 0*6*//) 

1006  20RMAT<//6  STEP  SIZE  ■ 6,20*4*10X«6  YSTART  • 6*26*4*10X«4  VDIE 


1 26*4tl0x*6  YEXIT 

2 26*4, 10X*6  REKIT 

3 |4*I0X«6  N1PTS  • 
1006  FORMATC/  6 STORAGE 

I NOT  ENOUGHA/) 


m 6«26*4*//6  YMIN  • 6*26*4«10X«6  RENTER  • 6« 

« 6,26*4«I0X*6  vexiT  • 6*26*4*//6  NTlMES  • •• 
6«t3«10X«6  NJPTS  • ••I3*10X«6  NMAX  • 6*14*/) 
SPACE  AVAILABLE  FOR  THE  STRAIN  CALCULATIONS  IS 


STOP 

END 


SUBROUTINE  PRELIM  ( R* Z*UR *UZ * CODE  * lELtNF , I JBC* P$BC*RR*NEL »NPC* 

1 NTIMES*N|PTS«YSTART*YMtN«flEXlT) 

COMMON/ GENC0N/NUMNP«NUN£L*H90<) 21* VOL*NEO*NS* I TERNOtlSTOP* 

1 Y|EL0*MBAN0*TEST*M0I AG«NB2, NUMPC* |TER*NB22*NSCALE*NPUNCH*NPR|NT* 

2 NCHECK aACFINl 

DIMENSION  R<1  >* Zfl )*COOE< 1 t*UR< 1 )«UZ< 1 I* lELf  NEL  •l)*N2<||* 

1 tJBC(NPC*l)*PSBC<NPC*l)*RR<ll 


66666666666666666666666R66666666«66*6666666666666»R«66666666666 6666666 

THIS  SUBROUTINE  READS  AND  PRINTS  ALL  CONTROL  INFORMATIONS  AND 
NODAL  POINT  DAT AS  AND  BQUNCARY  CONDITIONS* 

666 66666666 666666666626666666 66 666666 66 6*6666 66 666666666666666666666*6 


6R1TE<6*1000I 

RRITE  <6*2000)  HED*NUMNP*NUMCL 


66666666*66666*666666*66*66*66666666*6*6666666*66666666666666666666*66 

CALCULATE  INtTINAL  YIELD  STRESS* 

*6*66*6666**666666666666666*66*6666666*6666*6*6666*6666666666666666666 


CALL  NARO<0**YIBLOI 
6R|TEC4*2010I  YIELD 
«RITE<S*20I 1 IITER 


6666*6666*6*66**6666666666*6666 66*666 66 666666666P666666666666 6*666666* 

READ  AND  PRINT  02  NODAL  POINT  DATA 

66*66*66*66*6**6666 66*666666666**666* 66*6**** *6*66*6*66***6666*66*6 6*6 


L60 

I2<nprint«eo»o)  go  to  60 

WRITE  <B*11IA) 

•RITE  <6«2004) 

60  READ  <S«1002)  NtCOOE< N) «R< Nl »Z<N) 

90  t2<NUMNP«Nl  100*1 I0«60 
100  WRITE  (6«2009)  N 
CALL  EXIT 
110  CONTINUE 

YSTAPTa<Z<l)6Z<N|PTSA2l)/2«*0*001 

YMlRaZINUMNP) 

6EXtT*6<NUMNP| 

fPf RPPINT.EO.OI  GO  TO  120 

WPITE  <4*2042l  <K,CODE<WI*RlR),Z<Kl,K*ltNUMNP| 


95 


^MT^uoe 
f>  BTOgOff 

««Tnjr>ff 
rBT^iinfc 
e KTBIJDC 
rxTQOOf 
exrvuoe 

exr^t^e 

F«T»»JOf 

exTttuof 

CXTOUOC 

cxrvuof 

FxT*Uf>f 

exTAUDC 

exT»uoc 

exT9i|09 

fxTocme 

«XT«UD* 

rxritcioe 

exrxuDC 

exTvune 

txTnuoe 

exTwune 

exToune 

exrvuoe 

exTvuoe 

tXTMilOC 

exTRuoe 

exTvuof 

sxTffuoe 

exTttuoe 

exrmjoe 

exT^UDC 

exmuoc 

exTAMoe 

exTvuoe 

exTmioc 

cxn»(jD* 

exTRuoe 

cxrauDc 

exTxgoE 

EXTR'jOe 

EXTOUDE 

exrnuoe 

exrnuoe 

exnvuoff 

exTQuof 

exTsgo* 

exTAuoe 

EXTdUOC 

cxTffuoe 

exTRuoe 

exTHuoe 

exT^UDC 

exTAuoe 

EXT9UOe 

exTRuoe 

exTRuoe 

exTRuoe 

CXTRMDC 

exTRuoe 

extrude 

EXTRUDE 

EXTRUOE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

extrude 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUOE 

EXTRUDE 

EXTRUDE 


WEAD  AND  Ri;1nT  SPfCIAL  AUOnC  ROUNOARX 

NRE  • TOTAL  MJRRk*  D*  NCOAL  PDI^'TS  AT  xHlCH  xroCES  ARE  DESIRED 

NE(fl  • SLOUEMTcaC  KURRER  or  NCO^L  PCiKiTS  AT  amICH  ECRCCS  ARE  DESIRED 

RRI  llsB^COOltOtNArf  OP  rcOR  LINTS  TO  Ei  CONSTRUCTED  AT  ENTRANCE* 


DO  441  • NUMRC 

RRADCS* lOOSI  IjDCf  t tl > t f JDCC I tSI t <RSRC I I « J 1 1 J« t *4 1 
441  RRirCISf 100SI IJRCI  !• I )«  1 JBC ( 1 • , CPSflCC  f • « J«t  t4l 
440  CONTINUE 

■Rireio* loosi  i nei m«i»i •nrri 
READ  tSvtOOll  IRRin*  1*1*NT  IMES  I 
«R1 TEI4*1013)  |RR( M •laf  »NTI«ES> 


read  and  PRINT  QP  ELEMENT  RPOOEPTIES 


N*0 

fxo  READ  (S»t003l  M,(tEL(P*l> 
140  NaN^l 

fP  rM«>N|  170*  170*  ISO 
ISO  DO  ISO  J«t  *4 
160  ieLlN,Jl«tEL<N-l«J|Al 
170  IP  <M»N1  160«ia0«l40 
ISO  IP  INUMEL-Nl  100*100*130 
100  CONTINUE 


1P(NPO|NT«EO*OI  GO  TO  210 
PRITE  16*20011 
DC  208  Nal.NUMEL 
208  WRITE  16*20031  N*I1EL(N«I1 


210  Jmn 

00  240  N«|,NUNEL 
DO  240  I«t*4 
DO  230  L*l *4 

XXa lABSI 1EL(N« 1 I-1EL(N,LI 1 
fPIKX-Jl  230*230*220 
220  JaKX 
230  CONTINUE 
240  CONTINUE 

RBAN0a3*JA3 

MOl  AGal 

280  WRITE  (6*11221  NEO,m0AnD«MD|A6 


C P»4*«4«44PR*4p«PP*PP4***4»*«4P4««6**4*«*«P**tt*R«*R«RW««4464«0ORR6R0R40 

C POR  EACH  element*  ASSIGN  THE  MEAN  PRESSURE  VALUE  TO  THE  NODAL  POINT 
C OP  TM*  HIGHEST  NUMDER  AMONG  THE  PCUR  CORNER  NODAL  POINTS* 

C WRR6RRWRRR444RRWP0444A4W444444R444444444PP4W4P444PP4PP44PPPRR6RRRR6PP* 


NSTOPaO 

DO  370  JaitNUMEL 

MIDaNAXOflELl J* 11 « I ELIJ* 2 1 * lELI J*3l* IEL(J*6II 
CODE! NIDI aCOOEf  M|0 I ♦ 1 0* 

IP(COOE(M1D1*LT*20«1  go  to  370 
WR1TE(6*I0I21  MID 
NSTOPwt 
370  CONTINUE 

fp(Nsrop*eo*ti  STOP 


1000  POPMATI IHll 
loot  P0PMAr(SPl0*0l 

1002  PORMAT  ( t8*P8«0*8Pt0*01 

1003  PORMATf 16151 

1008  P0RMAT(2t8«4Pt0*01 

1006  PORPATI//  4 THE  NODAL  POINTS  AT  WHICH  PQRCE  CALCULATIONS  MIt  OSSIR 
tEO«  //  20ISI 

1007  PORMATI IHI *15X*  39M  LINEARLY  OISTRISUTEO  BOUNDARY  STRESSES/ 

1 / 4 NQoe-l**NOOe-J**«PRESSURE-l***PRES$URE-J* 

2 4 SHEAR  1 SHEAR  34| 

lOOB  PORMATI  IH  •2t8«4ei8*Sl 


jM 


exTAuoe 

eXTRUDC 

cktaudc 

exrauoc 

exT9uoe 

exrauDE 

CXTHUOE 

CXTAUOC 

CXTRUOC 

extrude 

EXTRUDE 

extrude 

extrude 

extrude 

extrude 

extrude 

extrude 

EXTRUDE 

EXTRUDE 

extrude 


EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 


lOlt  RORMATC/R  NODAL  POINT  ERROR  ROR  POINT  NO*  •« t 5 t *••••••• I T CONTAINS 

I MORE  than  one  element  INPORMAT |0N«/ | 

1013  PORMATf/R  THE  STARTINC  R-COOROI NATES  FOR  STRAIN  CALCUL AT  IONS • R/ / 

\ ltP9*A»/l 

1114  PORMATIlHlt  R NODAL  POIKT  IKPORMAT ION*, // | 

1122  PORMATI///  R NNUMPEP  OP  EQUATIONS  aR  • l*t/ 

1 R EANOvIDTH  aR,  14/ 

2 R DIAGONAL  ELEMENTS  , 14  I 

2000  FORMAT! |H  12A0/ 

1 SOHO  NUMBER  OP  NCOAL  PCINtS*— 13  / 

2 30HO  NUMBER  OP  ELEMENTS————  13  /\ 

2001  format  f|Hl,40H  ELEMENT  NO*  1 J X L I 

2002  FORMAT  ( 1 1 2 .P I 2* 2 • 2^ 1 2* 3 1 

2003  FORMAT  f 1 ItSf BIO, 1 1 

2004  FORMAT  (/  ,4  NODAL  POINT  TYPE  R-OROfNATE  ?-OOOfNAT6R| 

2009  FORMAT  C20HONOOAL  POINT  CARO  ERROR  Na  ISI 

2010  FORMAT!//  R INITIAL  YIELD  STRESS  a 4,  F1S*2//I 

2011  FORMATf///R  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWEOa** 1 3 ) 

RETURN 

END 


SUBROUTINE  PLAST !R, 2* UR vUZ, CODE , ICL , YV , STS, TEPS, EPS*NF «FPUR , 
t IJBC, PSBC , PR , FZ, sot E,RR,B, A, P ST, ALANOA,GVECTR, NO. NEL,NPC, 
1NB2,FL1M1T, tTCONT,NTlMES,NIPTS,N3PTS,NMAXI 
C 

C RRRRRRRR444R4R444444444RR4RRR4R444444444444444444444444444444RRRRRRR44 
C PLAST  IS  THE  CONTROLLING  SUBROUTINE  OP  THE  MATRIX  METHOD 
C RRRRRRRRRRRRRRRRRRRRR44RRRRRRRRRRR44R4R4444444I44RRRRRRR4R4444RRRRRM4RR 

C 

C0MM0K/GENC0N/NUMNP,NUF£L,HE0< 12I,VOL,NEQ,NS, ITERNO, ISTOP, 
lVfELO,MBANO,TEST,MOf AG,NeP,NUMFC,lTER,NBP2,NSCALE,MPUNCHf NPRINT, 

2 NCHECX, ACFINI 
COMMON  /WALL/  THETA ,PT , TANTH 

COMMON  /OIMEN/  M 1 ,M2, M3,M4, MS,M6, M7,M8, MQ , Ml Q , M 1 1 .M I 2 , M| 3 fMI 4 
DIMENSION  R!1  l,Z!l  ) ,UR!I  I ,UZ!  I ),COOE!  n,  lEL!  NEL  ,ll,YV!l) 

1,STSI 6, t ) .TEPSf  5,1) ,EPS!5, I ),NP< I ) ,PPUR!1 ) , IJBCCNPC,!  ) .RRll I, 
2PSBC!  NPC  , 1 ),PR! 1 ),PZ! 1 1.SOIE! 1 I, eC 1 I, A!  NO , 1 I ,P ST! NB 2 ,1) 
DIMENSION  ALANOAtl ) ,GVeCTR! 1 I 
C 

CRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR4R4RRR444R44R4444444RRRRRRRR4R 

C INITIALIZED  ALL  YARIABLES 

C GVECTRINI  IS  USED  POP  STORING  THE  AMOUNT  OF  VELOCITY  AND  MEAN  STRESS 
C TO  BE  MODIFIED  PROM  THE  RESULT  OF  ITERAT ION, THESE  VALUES 

C ARE  REDUCED  BY  DECCELARATIDW  COEFFICIENT  IN  ORDER  TO  ENSURE 

C THE  CONVERGENCE  CP  THE  SOLUTION* 

C TEPSf 5,N|arOTAL  EFFECTIVE  STRAIN 

C A,  6,  STAND  FOR  MATRICES  IN  EQUATION  AXaB,  AFTER  SOLVINF  THIS  EQUATION 
C THE  X VECTOR  IS  THEN  STORED  IN  8 

C YVCN)«  EFFECTIVE  STRESS 

C TEPS!1,NI«  TO  STORE  VHE  PREVIOUS  VALUE  OP  YY<NI  FOR  CHECKING 
C CONVERGENCE  FOR  XORK-HARDCNING  MATERIAL 

C TEPS!2,N|aOIFFERCKCE  OP  EPPCT !VE  STRESS  BETWEEN  THE  PREVIOUS  AND  NEW 
C SOLUTIONS*!FOR  WORK-HAROEN ING  MATERIALS) 

C ITSTa  INDEX  FOR  CONVERGENCE  OP  SOLUTl ONS! I TSTaS  INDICATES  THE  ITRATION 
C DOES  NOT  CONVERGENT! 

C FSORI,PFa  TOTAL  SQUARE  NORM  OP  EQUATIONS  IN  MINIMIZATION 
C FSORS,FPFFa  SQUARE  NORM  OP  EQUATIONS  CORRESSPONOING  TO 
C INCOMPRESSIBILITY  CONO IT  IONS ! VOLUMN  CONSISTANCE) 

C OlPPapP-FSQRl,  difference  OF  TOTAL  SQUARE  NORMS  OF  PREVIOUS  AND 
C PRESENT  ITERATIONS 

CRRRRRRRRRRRRRRRRRRRRRRRR4RRRRRRRRRRRRRRRRRRRRRRRRRRR»R«RRRR4«4RRRRRRRRR 

c 

00  99S  la|.Mt4 
998  AltlaO.O 

DO  999  Nai.NEO 
ALANOAINlaO* 

999  GVECTRINlaO* 

OO  442  Nat . NUMEL 
VVfNIaf.O 
DO  442  1«1  , 5 
442  TEPSf  l,N|aO. 

PSORIaO* 

PSORSaO* 

•FFFaO. 

iTSTal 

FFaO* 

OlFFaO* 


READ  THE  INPUT  VELOCITY  FIELD 


REAOfB, 10171  fURItl,  UZCl),  laf,  NUMNP) 


I 


extrude 

A21 

EXTRUDE 

A22 

extrude 

A2S 

EXTRUDE 

A2A 

EXTRUDE 

A25 

EXTRUDE 

A26 

EXTRUDE 

A27 

EXTRUDE 

A2B 

EXTRUDE 

A29 

E XTQUDE 

AtO 

EXTRUDE 

A3t 

EXTRUDE 

A32 

EXTR UDP 

A33 

EXTRUDE 

A3A 

EXTRUDE 

ATS 

EXTRUDE 

A36 

EXTRUDE 

AST 

EXTRUDE 

436 

EXTRUDE 

A39 

E KTR«JD£ 

A AO 

extrude 

AAl 

EXTRUDE 

AA2 

EXTRUDE 

AA3 

EXTRUDE 

AAA 

EXTRUDE 

AAS 

EXTRUDE 

AA6 

EXTRUDE 

AA7 

EXTRUDE 

AAe 

EXTRUDE 

AAR 

EXTRUDE 

A50 

EXTRUDE 

ASI 

EXTRUDE 

A92 

EXTRUDE 

A53 

EXTRUDE 

ASA 

EXTRUDE 

ASS 

EXTRUDE 

A56 

EXTRUDE 

AST 

EXTRUDE 

ASB 

EXTRUDE 

AS9 

EXTRUDE 

ASO 

EXTRUDE 

A6I 

EXTRUDE 

A62 

EXTRUDE 

A63 

EXTRUDE 

A6A 

EXTRUDE 

A6S 

EXTRUDE 

A66 

EXTRUDE 

467 

EXTRUDE 

A6B 

EXTRUDE 

469 

EXTRUDE 

A70 

EX TRUDE 

ATI 

EXTRUDE 

472 

EXTRUDE 

A73 

EXTRUDE 

47A 

EXTRUDE 

ATS 

EXTRUDE 

A76 

extrude 

477 

EXTRUDE 

476 

EXTRUDE 

479 

EXTRUDE 

ABO 

EXTRUDE 

ASI 

EXTRUDE 

AB2 

extrude 

4B3 

EXTRUDE 

ABA 

EXTRUDE 

469 

EXTRUDE 

466 

extrude 

467 

EXTRUDE 

466 

EXTRUDE 

469 

EXTRUDE 

490 

EXTRUDE 

491 

extrude 

492 

EXTRUDE 

493 

EXTRUDE 

494 

EXTRUDE 

499 

extrude 

496 

extrude 

497 

EXTRUDE 

496 

EXTRUDE 

499 

EXTRUDE 

900 

EXTRUDE 

SOI 

EXTRUDE 

902 

EXTRUDE 

903 

EXTRUDE 

S04 

EXTRUDE 

90S 

EXTRUDE 

906 

EXTRUDE 

SOT 

EXTR  JOE 

906 

c 

C OOURce  cnrcn  TH€  f»<«L.T  n*TA 

C9*9**90900*9090m9000099909m0900»0m000000»000mm09m9m0m0000000m0000009000 

c 

00  9A8R  l•|•NUMNP 

IP ( CODEC  1 1 •EO«i ••ofii •CODEC  It • :o«i I ••op^coocc 1 1 .CO* 3* •o*«cooeif I 

I •EO«M«t  UQ<ll«0«0 

IP( CODEC  It •fe0*9**0e .CODEC  1 1 .EO. 1ft* I URC It  aUZC I t«TANTH 
98A8  CONTINUE 

MRITZ<6«1020t 

VRITSCfttlOlTI  CUPCII,  UZCII*  1«1»  NUMNO  | 

c 

CAPAPAPEtAEAEAAAAPPPAOAtPEPPPAVPPAAPAPAPEPPPPPPPPAPPAPPAPPPPPPPAPPPPP*** 

C ITEParOTAL  NUNREP  OP  ITeOATIONS  AftStONEO* 

C ITERNO  > AN  INOSM  PQR  OUACiCftES  SUftRCUTINE  $T|PP  AND  COMMENT  OELOVl 

C**«**AA«**AP**A«R*AAA*PRPAA*RPAAA«R*R AAA ••••*«*••*• AAAMMMAMMRAMMMMMRMMM 

C 

ITEONOa  t 

DO  2000  K|ra| •ITER 
INOCONaKir 

c 

CAA*AA«AA*AAA«AA«AAA*E««A*««T»AAAA*TTAAA«AAAAAAA*A*AA**AAAAAAAAAAAAMAAAA 

C CHECK  IP  IT  IS  FIRST  iteration 

C IP  PIRST  ITERATION.  |TERNOa| , PERPORM  OUAOl  • OTtCRWlftC  JUST  MEAD 
C RESULTS  OP  OUADl  PROM  TAPE  U 

CAAAA««A*AAA«AAAAAAAAAAA**«««»«««««««*A*»RPR««RR**R«R*«P«*RR* 

C 

S79  iPCKtr*ii  sooffso.soa 

c 

CAAAAAAAAAAAAAA*RRAAAAAAA*PR**A*PRARAR*R*AAAAAAAAAAAAAAAAAAAA 

C NAKHRO  is  an  index  TO  SUBRCUTINE  STIPP 

C NWKHROal  JUST  CALCULATE  STRAIN  RATE  IN  SUftROUTtNE  STIPP 
C MAKHRDaO  DOING  AHQLE  THING  IN  SURROUT {ME  STIPP 

CAAAAAAAAAAAAAAAAAAA*AAAAAAAAA*AA*AAAAAAAAARAAAA*»RR«««PRR«iR*AAAAAAMAAMA 

C 

999  NVKHROaO 

IPf  ITCONT.NE*  It  NAKHROal 

CALL  ST IPP  CR«ZtUR«U?*COOEt lELt TV « STS.CRS«NP,PPUR ,PST« 

I I JSC . PS8C • A • R «NEL  t NO  « NPC • NB2 • AL  ANCA * N« K MRO  t 
ITERNOa2 

IPC ITCONT.EO*! t 60  TC  593 
C 

CAAAAAAAAAAAAAAAAAAAAAAA9R«AAAAAAAAAAAAAAAAAAAAAAAA«AAAAAAAAAAAMAMMMAA«A 

C UNNECESSARY  TO  CALCULATE  TOTAL  EPPECTIVC  STRAIN  FOR  NOM^MAROeNING 
C MATERIAI.S  EXCEPT  AT  PINAL  OUTPUT 

CAAAAAAAAAAAAARRAAAAAMAARAAAAAAAAAAAAAAAAAAAAAAAAAAARRAAAAAAAAAAMMAMAAAA 

c 

CALL  STRAINSCR*  ZfUR,U2« 1EL« EPS. TEPSt ACM  1 1 * AC  M2 1 , AC  MSI « Af MAI vAIMSI # 

1 ACM6t ,ACM7t  «ACMOt . AC M9t  « A C Ml 0 t . A( Ml  1 I » ACM! 2 | * RR* AfMl 3|, NiPTSf 

2 NJPTS.NTIMSStNEL.NUMNP.NMAXt 
C 

CRARRARRAARRRRRRRRARRARRRAARRRRRARARRRRRRRRRRRARRRRRRRAARRARRRRRMMMMMRMR 
C CALCULATE  NORMALIZED  EFFECTIVE  STRESS  CPOR  HORM^HAROENJNG) 

CRRRRAARARA «**»*«•« ARRRARRRRRRRRRAARRRRRRRRRARRRRRRRRRRRRRRRRRRRRMRMR AMR 

c 

DO  220  Nai . NUMEL 
CALL  HARD! TEPSCStNt .YVCNt I 
VYCNtaVVCNt/YIELO 
220  TERSf  I .Nt  aVVf  Nt 
NRKHROaO 

CALL  STIFF  CR.ZtUR.UZf CODE. lEL.TV.STStEPSvNP.PPURtPST. 

1 I JRC.PSBCfA.a .NEL.NO.NPC.NBa.ALANOA.NVKHROt 
993  CALL  MODIFY  C CODE* A. P. NUMNP.NCO.MBANO I 
990  CONTINUE 
C 

CRAARARRARRARRARRRRRRARARRRRRRRRRRRARRRRRRARPAPRRARRRMRRAMMMMRMMMMMMMMMM 

C SOLUTION  FOR  BANOCO  SYMMETRIC  MATRIX 

CRRRRRRRAAARRAARRARRAARARRAAAARRRRARARARRRRRRARRRRPRRRRRRMMRMMMMMMMMOMMM 

c 

CALL  TRIACNEO.MBAND.At 
CALL  BACKSCNEO.MBANO. A.Bt 

c 

carprararraarrarraaaaarrraraaarrrraarrrarrrrrrrrrrrabmbmrrrmrbbmmmmbmmbb 

C SET  CORRE8SPONOING  VALUE  POP  INCLINED  SPECIAL  •OUH'IAPV  COWOlTnUt 

CR»»R»»»AR»»»AA»»»»»R»»»»*»»R*»»»»»»RR»»»»R9»PRR»»»RPP»PPP»P»»»»»BMMMP»P 

C 

DO  769  Na|.  NUMNP 
rZa3«N*| 

IR-IZ-l 

769  lP<C00eiNt*E0«5*«0R»C0DECNl*eQ«lS»t  B( IPI ■Bl t Zl ATANITNCTAI 

C 

CRRAARARRRRRRRAARRRRRRRRRRRRRRRRA6RR6RRR6R6666RRPRR66R666PPPPMMBPBMBOOOP 

C IP  PIRST  ITERATION.  SET  INITIAL  STEP  LENGTH  ACCOPOING  TO 
C OECCELARATION  COEPP*  ASSIGNED.  OTMEPBISE  PINO  TNE  BEST  OCCCCLMAriON 
€ COEPP.  PROM  INPORMAYIOn  OP  PREVIOUS  ITERATION. 

CRAP*RAR*ARR*AAAAA*A0R«AR*AARA9A9AAAAAAAA99PAPPPP99APPPPPMMPPMM< 


98 


l\ 


I 


•XTHUOe 

COO 

907 

C4«0«0 

CHTMUDe 

CIO 

00  SOS  t«l vNUMNO 

CXTflUOC 

•II 

tf«9«t*l 

CXTttUOC 

Sit 

|H«II-1 

cx'muot 

C13 

SOS 

C4ac««ai izi«ai 1 zi»ai laiaai iri 

CXTHUDC 

914 

COHNOM«SOHTIC4I 

eSTMUDC 

CIS 

lP(KfT*CQ«l 1 00  TO  SSI 

eRTiiuoe 

SIC 

47 

ACOCP«I*0 

tXTIlUOC 

SIT 

910 

IPt CACOCPPCOHNOM|*LT*SrePLM|  00  TO  9S0 

cxtKuoe 

SIS 

saa 

AC0CP«0*S4AC0CP 

cirmioe 

SIS 

SO  TO  910 

fITWUOC 

SCO 

991 

AC0CP«ACP1M1 

CXTilUOC 

SCI 

STCPLMaACOCPSCOHNOM 

ixTuuoe 

set 

STCPOHsSTCPLM 

SXYWUDK 

Sts 

990 

PPaPSOHl 

KXTmiOff 

st4 

PPPPaPSOHS 

exT»u»c 

Sts 

c 

cxTiiuoe 

sac 

tXTHUOC 

SCT 

c 

CALCUWATC  HORCCS  M NOOAL  POINTS 

txTuuoe 

Sts 

CS4SS4«4444S4444«4444««44*«4044444*4C4444S44*4S«4«ttS4*444444444SSCSCC4«4 

exTiiuoc 

Sts 

C 

IXTMUDC 

S90 

CALL  CPaHCC<NPtPHtPt*PST«PPUH*0*M9AN0*NSP*NHPai 

cxTuuoe 

991 

HHITCICtlCOTI 

fxTmioe 

SS3 

SHITClCtlSICI  KIT 

exTMuee 

993 

e 

exr«uoe 

S94 

txmuDC 

S3S 

c 

STSlCfNIaMCAN  STHCSS  AT  PHCSCNT  SOCUTION* 

tXTMUOC 

93C 

e 

M.*NOAII«l*MCAM  STIieS*  PtIOM  MCVtCUS  SOLUTION.  IP  F inST  tTCMATION, 

■xrmioe 

937 

c 

TMCSC  VALUCS  AHC  THC  SAMC  AS  STS<C*NI 

KxnHiec 

93S 

CSC444S444SS44C44444444H4 44S444444444444444444S4444S44C4SC44C444HCC444S4 

CXTttUDC 

93S 

C 

cxthuoc 

940 

OO  133  N«l*  NUMCL 

CXTAUOC 

941 

laMMOl  ICLIN,  II,  irLIN.  tl,  ICLIN.SI  ,IEL<N.4II 

txmuoc 

94t 

tH«3«l 

tXTNUOC 

•«3 

111 

9TS<C*NI«0IIHI 

tXTAUDC 

944 

|P(KtT«CO«ll  ALANOAC IHlsOtlHI 

txTfluoe 

949 

OVCCTHI IH»«9I|H|-ALAN0a<|HI 

CXTNUOe 

S4C 

133 

Of tH|«0* 

tXTRUDC 

947 

c 

tXTIHlOC 

94S 

CSS4«S44444S«4C«S4*4444««««4«H««**0H*04*4HHSHHC0C*444444*«4HS«C444C«4H44 

CXTHUDt 

94S 

C 

CALCULATION  OP  SOLUTION  NOHM  POP  VCLOCITICS 

ixTiluoe 

990 

c 

SMOHHa  NOHM  OP  SOLUTION  VCCTOH  OP  VCLOCITICS 

CXTHUOC 

991 

c 

CNOHMb  NOHM  OP  CHHOH  VCCTOH  OP  VCLOCITICS 

CXTilUOC 

S9t 

CSS»»»»*»>S4S»S»44444»4—»»4>H»44»S»>4»44»44S»4S»44H4»H»4»»4S4SS<>9»— 4 

CXTilUOC 

993 

c 

CXTHUOC 

994 

SNOHM  • 0« 

CXTHUOC 

999 

OO  134  lal*  NUMNP 

CXTHUOC 

994 

tt«3Hf-l 

CXTHUOC 

997 

IH«IZ-I 

CXTHUOC 

S9S 

INaftTl 

CXTHUOC 

9SS 

OveCTHf lH|«Cf IHl 

CXTHUOC 

SCO 

OVCCTHI ItlvSflZI 

CXTHUOC 

SCI 

SNOHM  a SNOHM  ♦ UHfllCUHflf  ♦ UZClIAUZfll 

CXTHUOC 

sea 

c 

CXTHUOC 

9C3 

CXTHUOC 

SC4 

c 

MootPiv  tnc  vcLOcrrv  piclos  ano  mcan  sthcsscs 

CXTHUOC 

9C9 

CXTHUOC 

see 

c 

CXTHUOC 

947 

ALANOAI iNlaALANCAl iNI^OVCCTHf INIPACOCP 

CXTHUOC 

9CS 

UHf 1 laUHC 1 »♦9C IHIPACOCP 

CXTHUOC 

9CS 

UZC llaUZf 1I40IIZI4AC0CP 

CXTHUOC 

970 

134 

CONTI  NUC 

CXTHUOC 

971 

CNOHMaCOHNOM 

CXTHUOC 

971 

SNOHM  a SOHTISNOHMI 

CXTHUOC 

979 

C SNOHM aCNORM/SNOHM 

CXTHUOC 

974 

e 

CXTHUOC 

979 

c 

CXTHUOC 

974 

c 

PHtNT  SOLUTION  NOHM  OP  VCLOCITICS*  VCLOCITV  OtSTHISUTIONS  ANO  NOOAL- 

CXTHUOC 

S77 

c 

POINT  POHCCS 

CXTHUOC 

STS 

c 

CXTHUOC 

STS 

c 

CXTHUOC 

9S0 

IPINOP  •LC*  01  00  TO  129 

CXTHUOC 

SSI 

00  123  tafv  NOP 

CXTHUOC 

sst 

PHClIaPHCtlHVlCLO 

CXTHUOC 

993 

PZCtlaPZf flSVlCLO 

CXTHUOC 

994 

Its 

CONTINUe 

CXTHUOC 

SS9 

las 

CONTI NUC 

CXTHUOC 

9SC 

VHITCf 4*10191  SNOHM ,eNOHM*e SNOHM 

CXTHUOC 

SST 

«H|Tef4*l004l  KtT«ACQeP 

CXTHUOC 

sss 

00  499  la|*  NUMNP 

CXTHUOC 

9SS 

tZa3«I-l 

CXTHUOC 

9S0 

iHalZ-l 

CXTHUOC 

SSI 

■HITCf 4*10021  1*91 IH|,aCIZI*UHf ll*UZf t1*Rltl«Zf 11 

CXTHUOC 

sst 

430 

CONTI NUC 

CXTHUOC 

SS3 

VHlTC<4*f0t0l 

CXTHUOC 

SS4 

00  140  |al*  N8P 

CXTHUOC 

SSS 

SMCAHaPZf IIPCOSITHCTAI  ♦ PHI  I l«S INITHCT A> 

CXTHUOC 

ssc 

140 

•HtTCI4*10l2l  NPCtl,PR|tl*PZII)  * SHCAH 

99 


exTRuoe 

997 

exTRUrS 

592 

EXTRUDE 

590 

EXTRUOE 

600 

EXTRUDR 

601 

EXTRUOE 

602 

EXTRUOE 

603 

EXTRUOE 

60A 

EXTRUDE 

eoE 

extrude 

606 

EXTRUDE 

607 

EXTRUOE 

605 

EXTRUDE 

609 

EXTRUDE 

610 

EXTRUDE 

611 

EXTRUDE 

612 

EXTRUDE 

613 

EXTRUDE 

614 

exTRUOE 

615 

EXTRUOE 

616 

EXTRUOE 

61? 

EXTRUOE 

616 

extrude 

619 

extrude 

620 

extrude 

621 

EXTRUDE 

622 

EXTRUOE 

623 

extrude 

626 

EXTRUDE 

625 

extrude 

626 

extrude 

62T 

EXTRUDE 

626 

EXTRUDE 

629 

EXTRUDE 

630 

EXTRUOE 

631 

EXTRUOE 

632 

EXTRUDE 

633 

EXTRUOE 

634 

EXTRUOE 

635 

EXTRUOE 

636 

EXTRUDE 

63T 

extrude 

636 

EXTRUDE 

639 

EXTRUDE 

640 

EXTRUOE 

641 

EXTRUDE 

642 

EXTRUDE 

643 

EXTRUDE 

644 

EXTRUDE 

645 

EXTRUDE 

646 

EXTRUOE 

647 

EXTRUDE 

645 

extrude 

649 

EXTRUDE 

650 

EXTRUDE 

651 

EXTRUDE 

652 

EXTRUDE 

653 

EXTRUDE 

684 

EXTRUDE 

655 

EXTRUOE 

656 

EXTRUDE 

657 

EXTRUDE 

685 

EXTRUDE 

699 

EXTRUDE 

660 

EXTRUOE 

661 

EXTRUOE 

662 

EXTRUDE 

663 

EXTRUDE 

664 

EXTRUDE 

665 

EXTRUDE 

666 

EXTRUDE 

657 

EXTRUDE 

666 

EXTRUDE 

669 

eytruoe 

670 

EXTRUDE 

67 1 

E XTRUDE 

6T2 

EXTRUDE 

673 

EXTRUDE 

674 

EXTRUDE 

278 

EXTRUOE 

676 

EXTRUDE 

677 

EXTRUDE 

6*»F 

EXTRUDE 

679 

EXTRUDE 

660 

txTRtlDE 

651 

EXTRUDE 

582 

c 

c CONSTRUCT  TH«  ^LOll  tlKPS  «NO  CAQ6V  OUT  THC  INTCG^ATfOM  OF  THC  STAAIM 
C OATfS  TO  OerejlMlNE  the  total  ST»A1N  OXSTPIBUTION  FOA  WOAK-HAROCNIMC 
C VATFAIAL*  THt&  STSP  IS  PFSFORAEO  CNLT  AT  FINAL  ITEPATXON  FOR  NOH* 

C ^nPKHAQOENtNS  MATERIAL 

c 

843  IP<E$NORM«LE*FLXMtT*OR«K|T«SO«lTSA|  CG  TO  56 
IFf tTCONT*EO*ll  GC  TO  95 
86  tPf ITCONT.SO* II  GO  TO  56 
N«<»*RO«l 

CALL  STIFF  fF*Z»UP tUZtCCOE* 1C L* TV « ST$, EP$i NF,FPUR •FST t 
1 1 J9C*PS9C«A,e,NEL«N0«NP€«Ne2«ALAN0A«N«KHR0l 
SF  CALL  STRAINSCR#  ZvUfi.UZf lEL* £FS» TFPS* AIM | » « A f M?> , Af M3I « AC M4I ^ACMS) • 
] Af M6I  *A(M7I  »ACMei «ACM9I fA(M|OI*A<Mllt  » A C M 1 2 1 • RR« A C Ml  3 1 IPT9» 

2 NJPTS«NT tME$«NSL«NUMNPtNMAXI 
1F( ITCONT.EO*!}  GC  TC  95 
OD  230  NaltNUMEL 
CALL  HAR0CTEPSf9«Nl*VVfN| I 
VV(NI»VVCNl /VXELO 
C 


C *#*»F***AF*»*»»»P6P*****P*******AP»P*6***»#PP**AP»»*PPP*P*P*6RPRPRRPR* 

C TEPSfWNI«  PREVIOUS  EFFECTIVE  STRESS 

C TEPSf2*NI»r>|FFEREACC  OF  EFFCT  XVE  STRESS  BETWEEN  THE  PREVIOUS  AND  NEW 
C SOLUTIONS* f FOR  WORF^HAROENl NG  MATERIALS  I 

C •RA^RAVPRvAFWPAAFVVAPRAARPPAvAPPAAPPPMPPRRPvFPARPVPPPPFPPRPRRPRRPRPRR* 

c 

TEPSC  2#N|*VY<NI-TEPSC 1*NI 
230  TEPS<l*M«TV<M 

c 

C PRINT  THE  STRESS  AND  STRAIN^RATE  OXSTRIFUT IONS  AND  THE  EFFECTIVE 
C STRAIN  DISTRIBUTION* 

C A»««A«F»««*F«pF«#***«PVPF#FR««PFP6FPFA6»***PPPF*«PPF«P«FP4r*RR*RPPRPPpP 


55  CONTINUE 

RRIT£(6«IOOTI 
«RtTE<6tl009)  KIT 
on  222  N«l«  NUMEL 
C0NSTa2*MVV(N|/<3*AEPSC8tMI > 

SrS(t«N>«EPSl t«N>PCONST 
STS(2*Nl«EPS(?«KI*C0NSr 
ST SC  3«NI«5PS( 3«KI*C0NST 
SrSCA*N>«£PSf 6, NIFCONST/2* 
no  132  J«1  * 3 

132  STSC J*N>«STSC J,NI^STSC6«NI 

575f5«N>«£FS7RSfSrSf ] #N >* STS C 2* N I *STSf  3»N I * STSC 4*N 11 
WRX TEC6«I004I  N«  C EPSC I • Nl * 1 *1 «5 I *TEPS <5,H I , ( STSC I tN I * 1« 1*6  I 
222  CONTINUE 
C 

c R*****F4*4A**4**»*****»**A44AA4*»FA**A**»PP4*PFPFFAA*»6PPP**R»RRPRRRRR 

C CHECK  THE  CONVERGENCE  OF  NE«  SOLUTION,  IF  IT  IS  NOT,  THEN  REDUCE 
C ACCEL*  COCFF*  and  corresponding  STEP  LENGTH  ANO  CHECK  IT  AGAIN 
C FOR  FIRST  ITERATION,  THIS  STEP  IS  OMITTED  AND  JUST  COMPUTE  THE  TOTAL 
C SQUARE  NORM  OF  EQUATIONS  ANO  I NCOMPPESS I BIL ITT  CONDITIONS 

C A4«ai»FP*AA»**6V«*  A*44*4F4*«4««F«*AR*4FF»44*«44*4P«P*MP***«***PR6*P*RF6 

C 

564  NWKHROaO 

CALL  stiff  (R,2»UP»U?,C00E,IELtVT,STS*EPS«NF,FPUR,FST« 

1 1 J5C«PSeC,A,e,NEL*N0,NPC,Nn2, alanda«nvkhroi 
CALL  MODIFY  I CODE • A ,6 ,NUMNP,NEQ •MBANOI 
C 

c *ft4P«*R*«***««*«4«4««*4R«»FF«644F***P4«4*«*«*F*«F»FFF«A«««P*44466R«664 

C COMPUTE  TOTAL  SQUARE  NORM  QF  EOUAT |ONSCF$OP 1 ) ANO  SQUARE  NORM  OF 
C INCOMPRESSIBILITY  CONOfTtCNSCFSORSI 

C «4R«R*R*«**4«4««*4PR*44«R««i*««**«4«**P*4:«ARR**»RR««P*RR4R*4R*R6«*PRRP6 

C 

FSOR5«0*0 
FSQRta0«0 
00  501  N«1 ,NUMNP 
1Ir2MN 

FSOClaFSORlABC 1 1 >R6C  XX  I 

FSQRSaFSORSTBCXIIRFCIt I 

DO  501  JJa2,a 

IXalX-l 

BLANf  faBC I I I 

JAaJJ 

JKatt*- JA*1 

655  IFf JA*GT*MBAN0I  GO  TO  669 

5LAN1  |a  plant  WAC  I I • J«  I A ^LANDAt  JK  I 
JAaJF^S 
JKa  JKAS 
GO  TO  655 

550  |FCN*GT*1I  GO  TO  TQ2 
GO  TO  591 


100 


•MTqgor 

683 

792  00  893  J«3»ll t3 

exTPuoe 

684 

Eifmuot 

688 

IRf XK.GT.MBANOI  go  to  891 

CXTPUDE 

686 

893  BLANl laRLANl X*A( J«XK)«ALAN04( Jl 

EXTOOTE 

68? 

891  8S0R1 •RSORl ♦etANI IRELANf I 

EXT0JOE 

698 

E0|  continue 

EXTQUOE 

689 

I8CKlT*EO«ll  GO  TO  788 

EXTOUC2 

6R0 

C 

ExTAUnE 

6Rl 

c 

EXTOOOE 

69? 

c 

CCMPAR3  SQUARE  NORR  08  RtNCT|ON$  RlTM  PREVIOUS  ONE* 

£xTQUr>E 

693 

c 

'AAAP9AA99AAA 

EXTAUHE 

694 

c 

EXTAUOE 

69* 

Of  8RsRR.*50A| 

EXTRUDE 

696 

I8tFS0R|.FR|  t?;, 871,672 

extrude 

697 

c 

EXTRUDE 

698 

c 

A •«««•«•*••*•  »»**••«•  ••4«R*««RRRR*ARA«A«R«ARAAR«««RRAtR4RR4Q«49»RRR 

EXTRUDE 

690 

c 

CASS  FOR  SQUARE  NORN  OF  NER  SOLUTION  LARGER  THAN  PREVIOUS 

ONE,  REOUCE 

EXTRUDE 

TOO 

c 

ACCEL*  COEFF,  RV  half  AKO  CCRRARE  AGAIN* 

EXTRUDE 

?0I 

c 

'AAAP«PA«A9AA 

EXTRUDE 

70? 

c 

EXTRUDE 

T03 

572  IF( CARSlDtFFI /FFI .LE*! *C-3l  GC  TO  44A 

EXTRUDE 

▼04 

STERLH«0*5FSTEPLM 

EXTRUDE 

TOS 

ACCCF«AC0CF«0«5 

EXTRUDE 

706 

1NOCON«|NOCONF1 

EXTRUDE 

TOT 

1F<  tNOCCN-< KIT 451 168 1*681 *59? 

EXTRUDE 

Toe 

591  DO  595  I«1*NUMNR 

EXTRUDE 

709 

IZ»3R 1-1 

EXTRUDE 

710 

1R«IZ-I 

EXTRUDE 

711 

IN-TZ^I 

EXTRUDE 

71? 

ALANOAC iNlaALANOAl IN l-GVcCTRC IN  IRACOEF 

EXTRUDE 

713 

UR( tl «UR( 1 l-GVECTR( ICIAACOEF 

EXTRUDE 

714 

895  U7(  llaUZdl-GVECTRCIZlAACOEF 

EXTRUDE 

71S 

FR1TC(6«11SS1  STERLH*r50Rl«FF*0|FF*FS0RS«FFFF 

EXTRUDE 

716 

llR|Tei6*1030l  ACOEF 

EXTRUDE 

TIT 

DC  772  I«i*nUMNP 

EXTRUDE 

718 

lZa3«I-l 

EXTRUDE 

719 

|R*tZ-l 

EX^RVAE 

720 

«eiTE(6*l0021I*B(  IP  1,8  ( 1Z1*URM1*UZC  II 

EXTRUDE 

T?1 

772  CONTINUE 

EXTRUDE 

??? 

IKIITCONT.EO.il  GO  TO  964 

EXTRUDE 

7?3 

DO  861  N«1 ,NUREL 

EXTRUDE 

7?4 

861  VV<NlaTEPS< l*Nt-T£R5f 2*N>*0*S 

E XTRUDE 

T25 

GO  TO  864 

EXTRUDE 

726 

c 

EXTRUDE 

??▼ 

c 

EXTRUDE 

728 

c 

IF  CONTIMUOtS  REDUCING  CECCEL*  CCEFF,  OF  FIVE  T|HfcS*  TK  SOLUTION 

EXTRUDE 

72R 

c 

STILL  NOT  CONVERGENT*  TM«n  PUNCH  AND  PRINT  THE  RESULT  AND 

STOP  THE 

EXTRUDE 

730 

c 

PROGRAM* 

EXTR UOE 

731 

c 

« ««  A*  A*  ***«•«*«•****  A«««r«*«A«*AAAA*»A«AA«iA«»A«R 

EXTRUDE 

732 

c 

EXTRUDE 

733 

58?  «PIT£C6*1155I  STEPLN*F£ORi * FF, D IFF* FS0P5* FFFF 

EXTRUDE 

734 

«RtTEC6*10561 

EXTRUDE 

735 

1TST»2 

EXTRUDE 

736 

GO  TO  600 

EXTRUDE 

737 

c 

EXTRUDE 

738 

c 

extrude 

T3R 

c 

CASE  FOR  SQUARE  NORM  OF  nEA  SOLUTION  LESS  THAN  PREVXCUS  ONE*  IF 

EXTRUDE 

740 

c 

CONVERFENCE  CHARACTER  IS  IN  GOOD  SITUATICN*  THEN  INCREASE 

STEP 

EXTRUDE 

741 

c 

LENGTH  FOR  NEXT  ITERATION* 

EXTRUDE 

742 

c 

EXTRUDE 

743 

c 

EXTRUDE 

744 

871  FS0R2*I $TEROR/STERLH|A|FSOR|-FF|aFF 

EXTRUDE 

745 

1^(49042-0. 96S4KPI  574.974.979 

EXTRUDE 

746 

975  STERLH«0*8A$TCPLH 

EXTRUDE 

T4T 

GO  TO  444 

EXTRUDE 

746 

S74  |FCFSOR2-0*900AFF1  576*576*444 

Extrude 

749 

876  STERLHvUaASTEPLH 

EXTRUDE 

750 

444  continue 

EXTRUDE 

751 

STEPORmSTEPLH 

EXTRUDE 

752 

AR1TE<6* 11551  STEPLH*F$0Rl*FF*DrFF*F80RS*FFPF 

EXTRUDE 

783 

788  CONTINUE 

EXTRUDE 

754 

c 

Extrude 

788 

c 

exTR«iDE 

756 

c 

STOP  ITERATIONS  4N0  PUNCH  RESULTS  IP  (ERROR  NORRI/ ( SOLUT ICN  NORMI 

Extrude 

757 

c 

REACH  THE  ASSIGNED  ACCURACT* 

extrude 

758 

c 

PRPRRRRPRRRRR4RRTRPRRRRRR4R4RRR4RRRRRR4PPR4R4R4RRPRRR4RPR4RRPPPPPR4RPR 

EXTRUDE 

759 

c 

EXTRUDE 

750 

rF<E8N0RM*LC*FLrMrT  1 GO  TO  600 

EXTRUDE 

761 

c 

EXTRUDE 

762 

€ 

EXTRUDE 

763 

c 

PUNCH  SOLUTIONS  AT  EVERT  FIFTH  ITERATION  BUT  DOES  NC T STOP 

THE 

EXTRUDE 

764 

c 

ITERATIONS.  PUNCH  SOLUTIONS  AT  PIKAL  ITERATION. 

EXTRUDE 

768 

c 

EXTRUDE 

766 

c 

EXTRUDE 

767 

K|«K|T/5 

EXTRUDE 

768 

K2«K|99  2^01 

PT' 


EXTRUDE 

769 

FXTRUDE 

770 

EXTRUDE 

771 

EXTRUDE 

7T2 

600 

EXTRUDE 

773 

EXTRUDE 

774 

EXTRUDE 

775 

EXTRUDE 

776 

EXTRUDE 

777 

EXTRUDE 

77B 

EXTRUDE 

779 

EXTRUDE 

7B0 

100 

EXTRUDE 

781 

660 

EXTRUDE 

782 

EXTRUDE 

783 

EXTRUDE 

7R4 

2000 

EXTRUDE 

788 

c 

EXTRUDE 

786 

C •**« 

EXTRUDE 

787 

C PP 

EXTRUDE 

796 

C IN' 

EXTRUDE 

789 

C ***' 

EXTRUDE 

790 

C 

extrude 

791 

2205 

EXTRUDE 

792 

EXTRUDE 

793 

EXTRUDE 

794 

EXTRUDE 

795 

500 

extrude 

796 

EXTRUDE 

T97 

C 

EXTRUDE 

798 

1002 

extrude 

799 

1003 

EXTRUDE 

800 

1004 

EXTRUDE 

801 

1005 

EXTRUDE 

802 

EXTRUDE 

803 

EXTRUDE 

804 

EXTRUDE 

805 

EXTRUDE 

606 

1006 

EXTRUDE 

607 

E XTPUDE 

608 

EXTRUDE 

809 

EXTRUDE 

610 

EXTRUDE 

811 

EXTRUDE 

812 

1007 

EXTRUDE 

813 

1010 

EXTRUDE 

614 

1 

EXTRUDE 

915 

EXTRUDE 

616 

1012 

EXTRUDE 

817 

1015 

EXTRUDE 

816 

EXTRUDE 

819 

EXTRUDE 

820 

1016 

EXTRUDE 

821 

1017 

EXTRUDE 

922 

1019 

EXTRUDE 

823 

EXTRUDE 

824 

EXTRUDE 

826 

EXTRUDE 

826 

1020 

EXTRUDE 

827 

extrude 

828 

EXTRUDE 

829 

1030 

EXTRUDE 

830 

EXTRUDE 

631 

EXTRUDE 

832 

EXTRUDE 

833 

1 

EXTRUDE 

834 

1 

EXTRUDE 

835 

EXTRUDE 

836 

EXTRUDE 

837 

1056 

EXTRUDE 

838 

1155 

EXTRUDE 

839 

1 

EXTRUDE 

840 

EXTRUDE 

641 

1 

EXTRUDE 

842 

EXTRUDE 

943 

EXTRUDE 

844 

1 

EXTRUDE 

84S 

EXTRUDE 

846 

2222 

EXTRUDE 

847 

C 

EXTRUDE 

848 

EXTRUDE 

849 

IF|K2«eir«KtT)  GO  TO  600 
t^(K1T«E0«tTeP I GO  TO  600 
GO  TO  660 


NN«36MiPNP 

1F<NOlP4CM  •EG*  01  GO  TO  660 

OUNCH  |0I7«  (TEPSf StM)»Nsl«NUMELl 

PUNCH  ioi7«  ( <epsc itNitivitSi*  n«i,numcli 

PUNCH  |0lTt<<STSII«N),I«l«6)  •N^l  .NUMEO 

OO  too  IsttNBP 

t(«NF<  J > 


|P( tTST.EO.2)  GO  TO  2205 
IFCCSN0PM«C£«PLIMIT t 60  TO  2206 


WPITEf 6»t0l9»  KIT 
00  600  K»I»NUMEL 

SUMvEPSf  t tN|4EPS(2tN)^EPSf 3«NI 


«PlT6r6«2222>  STEPLH 


10X«  2P13.7I 


FORPAT<tHt««  STRAIN  RATE-'STRESS  SOLUTION  AT  fTER«  NUMR6R  ■P»f4// 

1 « EL*  N0**R*STRAIN*«  2-STPAXK*T»’-STRAtN*RZ-STRAlN*EF-STRAlN«TOT*E 
2FF....R-ST  ReSS**Z*STPESS*TffoSTR£SS*PZ-5TRES5*EF-STRE  SS«AVG-STRESS« 
3**«/* 1 1 XfPRATE* •6X*«RATE**6X»PRATE6«6X ••RATEP«6X*6RATEP*8X««STRAtN 
46/1 

FORPATf///  30X«  * VELCCITV  SOLUTION  AT  ITERATION  NUMBER  aPtlA 

l/«30X,«  06CELARAT10N  COSFFl CIENTs*  Pt3*6* 

1///  15X«*CORRECTEO  AMOUNT  OF  VEL*Ft t6X««NE«  SOLUTION4,2l XtPCOCROlN 
1ATES«/«24X ,*1N* •»«29X**OF* 

2 /*  nodal  NO*  R-VELOCITt  Z-VELOCITV  R-VELOCITV 

3 Z-VELOCITY  R-COORO  Z*C00R04» 


iMAT<  //  4 nodal  point  FORCE*// 

Z-FORCE* 

2 surface****) 


*5HEAR  FORCE  ON  DIE 


50X, 

60X« 

60X* 


NORM  OF  SOLUTION  VECTOR  •**  Ft3*a 
NORM  OF  ERROP  vector  ■« • F|3«6 

fractional  norm  «*«  F13*6) 

VELOCITY  SOLUTION  AT  ITERATION  NUMBER  •*«14l 


1  * EL*  no***r*strain«*«z-strain**th»strain**pz-strain**ef-strain** 

2*****  TOT^EF^STRAIN**/*  t2X»*RATE**7X**RATE*7X»  *RATE* 

3,7X**RATE*t7X**RATE*/) 


lTtON*t/«*  FROM  FIRST  NODAL  POI NT (R*V6L* • Z-VEL*I  TO  LAST  NQOAL 
2 P9INT(P-VEL**  Z*VEL«I*I 


1/f 2SX«*REOUCS  OECCELARATION  COEFF*  RiTH  THE  FOLLOMING  VELOCITY  CIS 
tTPtBUTtONS*/*26Xt«ANO  CALCULATE  ALL  NORMS  AGAIN* 

I/«25X«*NEV  06CCELARAT1CN  COEFF I C I ENTF**Ft 3«B« 

1 //  t5X**CORRE6TEO  AMOUNT  OF  VEL*»« tOX«*NEV  SOLUrlON*/*24X«*tN* •*« 
l29X«*0F**/« 

2 * NODAL  NO*  R*VEL0C1TV  Z*VELOCtTY  R^YBLOCITY 

3 Z-VSLOCITV*) 


) /•26X«*T0TAL  SQUARE  NORM  OF  ALL  EQUATIONS  AT  PRESENT  SOL 

lUTtON»«EI6«8t 

I /t2SX**T0TAL  SQUARE  NORM  OF  ALL  EQUATIONS  FROM  PREVIOUS 

1S0LUTI0N**£16*8t 
I /«25x««CTFFER6NCE«*EI6*Bt 

I //«2SX**S0UARE  NCRM  OF  PRESENT  VALUE  IN  VOL *CONS |STANCEF*Bl6*Bt 
t /*28X««SOUARE  NORM  OF  PREVIOUS  VALUE  IN  VOL*CONS 1ST ANCE«*EtS*B I 


RETURN 

END 


102 


\ i 


I 


Y 


IITTRUOC 

ESI 

soanourtNE  HAaofSTaAiM.vittDi 

CMTHUOC 

•fit 

c 

fXTVUCt 

• fiS 

c 

EXTftUO! 

•S4 

c 

saeciar  the  matcaial  pAo^eArirs  »ATeAfA4,  ftraess* 

tKtttMOf 

• Sfl 

c 

STB*  IN  BXLBTIONIHIPt  *Ot>  N0N-M*n0eNIN6  MATEItI*k>  BCPtACC  •VltLM 

CARD 

exTiluoe 

c 

ANO  set  VIECOsTHC  constant  VIEL9  STRESS* 

EMTQ^e 

•ST 

f 

»R««RR 

e«Ti»upe 

asA 

c 

EKTffuoe 

•so 

EVTiiuOE 

•so 

AN«0*  OS 

EKT41POE 

••1 

VlELOvI 1 •♦STRAIN/B)R*AN 

f xTaoOE 

• 62 

RETURN 

CXT9UOE 

•as 

ENO 

e«Toiir>S 

••S 

SURROUTlNE  STSER  < R • 2 tUR *U2 *COCE • 1 EL*?? tSTS tERS •NE«E8UR*ESTt 

EXTRUOE 

•as 

1 12BC*RS8C*A*8*NeL*N0,NRC*N8  2*  AI.AND4*N«RHR0| 

EXTRUOC 

•67 

C 

EVTRUOE 

•68 

c 

EXTRUDE 

•6R 

c 

CALCUtATION  STIERNESS  RATRlX  FOR  ENTIRE  STSTfiR* 

ExTRure 

•TO 

c 

9999FR 

EXTRUDE 

• 71 

c 

EXTRUDE 

672 

C0«M0N/6eNC0N/NUNNP,NUPfCtHEDI 1 2t iVCL<Neo.NS« ITERNO« ISTOP t 

EXTRUDE 

•T3 

IVIEUD.MP*SO»TEST,WDTAGtNRF,NI.'«BC»ITSP«NRP2,NSC*LE,NPUNCH,KP«|NT. 

extrude 

• 74 

2 NCHECK.ACPINI 

extrude 

•75 

C0NR3N  /XALL/  THET4*Rt,t4NTH 

ExTmjDE 

876 

CQNRON  /STRRAT/  Rr9*9l«M(9l 

EXTRUDE 

877 

DIMENSION  4t4ND4CII 

EXTRUDE 

• 78 

DIMENSION  R(  11*21 1>  •CODE!  MtURf  1 1 *U7I1  1 tlELC  NEL  •ll*8(l>*Af 

NO*  11 

EXTRUDE 

8TR 

l.EBSIS«ll*S?SI6*tl»NFI||,FSTf  N82  * 1 1«ERU* I 1 1 * 1 JBCf  NRC  *11* 

EXTRUDE 

S80 

2 PSSCfNRC* t > »RR(4I*22(6) •1XVC9I •UUIOI*TT( j I • 10121 

EXTRUDE 

• 81 

c 

E XTRUDE 

882 

c 

EXTRUDE 

8*3 

c 

lNlTt4Ll2E0  A AND  6 MATRICES  IFOR  EOUATION  AK»8t 

extrude 

•84 

c 

•ECAUSfi  BANDED  STMMeT»tC  PROPERTY  QX  THE  »T|FrNESS  MATHIX  A*  THE 

EXTRUDE 

•as 

c 

STORAGE  OF  A 18  IN  A SQURE  ARRAY* 

EXTRUDE 

886 

c 

IR9RR9 

r XTOUns 

887 

c 

extrude 

888 

NSTORsO 

extrude 

89R 

00  50  N«|*  NSC 

EXTR UOC 

aRO 

AINlftO* 

EXTRUDE 

891 

DO  fiO  M«|«  NBANO 

extrude 

•R2 

50  ACN*M1«0* 

EXTRUDE 

693 

c 

EXTRUDE 

804 

c 

»*RRR4 

EXTRUDE 

aos 

c 

CONSTRUCT  ELEMENT  LEVEL  MATRICES  FOR  ST|RRNCSS 

EXTRUDE 

6R6 

c 

fP  ITEPNO*Or*l  MEANS  SURROUTINE  OLADI  HAS  SEEN  PfRROMEO  ONCfi  8CF0RK 

EXTRUDE 

8RT 

c 

AND  IS  STORED  IN  TARE  1* 

EXTRUDE 

898 

c 

•••949 

EXTRUDE 

899 

c 

E XT9UDE 

900 

IFttTERNO.GT.il  GC  TO  20 

EXTRUDE 

901 

REMIND  1 

EXTRUDE 

202 

DO  SIO  N«1«NUMEL 

EXTRUDE 

903 

00  100  1«t*4 

EXTRUDE 

904 

ii«iCL(N,n 

EXTRUDE 

90E 

RRCii*Rf rt I 

EXTRUDE 

906 

100  221  t»s2C  tn 

extrude 

90? 

CALL  OUADI tAR*22*VOLl 

EXTRUDE 

908 

A?  IFC VQL*GT* 0. 1 GO  TO  510 

EXTRUDE 

909 

RRtTEtO. lOOSI  N 

EXTRUDE 

910 

NSTCRal 

EXTRUDE 

911 

510  continue 

e XTRUDE 

912 

20  REWIND  1 

EXTRUDE 

913 

IFINSTOP  .80*  11  STD® 

extrude 

914 

c 

•XTRUDE 

91S 

c 

•••••• 

EXTRUDE 

916 

c 

CONSTRUCT  R AND  N MATRICES  At  ELEMENT  LEVEL* 

E xTRuce 

917 

c 

extrude 

918 

c 

EXTRUDE 

919 

00  1000  N«lt  NUMEL 

EXTRUDE 

920 

DO  60  I«l*4 

EXTRUDE 

921 

n«13LfN*ll 

EXTRUDE 

922 

X2«2Rt 

EXTRUDE 

923 

llalE-l 

extrude 

924 

fXYf 

EXTRUDE 

925 

IXV<I2>*3»II-1 

EXTRUDE 

926 

UUI  12l«U2t  1 1 1 

E XTRUDE 

927 

60  UUI fl l•UR(XII 

EXTRUDE 

928 

1 laMAXOt  tELIN.l 1 • tEL(N«2l * lEL <N, 3 1 . I EL t N«4 1 1 

EXTRUDE 

929 

IXYI9 )«3R| I 

EXTRUDE 

930 

CALL  0UAC2lUU«|RLKAX,TEx.TeY tTE2*Texvi 

103 


CXTQUOr 

RSI 

C 

eXTBUOC 

9S2 

c 

RRRi 

*6*9*9699«**9 

fKTnupe 

933 

c 

NKKMtD  IS  *N  INOex  TO  INOICATC  UNCTMe*  ONCV  ST**IN  •ATE  CALCULATION  It 

exT»uoc 

939 

c 

nEOUIEBC  OA  THE  ••-ole  STITENESS  HATBIX  need  to  EC  CONETOUCTEO. 

tXTBUOC 

939 

c 

R«6< 

A***** ****••«» *• ************ •«•••*****•****•**•****•** 

****** •*6**** 

exTPuoff 

936 

c 

exTPuoe 

93T 

|R( AMMMHO.CC.I t GO  TO  EOt 

exTPuof 

990 

c 

tXTPUf>E 

939 

c 

**#■ 

exrxuoe 

960 

c 

PERXOAH  THE  ASSEMSLT  OPERATION.  BECAUSE  MATRIX  A IS  SYMMETRIC.  ONL* 

EXTBUOC 

94t 

c 

UPRER  H6LP  OP  THE  MATRIX  IS  CREATED*  AND  THE  STORaOE  POR 

matrix  a is 

eXTRUDE 

962 

c 

A SQUARE  ARRAY  OECAUSE  OP  EANOEO  SYMMETRIC  PROPERTY* 

EXTRUDE 

963 

c 

•*R| 

••*«•*•***•** 

EXTRUDE 

966 

c 

EXTRUDE 

965 

DO  62  iRltO 

EXTRUDE 

966 

H(  1 laHI  t l*VVfNl 

EXTRUDE 

967 

00  62  Jal»6 

EXTRUDE 

QRR 

62 

P<  I tJIvPIl • JI*VVIN) 

EXTRUDE 

969 

DO  200  lal*  9 

EXTRUDE 

950 

If atXYfl) 

EXTRUDE 

951 

SCtMaOf  ltl*H<|) 

EXTRUDE 

992 

00  200  9 

EXTRUDE 

983 

JJ»IXT|J»-|IAl 

EXTRUDE 

956 

IPfJJ  vUT*  1)  60  TO  200 

EXTRUDE 

958 

A«I I,JJ«MA<II,JJtAR(l«J> 

EXTRUDE 

956 

200 

CONTINUE 

EXTRUDE 

987 

201 

continue 

EXTRUDE 

980 

EPSMfNIaTEX 

EXTRUDE 

959 

S»S<2tN|aTeT 

EXTRUDE 

960 

EPS<3«N)aTB2 

EXTRUDE 

961 

EPS! 6«N)«TEXY 

EXTRUDE 

962 

1000 

CONTINUE 

EXTRUDE 

963 

663 

DO  ISO  N«1,NUMEL 

EXTRUDE 

966 

150 

PPSrStNIaReARIEPSC  1 «N1  tEPSf  2*N»  tCPSf  3*N1  «EPSU  «N|  1 

EXTRUDE 

965 

|P| N«XHRO*EO«1 1 RETURN 

EXTRUDE 

966 

c 

EXTRUDE 

967 

c 

Y*6i 

*•***«*««*»•* 

EXTRUDE 

960 

c 

PREPARATION  POR  PORCE  C A4.CUI.6TION 

EXTRUDE 

969 

c 

*6*1 

EXTRUDE 

970 

c 

EXTRUDE 

971 

tPfN8P  .te*  01  60  TO  602 

EXTRUDE 

972 

M0AN02«2PMeANC-l 

EXTRUDE 

973 

00  330  f •! • NOP 

EXTRUDE 

976 

I2«2*I 

EXTRUDE 

975 

tR«IZ*l 

EXTRUDE 

97€ 

DO  930  J«l»  MBAh02 

EXTRUDE 

977 

PSTC  tR* J1«0* 

EXTRUDE 

976 

330 

PST(12*JI«0* 

EXTRUDE 

979 

DO  600  l«lt  NCP 

EXTRUDE 

900 

1 T«NP< ! 1 

EXTRUDE 

901 

17«9*lt-l 

EXTRUDE 

902 

IRaff-l 

EXTRUDE 

903 

I1Z«2*1 

EXTRUDE 

906 

TlRaf 17-1 

EXTRUDE 

908 

00  601  J«M0ANC«  M8AN02 

EXTRUDE 

900 

JJsJ^MOANO^l 

EXTRUDE 

967 

PSTC11R«J1«A< iRf JJl 

EXTRUDE 

900 

601 

PSTCI12«JI«A<tZ« JUl 

EXTRUDE 

909 

00  603  Jal*  M8AN0 

EXTRUDE 

990 

NR«IR-J*1 

EXTRUDE 

991 

N7*I7-J*I 

EXTRUDE 

992 

JJsMRANO-JYl 

EXTRUDE 

993 

IPCNR  bLE*  01  60  TO  606 

EXTRUDE 

996 

P$Tf ITRf JJ1«A(NR,J) 

EXTRUDE 

098 

60* 

IPIKZ  *Le*  01  60  TO  603 

EXTRUDE 

996 

PSTIlI7*JJI«A<NZf 

EXTRUDE 

997 

603 

CONTINUE 

EXTRUDE 

990 

PPURItf RlaailRl 

EXTRUDE 

999 

600 

PRURI llZlaRltZI 

EXTRUDE 

1 000 

602 

CONTINUE 

EXTRUDE 

tool 

c 

EXTRUDE 

1002 

e 

*«*i 

************* 

EXTRUDE 

1003 

c 

ADD 

PRESSURE  80UNC6RY  CONDITION*  fPOR  CONSTANT  PR  fCT fONAL 

STRESS  DIES* 

EXTRUDE 

1006 

c 

THE 

SHEAR  TRACrtONAL  STRESS  IPRfCTlDNAL  STRSSSI  IS  A001N6 

HEREl* 

EXTRUDE 

I 006 

c 

***j 

• 

: 

t 

EXTRUDE 

tooo 

c 

EXTRUDE 

1007 

I*  ( NUMRC  .CO.  0 I CO  TO  AIO 

extrude 

1000 

tOf 1 laz 

EXTRUDE 

1 009 

IDf  21*1 

EXTRUDE 

1010 

Rt«l*0 

exTRUD* 

1011 

RJal*0 

EXTRUDE 

1012 

DO  420  LaUNUMPC 

EXTRUDE 

1013 

1*1  J8C(L*I  1 

EXTRUDE 

10)6 

JpfJRCf L»21 

EXTRUDE 

1015 

OR«R(tt-R(JI 

extrude 

1016 

02«Zlll-IIJt 

i 


eiTauoe 

ev^uDc 

ektruoc 

evTRuoe 

ruTRue? 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

extrude 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 


EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 


r 


tOlT 

toie 

lOtR 

lOEO 

lOEl 

lOEE 

10E3 

lOER 

tOES 

lOEt 

1027 

lOES 

10E9 

1030 

1U31 

1032 

1033 

1034 

1039 

1036 

1037 
1030 
1034 

1040 

1041 
104E 


DC  E40  RRltE 
N«I0<R| 

IxlJOCILtRl 

J«1 JRC<L«Nl 

RtaRCn 

RJaXCJl 

tE«X4t-l 

ll«tE«l 

Pl«RSRClCtM» 

PJ«PSRC<L«NI 

SIxPSBCfLfM^EI 

SJ-PS0C<ttN^2» 

PMsfRIRC  3»0RPl4Rj|»SJR<PtRPJI l/tE«0 
SM« <R 1R< 3*0  49 1XSJ  14RJ4<SI4SJII/1E*0 

•1 aDERPtfRORRSM 

Rj^a«OR4Pftl»0249M 
0(tll«efII»4Rt 
0CI2I  aRIIEIRRE 
240  CONTINUE 
420  continue 
410  CONTINUE 


tOOS  PORNATI//  24H  ELfiNENT  XtTM  NEGATIVE  AREA  •*  ID) 


RETURN 

END 


1044 

1045 

1046 

1047 
1040 
1044 
lOSO 

1091 

1092 
1 093 
1094 
1099 

1096 

1097 

1090 
1094 
1060 
1061 
1062 

1063 

1064 

1069 
1066 
1067 
1040 
I 064 

1070 

1071 

1072 

1073 

1074 
1079 

1076 

1077 
1070 
I07R 
1060 
lOEl 

loot 

1063 

1004 
1069 
1094 
1067 
1009 
I 0r4 
1040 

1091 

1042 

1043 
1 144 
1246 
I 006 


SURROUTINS  0UADlfRR»E2,V0tl 


CONSTRUCT  9 MATRIX  AND  X MATRIX  iKaC  TPANPOSEI 9)  I COM  8t  1 • AND  Q VECTOR 
AT  ELEMENT  LEVEL  POR  ALL  ELENENTS  AND  WRITE  THESE  NATRtX  ON  TAPE  t* 


COMMON  /CUAO/  6«4«9)*XXI9*9  )*Cf 9) 

OTMENStON  RRI41*22|4),SSC4I«TTC4» 

DATA  SS7«1 •»4l««4l«»*l«/«  TT/*1*«-1«*41**41*/ 
S0RT2«1*41 421 3962373042 
00  1 t«l*9 
OC  1 laO* 

DO  1 4«) *4 
91 J.tlaO* 

R1«RRC II 
REaRRCEl 
R3aRR|3) 

RAaRRIAl 

E1«72C1 I 

r2«Z2(2l 

23x22131 

24a2244l 

Rl2xRl«R8 

R13xRI«43 

R|4aRl«R4 

R23xR».R3 

R66 aRt-PA 

R34XR3-R4 

2l2x21«22 

213x21-23 

214x21-24 

223x22*23 

224x29-24 

234x23-24 

VOLxRI  3R224-R24921  3 
|R|^«Lt*0*l  RETURN 


RRRRRRRR»R>RRR4»»R4RRRRR4R4»RRR4»XR»XARW»»»R»»RRR»»MA>»»»>M9R60»»RR»WW 

CALCULATION  OP  TOTAL  VOLUME  PEA  RADIAN  OP  THE  ELEMENT* 

R*R*RR4RR44444R4RR44R44R4R4R494R*RR4R9«R4R4RRRR«RRR44RXRRR«RR606RR0RRR 


V0L1x|2|-24|R|R14A1aR9«R44R19R4| 
V0L2x| 23-1214 IR34R34R2RR24R3RR2 I 
VOL  3x129-21 lOIRlAR I ARtR* EAR IRRtt 
V0L4xf24-23)4|R3RR3AR44R4AR3RR4l 
VOL Lx I VOL 1 A VOL 24 VOL 34 VCL4 1 / 9 * 


R4RARRRRRRRRRRRRRR44R4RRXRRRX4RRRR4RRR 

THE  P0LLOWIN9  CONSTRUCT  9 AND  K NATRIX 
RQINTt  e«  the  BLEMENT* 

RRR4«444RRXR44R0RRRRRRRRWRRRRRRR44RRRR4i« 


XRRR4RRRRRRRR0RRRRRR96R44 

TAWIND  POUR  INTtSRATION 


! < 


iRRRRRRRRRRRRRR—O 0ORO69D096DDD4 


105 


I 


CVTDUOE 

evrituoe 
eiTHUDC 
EnTMuee 
CXT^UPE 
f VTRilCE 
CHT^UOS 
CITSUDS 
exTttuoe 
evT«uf>e 
FXtVUDC 
evTvuoc 
exTituoe 
evTxuos 
exrmioE 
evTRgoe 
exTQUDC 
EXTttUOE 
exTffuoe 
EXTRUDE 
EXTRUDE 
extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRI/>E 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 
EXTRUDE 
EXTRUDE 
E XTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 

extrude 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 


DO  IE  tX*l«4 

S«SSf III*0«S77390aeRI6E6a6 
TaTTIt  t )R0«57T350E69 169426 

XJAaVnL^SvCRaRAZta-RiaRZaAI^TRfRXSRZlR-^RIRRtXSI 

XJaXJE/A. 

SP«|«^S 
TMal.-T 
TP* I • ♦T 

H1«0.25*SMRPTP 

MTaO«2SPSP«TM 

h3«0«29PSPPTP 

H4»0«2S«$MM*TP 

ff«  HI  PR  I ♦H2PR24H3PP  3PM4PR4 

Gl«Ht /R 

62«H2/R 

63«H3/R 

G4«H4/R 

PPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPRPRRRPRRP 

CONSTRUCT  R MATRIX  AT  EACH  POINT 

pttPpppRRpRARpRVRR PPPPPPPPPPPPPapPPPPPPPPRPPPPPPPpRRPMPRPRRRRPpRRRRRRMR 

X|a(<-R24*R34PSPR23*TI/XJa 

X2«(R13-R34P5-RI4PTI/XJ« 

K3«<R24«P|2PSTR|4PTI/Xje 

X4«C-R13pR|2PS«R23PT)/XJ6 

V|«IZ24«Z34PS-Z23«T|/XJ9 

T2«(«ZI3PZ34PSTT14PTI/XJ8 

T3a<*Z24«Zl2P$-Z14*TI/X3e 

T4s<213-ZI2PSpZ23*TI/XJ6 

Btlf  llaVl 

Bl  l«  3)«T2 

Rfl*  S)«V3 

Bfl t 7>aV4 

Bt2«  2)«X1 

Bf2«  4«aX2 

B(2«  BI«X3 

612*  BI«X4 

a<3«  lt«6l 

Bf3«  3laG2 

eC3«  SI«G3 

RI3,  7|sG4 

RI4«  1)«Xl/$ORT2 

Bf4*  2>«V|/S0RT2 

61 4»  3l«X2/SORT2 

B<4,  4|aY2/S0RT2 

9I(A,  SlaX3/BORT2 

R<4,  6t«Y3/SORT2 

e<4*  7I«X4/S0RT2 

Bf4«  6I«T4/S0RT2 

RPpRpRRRRpppRPBRRPPRRpAPPPPpRpRMpppPRRpRftRRRRRRRR^RRRRRppppRRpPBBRRRpR 

CONSTRUCT  X MATRIX  AT  EACH  POINT  flE*  XXCItUII 

PPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPBPPPPRRRRRRP 

00  6 lalfB 
DO  6 JaZ  «B 
0UMl«0* 

00  S Xa|,4 

s OUMlaDUM|PB€K«J>PBf K«l> 

XXC  tfJIaOUMt 
XXf JtllaOUMI 
6 CONTINUE 

PPPPPPPPPPPPPPPPPPPPPPPPPPPPRPPPPPPPPPPPPRPPRRPBPPPPPPPPPRPBRPRRRRRRRR 

CONSTRUCT  0 VECTOR  AT  EACH  POINT*  AND  ADO  ALL  POUR  INTEBRATION  ROtNTS 

PPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPBPPPRPPPP 

XJRaXJPR 

OC  tlaOlllPCYlPBlIPXJPR 
0<2t«0<2lP  XI  PXJPR 

Of 3|aOf3>PfV2PB2>PX^PR 
0<4|aOI4l4  X2  PXJPR 

OI9|aOfSlPfV3PB3l4XJPP 
OfBlaOCBlP  X3  PXJPR 

OfT|aO<T|P(V4PG4IPXJPR 
OCBIaOfBlP  X4  PXJPR 


c«r«uoe  It*! 

CMTVUOC  \l%9 
rvTvuoc  HAS 
FvT«Uf>C  1 144 
eiTftgDC  It 49 
eirrituos  1196 
I1IT4UDC  1147 
9474007  1149 
9479009  1149 
9 479009  1190 
9479009  1191 
94T900C  1199 
9479009  1199 
94T900C  1194 
9479009  1199 
9479009  1194 


e •••••••••••••••••••••9469944«***«999«»494*«9«4499999999999999999999999 

e 99IT9  O 99^709  I44T9R  INT964atION«  1*C«  400tN6  4009  90tNT6l  4M0  TMK 
C r074L  VCC0M4  04  TmC  9L949Kr  OM  7449  !• 


U4I49I3969979092 

I«9t474497l99t996 


CXT4UDC 

lacA 

c 

CXTAUOC 

lacs 

MCAOl IIIEI J,4  1 

EXT4UD8 

ia«6 

00  IT  avitS 

CXTMUOC 

tCMT 

IT 

HIM  l■Mf4  lOUUIJt 

BXTAUOC 

ia«s 

MCAO  fit  VOL 

CXTMUOC 

ia«M 

DO  IS  l•l«M 

CXTMUOC 

ICT4 

00  18  JttltS 

CXTMUOC 

laTi 

IS 

E(J,|I«E<|,2I 

CXTMUOC 

laTz 

TCX*TCX/V0L 

CXTMUOC 

lars 

TCV«TCV/VOL 

CXTMUDE 

taro 

tcz«tcz/vol 

CXTMUOC 

lars 

text* TC XT/VOL 

CXTMUOC 

tZTA 

c 

CXTMUOC 

larr 

MCTUMN 

CXTMUOC 

lars 

END 

iKTmioe  t2«o 
iiTituoc  ta»t 
txTwuoc  laat 
cxrmioc  ttas 
txTAUDe  ia«« 
cxrmioe  tats 
txT«uoe  laM 
exTwuoe  laat 
txTAuoe  laaa 
CKTAUDC  ia«« 
exrmjDc  ia«o 
cxTmioc  ia«i 
ixtuuoc  la^a 
exTmjoe  ia«9 
axTiiuDt  ia«4 
cxr»uoc  ia«9 
txTiHiDe  ia«6 


^UNCTION  HaAII  |«X»RV«a2fRXV) 


c 

C •««a*«a««**a«**9***«a***«**«*«««*****aa*»a«»**«^a«»««**«a*««a««a«*aaaa 

C CXCCULATC  TNC  C^^CCTtve  STAAfM  XATC 

C 

aaAA«(RX*llV  iaA2Aiav«Aai#A2^lAZ-axlA*2 
AaAA«a«ASOIIT<Q«8AAeAAA3«*AXVAAa/A«|/3« 

c 


C tW  TmC  tAACCTXve  XTAAIN  AATC  IS  TOO  SMALL  COHACSMONOIMG  TO  A AIGIO 
C rLfMCNTt  SET  IT  TO  A LIMITING  VALUE* 

C ••••••••••##«•••••••«•••*««•••« AAAAAAAVAMMAAAAAVaMVAAMAAAAVMAVMMMMMMMM 

C 

IE|AOAM*LT*t*E*OAI  MOAII«l  *8-04 

MBTUAN 

ENO 


EXTIHK>8  ia«« 
cxTAUDC  taoa 
EXfAUOe  tsoo 
ExTmiOC  I SOI 
EXTAUOE  isoa 
EXTAUOC  ISOS 
BXTMUOC  tS04 
CXTWUOt  ISOS 
EXTMUOt  tSOO 
EXTMUDC  I SOT 


FUNCTION  EESTAS<SM*SZ«ST*SMZ1 

C 

C ••••••••••40«0040»M#«MOO«444««4*0«0*«*44*040M0004*0««000*0000000000000 

C CALCULATION  OE  EEEECTIVE  STOCSS* 

C •4000000000004M004440MM044VA444MM4M4444A4M4000M0004900M4MOMMOOOMOOMOOO 

c 

EESTMS* i SM^S  Z lOOZAl SZ«ST I0424C  ST-SM lOOa 
CrSTMS«SOMTI0«8OEESTOS4S**S0Z4«a) 

Mf TUC  4 
ENO 


EXTMUOC  1300 
BSTOUOC  IStO 
EXTMUOE  IStl 
EXTOUOC  iJia 
ExTouee  ISIS 
fXTfMlOB  ISI4 
EXTAUM  ISIS 
tXTOUOB  IStO 
ffXTAUDE  ISIT 
fXTAUM  IStO 
BXTAUOC  tSIO 
BXTAUOC  isao 
BXTOUOt  isai 
BXTAIM  I sat 
tXTOUDt  tsas 
BxTAuoc  isao 
txtouDt  I sat 

tXTOUOC  IStO 
EXTVUM  ISaT 
tXTWUM  isao 

txrouoi  isao 

tXTOUDt  ISSO 
tXTOMM  ISSt 

txTiHM  isaa 

tXTAUOC  ISSS 
tXTOUOt  ISSO 
IXTOUM  ISSO 
tXTOUOC  ISSO 


c 

c 

c 

c 

c 

c 

c 


< 

c 

c 

e 

c 

c 

c 

c 

c 

c 


SUOAOUTIfiC  MOOIEV  « CODE*  A,0«NUMNM,NC0»M0AN0  I 


0MO4O»<»O0O4O»4M0»O— 40— 4000»00»—»0»0»»»»000»»»»»0»<I»»»»MO»»00>OOOM 

OCTtOMINC  THE  OOINTS  EOO  XMlCM  COOOESOONOINO  COMOONCNT  OE  H I0 
OOCCtEtCO  EOUAL  TO  ZCOO  |M  AX«Ot  AND  CALL  CONOCN  EOO  MATOIX 
CONOCNSAT ION 

444— 4— 44—4  440  44  44»»»04»0»4»E4»»4»»»44»»»»4»4E»—»0440»»44  044444  40> 


COMMON  /MALL/  TMtTA*ET«TAMTM 
OIMCNSION  COOCIIItAINIOtlltOfll 

00  tat  141*  NUMMM 

tL4SMt 
IZ«tL*l 
|0*IZ*S 
C4C00CC 1 1 

4444444 4444444444444444 44444444 44 4444444444444444444444444444444444404 
ei«CM  IE  THIS  OOIMT  CONTAINS  NO  INEQNNATION  AOOUT  MtAN  StlltSt 
4444444444444444444444444444444444444444044444444444444444444444444444 


tE(C«LT*|0*l  CALL  CONOtM  A,0*MCO«MO4NO«  tL»0«l 

4444444444444444444444444444444444444444444444444444444444444444044444 
C4CCK  IE  TNC  11-Va.OCfTT  IS  SEfCIEtCD 
44444444444444444444444444444444444 44444444444444444444440 444444444 444 

|ECC*C0*l«*00*€*80*ll**04*C*CC«9*»OE«C*ta*IS*l 
I CALL  COMOCNf AtO*NEO*aOAMO*t4*0*| 
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sxrffuof 
e«T9UOE 
EXTOUOE 
Exrauos 
EXT9UDE 
cxrauoE 
EXTQliOE 
eXTQUOE 
EXTPUOE 
EXTRUOe 
EXTRUDE 
EXTR UDS 
EXT®UDB 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 


EXTRUDE  I 
EXTRUDE  I 
EXTRUDE  1 
EXTRUDE  I 
EXTRUDE  1 
EXTRUDE  1 

extrude  I 

EXTRUDE  1 
EXTRUDE  I 
EXTRUDE  I 
EXTRUDE  I 

extrude  1 
extrude  I 

EXTRUDE  I 
EXTRUDE  I 

extrude  1 

EXTRUDE  I 


EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 


CHECK  tR  the  Z-VELOCITV  1$  SRSCiriEO 


t*  (C«EC.2**0X*C*E0*I2«*0F«C*£Q.3«.CR*C«ED*19* I 
1 CAUL  CONOEN(A,B,NEO«KeANO*t7*0*» 


?SE 

356  C 
3sr  C 

359  C 
35R  C 

360  C 

361  C 

362 

363 
36A 

365 

366 
36T 
366 

369 

370 

371 


CHECK  IR  THS  RQIHT  tS  ALONG  THE  INCLINED  BOUNCARV 


XR(  C«SQ*f  ••OR«C«EO«  1S«  I CALL  BCRI  X(  A •& t NE0«  RR.ANOf  1 « THCTAI 
121  CONTINUE 

RETURN 

END 


SUeRDUTXNS  CCNOENf A,B»NE0tH9AN0»N,U» 

OERRORM  THE  matrix  CONDENSATION  RHEN  THE  VALUE  OR  A COMPONENT  X IN 
AX«D  IS  SPECIFIED  EOUAL  TO  2ERO« 


dimension  B(NED|, AiNfO, II 
DO  235  Ms2*MBAN0 

K«N*M^| 

XRfKI  235«23S»230 
230  A(K •MIsO.O 
235  AfN«M|a0«0 
A(  N,1  I ■!  ,0 
Bf N >«U 
RETURN 
END 


subroutine  BCRIx (A«etNEC«MBANO«N«THETAI 


This  subroutine  is  to  ensure  THAT  THE  velocity  along  THE  DIE  IS 
tangential  to  the  DTE  ROR  CONICAL  DIES 


CtMENStON  BfNEOI* ACNEC* t» 


SINCE  UR«U7*TAN<THETA|  ALCKG  THE  OfSt  A C0RRESO9ND INC  CHANGE  IS  MADE 
IN  THE  stiffness  equations  FOR  OCRS  AND  COLUHNS  CORRESPONDING  TO 
THESE  COMPONENrS*««THEN  THE  EQUATIONS  CONTAINING  UR  ARE  ELfMINATEO 


N2a3«N«l 
NR «N2-i 

ALRA«  tanithetai 
DO  3B0  M«I«M9AN0 

350  A(NR«Ml«AfNR*MI«ALPA 
A|NR«| lAAINR^l IFALPA 
ACNR«2I«AINR«2|*2« 

00  340  M«2,M9AN0 

KPsNR«M41 

|RfKP*LS«0)  GO  TO  341 

340  A|KR«M|«4|KRtR|RALPA 

341  CONTINUE 

ncNPt  • bcnriralpa 

DO  351  MaEfMBAKO 

K2«N2-MA] 

IRfK7*LS*OI  GC  TO  352 
A(K2«M|aACK2«M|«.A(Kr,M-l| 

352  CONTINUE 

TRI m*EQ*mbaN01  so  TC  353 
K2«NZH«-^t 

IRf KI«GT •NC91  GO  TO  353 
A<N2«M) •ACNZ«M1 ♦Af NRfPRl ) 

353  CONTINUE 

351  CONTINUE 

AC NT. llMACNZt llHAf NR*2I 
ACNRtllMl.O 
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CKTilUOe  1«1« 
CKmuoe  i^is 
CKT^UOC  14l€ 
exTftuoe  1417 
exTRUOP  1418 
exTAUOE  1418 
CXT8UOE  1420 
exTOUDE  1421 
ev74UOS  1422 
EXT4UDE  1423 
EXT4U0E  1424 
EXTRUDE  1429 


EXTRUDE  1427 
extrude  1428 
EXTRUDE  1420 
EXTRUDE  1430 
EXTRUDE  1431 
EXTRUDE  1432 
EXTRUDE  1433 
EXTRUDE  1434 
EXTRUDE  1439 
EXTRUDE  1436 
EXTRUDE  1437 
EXTRUDE  1438 
EXTRUDE  1430 
EXTRUDE  1440 
EXTRUDE  1441 
EXTRUDE  1442 
EXTRUDE  1443 
EXTRUDE  1444 
EXTRUDE  1449 
EXTRUDE  1446 
EXTRUDE  1447 
EXTRUDE  1448 
EXTRUDE  1440 
EXTRUDE  1490 
EXTRUDE  1491 
extrude  1492 


EXTRUDE  1454 
EXTRUDE  149S 
EXTRUDE  1496 
EXTRUDE  1497 
EXTRUDE  1498 
EXTRUDE  1490 
EXTRUDE  1460 
EXTRUDE  1461 
extrude  1462 
EXTRUDE  1463 
EXTRUDE  1464 
EXTRUDE  1469 
EXTRUDE  1466 
EXTRUDE  1467 
EXTRUDE  1468 
EXTRUDE  1469 
EXTRUDE  1470 
extrude  1471 
EXTRUDE  1472 
EXTRUDE  1473 
EXTRUDE  1474 
EXTRUDE  1479 
EXTRUDE  1476 
EXTRUDE  1477 
EXTRUDE  1478 
EXTRUDE  1479 
EXTRUDE  1480 
EXTRUDE  14M 
EXTRUDE  1482 
EXTRUDE  1483 
EXTRUDE  1484 


DO  399  m«2,M0AN0  ] 

XRsMR*M4S  ^ 

1RCXR«I.E«0I  6C  TO  360  1 

ACKR«M|«0*0 
360  CONTINUE 

4fNR«M|«0*0 
395  CONTINUE 

6<NZI*8(NR>4B(N2I  1 

9fNR1«0*0 

' j 

RETURN  1 

END  i 


SU8R0UT1NC  TRIA(NK«MM«A1 
C 

C ***44R<r*4R*4**44RRR4RR4r<rR4RR*RR*4RR*R*RR*RRR4*RRRRRR»4R*4444RRRR**RRR* 

C TRlANCUtARIZATXON  OR  GUASSIAN  ELIRINATIOK  FOR  THE  SOLUTION  OF  6AN0ED 

c symmetric  matrix* 

c »««RRRRRR**4R444«444««R»**R4R4R44RRF««4R**4449*«4»F*R4«RR«R4R44R4RRRR4 

c 

OXMENSXCN  AINA, 11 
1000  Nao 
too  NaN^t 

IF(N«C0«NN» RETURN 
|F(A<N,ll*eO«0«01GO  TO  100 
laN 

M8«MlN0CMMtNN-N4>t  1 
DO  260  L«2,MB 
Iat4l 

CaA(N,L l/A(N, I 1 
IFCC«eQ*0«0)6C  TO  260 
JaO 

DO  290  Mac, MB 
Jaj4l 

290  A(  I,J|aA( 1, Jt*C4AfN,K1 
ACNtClaC 
260  CONTINUE 
GO  TO  too 
END 


SUBROUTINE  BACKS f KK, MM, A, Ml 
C 

C ft»4R*RRR*«*****R4RR«RRR4*44RRR*4«4**»«*4R«4R**4*»4R*FR4RORR4R4R*«R4*44 

C BACK  SUBSTITUTION  FOR  SOLUTION  OF  PANOEO  STNMETPtC  MATRIX 

C 4*R«**4«R»4RRR«4*444R0«*4t*««4**4*R«*444Ra«4«4«*«R4*««»*R*R***R4R*4iR«R« 

C 

DIMENSION  A< 11,81 II 

MMMsMM-1 

NaO 

270  NaN4t 
CaBCNI 

|F<  A(Nl •NE*0*0IBINtaB<KI/AfN| 

XF<K«C0«NN»G0  TO  300 
iLaN+l 

lMaMtN0fNN,N4MMM| 

M«N 

DO  289  Xa|L,tH 

M«M4NN 

2ns  B( X laRC X |-A<MI*C 
GO  TO  270 
300  XLaN 
NaN-1 

IFfN.EO.OI  RETURN 

I Ha Ml NOC  NN,N4MMM| 

MaN 

DC  400  IatL,lH 

M«M4NN 

400  B<N|a8<NI«>A(M|*8<  tl 

GO  TO  300 
C 

END 

no 


exTPuos 

|4BE 

SU8RDUTINE  ePDRCe(NR«RR«P2«R5T*PPUR*B«Ma4N0tN58«NBP21 

CrTouoc 

14B7 

C 

EJfTPar'g 

t49e 

C 

RRRI 

ff XTRUOf 

1480 

c 

C4tCUL4TlCN  CX  RDRCES  OK  BCUNDART  KOOAL  POINTS 

CXTauof 

)490 

c 

*■•1 

PXTRUDC 

1491 

c 

e*T»jc? 

149E 

OIMEKSICN  N»« 1 » •RST1N8X2, l I • RPUR(Ka«21 . PR ( 1 1 • PZ 1 1 l*Bt 11 

gXTBUf5£ 

1493 

l*fNRP  .LE*  Cl  «0  TD  124 

£*TS»U?5 

14  04 

*'P*KD2«2*k84KC*  1 

■ifTa  U05 

1495 

DO  121  1»I.  N»X 

■TXTQijn* 

1496 

12«3PN«< I*«l 

E KTRU^E 

1497 

|R«12-1 

f KTQUOE 

1498 

I IZ«2Rl 

EXTRUDE 

1499 

1|R»1I7-1 

e*T»UOE 

1500 

SUP9«0« 

ExTbuas 

1501 

SUN  z«  0 • 

EXTRUDE 

1502 

DO  122  J«1 • M8AN02 

extrude 

1503 

JRsIP^J-MB4KD 

EXTRUDE 

1504 

JZbJS^I 

extrude 

1505 

txf  JR  •!»£•  01  <iC  TC  12o 

EXTRUDE 

1506 

SUNRr SUM0«PST( 1 IR • J tP8f  JR  1 

EXTRUDE 

1507 

126 

IPfJZ  .LE*  01  CC  TO  122 

extrude 

1508 

$UM2«SUMZ«.PSTCI  1Z*J)»5(  JZI 

EXTRUDE 

1509 

122 

continue 

EXTRUDE 

1510 

PR( llaSUMR^pPURCllOl 

EXTRUDE 

1511 

121 

pz(  n«suNz-ppuR(i  in 

EXTRUDE 

1512 

124 

CONTINUE 

EXTRUDE 

1513 

RETURN 

EXTRUDE 

1514 

END 

EXTRUDE 

1516 

SUSSOUTINE  IKTftPOl.  1 X.  T.  U.  X«  . YT.UUt 

EXTRUDE 

1517 

c 

EXTRUDE 

1518 

c •Rli*Ra«*««*P»«RP*»»P«PP»«PPRPPPPPPP4P«49«RR*RR*R**R*RRRRR««RRRRRR»«RR« 

EXTRUDE 

1319 

C INTERR0U4TE  THE  VALUE  DP  UU  AT  COORDINATES  C XX«VT)  6Y  KN0«IN6  THE 

EXTRUDE 

1 520 

C values  OP  U AT  PCUP  SURRCUKDXN6  POINTS* 

EXTRUDE 

1321 

c RRRRRRPVRRRPbARPPRRVAPbRRRRRVPRRRRRPARRRARRRRRRRRRRR* 

EXTRUDE 

1522 

c 

EXTRUDE 

1523 

OtPSKStCN  X(4I vXlAI ,u(4  » «A(4t4 I .COEPIA 1 

EXTRUDE 

1524 

C 

EXTRUDE 

1 525 

C P«»«PPPRPP*R««*«P*4#PPRR«9  4»PB«R«»««RM*»R«RRR*«*t«»»R««RR«RR«RRR*RRRRR 

EXTRUDE 

1528 

C DETERMINE  the  matrix  a,  TME  INVERSE  CP  The  INTERPOLATION  MATRIX... 

EXTRUDE 

1527 

C IP  IPLA6«1*  THE  A MATRIX  IS  ALREADY  KNOWN 

EXTRUDE 

1528 

EXTRUDE 

1529 

C 

EXTRUDE 

1530 

IP< IPLAC.EC.I ) 60  TC  100 

EXTRUDE 

1 ESI 

X|aX( 1 1 

EXTRUDE 

1532 

X2aX(2 1 

EXTRUDE 

1533 

X3>X| 3) 

EXTRUDE 

1534 

X48XC4) 

EXTRUDE 

1 535 

riaYd  > 

EXTRUDE 

1535 

Y2*Yf  21 

EXTRUDE 

1537 

Y3aV(3» 

EXTRUDE 

1538 

W4aV(4) 

EXTRUDE 

1539 

X12«X1-X2 

EXTRUDE 

1540 

X|3«X1*X3 

EXTRUDE 

1541 

X14*Xt-X4 

EXTRUDE 

1542 

X23«X2*X3 

EXTRUDE 

1543 

X24aX2*X4 

EXTRUDE 

1544 

X34aX3»x4 

EXTRUDE 

1545 

Y12»YIRY2RX12 

EXTRUDE 

1546 

Y13*YIPY39X13 

EXTRUDE 

1547 

Y14*Y|RY4»XI4 

EXTRUDE 

1548 

Y23«Y2«Y3wX23 

EXTRUDE 

1549 

Y24aV2RT4RX24 

EXTRUDE 

1550 

V34RY3Ry4«X94 

extrude 

15S1 

4(  l«||aaX2PV34.X?*Y244X4RY23 

EXTRUDE 

1552 

4(1 .2>*-X14Y344X3RYI4-X4*Y13 

EXTRUDE 

1553 

41 1 *31 •♦X|RY24*X2RYl44X4wyi2 

EXTRUDE 

1SS4 

A(  l«4la.xi«T234X2PY13->X3PV|2 

EXTRUDE 

1555 

A(2«1 la*y23«Y24-V34 

EXTRUDE 

1556 

A( 2»2)«4Y|3»T144Y34 

EXTRUDE 

1557 

A(2t3)«*Yl24Y|4-Y24 

EXTRUDE 

1 588 

A( 2 *4 1 m*yi 2-Y134Y23 

EXTRUDE 

15SR 

ZIIpYIR/23 

EXTRUDE 

1560 

Z12»y1RX24 

EXTRUDE 

1561 

Zl3aYlPX34 

EXTRUDE 

1562 

22IpX2*x13 

fXTRiJOE 

1 563 

Z22«Y2«X!4 

EXTRUDE 

1564 

Z23«Y2*X34 

EXTRUDE 

1565 

?3l*Y3RX|2 

Extru«5S 

1544 

Z32«T9*X14 

extrude 

1567 

ZS3«Y3«X24 

EXTRUDE 

1568 

Z41 ar4PX12 

EXTRUDE 

1569 

?42«Y4PX|3  ^21 

EXTRUDE 

1570 

ZA3RY4RX23 

6*T<»U0» 
EKT9UOC 
EXTQUOE 
EKTRUOE 
EXTRUDE 
EKTR unE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRliDE 
EXTRUDE 
E XTQUOE 
EXTRUDE 
EXTRUDE 
E XTRUOE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 


EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 


A(  3, 11  •♦«2*223«’X3«7  33«^X4«749 
A<3*2I»*K|«/| J4XJR222-X4R242 
A(3«3»«4X|*2tE-X24222^X44241 
AC  3,4|B-X|«2t  1AX2«Z2I<-X3«73t 
A(4«l >«-223AZ33«243 
4f  A,2)  •♦?ll-232«^242 

A<  4, 3>b-?12^/2?*Z4| 

AI4«4)«^21 1-221 ♦Zll 

DETER«Af 1 • 1 I^AI  1 • 2 ) ♦A  I 1 , 3 1 ♦ A C 1 ,4  I 
100  CONTINUE 
|FLA(S«0 

t^( DETER •E0*0« I GO  TO  150 
DO  130  J«lt4 
COE^IJI *Q« 

DO  120  !«1»4 

120  COE^l Jl bCOEFI J|AA< J,1 |*U(  I I 
130  COE^fJ*«COEFf Jl/OETEO 

UU«C0EF( I |♦C0EF(^|•XX♦C0EF4  3 l4VV«C0eF( 4 I4XX4VV 
RETURN 

|«0  UU-|U(1  )•( YV-V<  21 lAUl 2 14ITI  1 l-VYll/C VC  1 l-VC  21 • 
return 
END 


SUBROUTINE  b I ''A  I NS  ( R • 2 • UR«  U2  1 1 EL  • CPS  *T  EPS  t ST  R*  X « V • X ARDfUf  V»  XX*  YT, 
1 A*B« AA*RS*C«ECtNI*NJ*NTIMES«NEL  aNNP.NMAXl 

44 44444 

bUBROUTINE  FOR  DETERMINATION  OF  FLO*  PATTERN*  NET40RK  OF  GRID 
DISTORTIONS*  THE  EFFECTIVE  STRAIN  DISTRIBUTION  AND  THE  VALUES  OF 
effective  STRAIN  FOR  EACH  ELEMENT. 


dimension  RC 1 I*  Z( 1 >*URC 1 1.UZC II* lELCNEL*  t I •FPS< 5 « II «TEPS<  8tll« 

1 STRCNI *11 *XCNI *1  I * YC I I *XARB<|  • •UC N I » 1 I * VC N f * | l*XX(NTtMeSt lit 

2 YVCNTIMES*! I* AC  I I *BC I I .A A C N I * 11 *084  1 1 tCC 1 1 *CE f NTI MCS* 1 I 
DIMENSION  RRC4I * 22 44  I * URR f 4 I * UZ2C 4 1 

COMMON  /STRPATH/  STEP * Y$TART * YD I E * YE Xl T* YMI N.RENTER fRCX ITtVCX IT 

COMMON  /WALL/  THETA*FT*TANTH 

NIMlsNI-1 

NJM1«NJ-1 


FIRST  ARRANGE  THE  AVAILABLE  INFORMATION  IN  PROPER  FORM 


K«l 

L«0 

too  LvL^l 

VCL>«2(KI 
DO  110  IslfNtMl 
XC  I *LI«R4K| 

UCt  *LlsUR4K| 
V41*L>«U2<KI 

110  K«K41 

XCN|*LI«XCN|M1,LI 
U4NI*L)«U4NIMt,LI 
VCNt  tLI «V4N|M| tLl 
IF  4K.GE.NNPI  IS0«|00 


DETERMINE  THE  COORDINATES  OF  THE  CENTERS  OF  THE  CLEMENTS 


ISO  DO  160  N«1 *NCL 
I1«1ELCN«1I 
I2«ICLCN*2I 
l3«lCLrN*3t 
14«lCLfN*41 

A4NI«4Rf II I4R(I2I4R4 l3lARf 141 1/4. 
160  BCNIaC  24 1 1 1 ♦Zf I 21 42 4 1 31 42< I 41 l/4« 


ARRANGE  THE  NEVLY  OCTCRMINEO  VALUES  IN  A PROPER  FORM 


K«t 

L*0 

I6S  LFL4I 

BBCLIaSCKl 
AAI  ItLI'O. 
STR4I*LI*CPS48*KI 
DO  170  l>2tNIMl 
AAf  I«LI«A4KI 

STRCf  *Ll«CPSf8tKl 


CKTQUOe 

1695 

170  X«X4l 

exTituoe 

t6S6 

44  (N1  «(.1«REX1T4B8ILI«T  ANTH 

CXT9UDE 

1637 

lRCBB(LI*CT»V0iei  AAlKlfLlPt* 

exTPuor 

1658 

1R<e8Ct )«LT*VEXtTl  A4INI»t)«REXTT 

E VTQUOC 

1659 

8TRfNl «L)«EPSf5«K*t I 

BXT^Uf'C 

1660 

ZR(K«UT«N£L>  GC  to  165 

EXreuDC 

1661 

C 

E XTOUOC 

166? 

C 

444*44 

exTQliOE 

1663 

c 

MO«  DO  THE  INTERPOLATION  AND  DETERMINE  THE  RLOR  PATTERN 

EXT9UOE 

1 664 

c 

444444 

EXTRUDE 

1665 

c 

•XTQUOE 

1666 

DO  500  N«1 tNT  IMES 

EXTQUOE 

1 667 

c 

EXTQUDE 

1666 

c 

4*4444 

extoupe 

1669 

c 

SET  the  COORDINATES  OR  THE  STARTING  POINT 

EXTOUOE 

1670 

c 

444444 

EXTQUOE 

1671 

c 

EXTRUDE 

1672 

XX( Ktl 1 PCfNl 

EXTOUPE 

1673 

VVCNf IIsYSTART 

EXTRUDE 

1674 

c 

extrude 

1675 

c 

PRROP4 

4*4444 

EXTRUDE 

1676 

c 

determine  the  location  or  the  present  coordinates  or  the 

POINT 

IN 

EXTRUDE 

1677 

c 

TERMS  OR  THE  FOUR  SURROUNDING  NODAL  POINTS 

EXTRUDE 

1678 

c 

extrude 

1679 

c 

EXTRUDE 

1680 

iplag«o 

EXTRUDE 

1681 

NC0UNT«0 

EXTRUDE 

1682 

MOLOsl 

EXTRUDE 

1683 

LOLO«I 

EXTRUDE 

1684 

?10  NCOUNT«NCOUNT4l 

EXTRUDE 

1685 

XXXaXXl N«NC0UNTI 

EXTRUDE 

1686 

VVVaVY  t N tNCOUNT 1 

EXTRUDE 

1 687 

IRC N*£O.NTtMeS*ANO«YYY«Le*YEX ITl  GO  TO  300 

EXTQUOE 

1688 

IR< VYY«LE«YMtNl  GO  TO  300 

eXTBUDE 

1689 

DO  220  l«MOLD«NJMt 

EXTRUDE 

l€90 

220  |P( VVV.LE.rl I I.AND. VTV.6E.VI 1*1 1 1 M«| 

EXTRUDE 

1691 

DO  230  1>1«NI 

EXTRUDE 

1692 

2 30  XAPBl  I lax  Cl  •P>*(Xf  I«M>-X|  t,M41  1 >«  < Y f M l-YVYl/f  Y(N  |•Y(M•• 

111 

EXTRUDE 

1693 

00  240  tat,NlMl 

EXTOUPE 

1694 

240  IRf XXX*CE«XARe<Il*AND«XXX,LE*XAR8( 14111  Lal 

EXTRUDE 

1695 

IRC NCaUNT*EO* 11  GO  TC  250 

EXTRUDE 

1696 

IRC M«E0«M0L0*AN0«L*F0*L0L0I  go  to  260 

EXTR  UDE 

1697 

250  RRCllaKfL«Nl 

EXTRUDE 

1698 

PRC  2)*XCLtM4tl 

EXTRUDE 

1699 

RRC31aX(L4t»M411 

EX  TP  UDc 

1 700 

0PC4|ax<t41«MI 

EXTRUDE 

1701 

Z2l ll«?2C4|aY<NI 

extrude 

1702 

Z2C2|s22C3l«YfM4t 1 

EXTRUDE 

1703 

URRC 1 laUfLtMl 

EXTRUDE 

1704 

URR(2laU<LtM411 

EXTRUDE 

1705 

URRC3>«U<L4|»R411 

EXTRUDE 

1706 

URRC4|aU(L4ltNI 

EXTRUDE 

1707 

UZZU  1«V<L*M) 

EXTRUDE 

ITOb 

UZ?C21av<LtM4ll 

EXTRUDE 

1709 

UZZC3>aVfL4l«M4l 1 

EXTRUDE 

1710 

UZZ(4|aV(L4| «P| 

EXTRUDE 

1711 

60  TO  270 

EXTRUDE 

1712 

c 

extrude 

1713 

c 

444444 

EXTRUDE 

1714 

c 

the  POINT  IS  LOCATED  IN  THE  SAME  ELEMENT  AS  IT  MAS  ROR  THE  LAST 

STEP* 

EXTRUDE 

1715 

c 

HENCE  THE  INVERSE  OR  THE  INTERPOLATION  MATRIX  |$  ALREADY 

KNORN 

EXTRUDE 

1716 

c 

l•44R•4 

444444 

EXTRUDE 

I7»7 

c 

EXTRUDE 

1718 

260  IRLAGal 

EXTRUDE 

1719 

c 

E XTRUDS 

1720 

c 

EXTRUDE 

1721 

c 

DETERMINE  U AND  V AT  THE  POINT  OK  RLOM  LINE  BY  INTERPOLATION. 

EXTRUDE 

1722 

c 

•••444 

EXTRUDE 

1723 

c 

EXTRUDE 

IT*» 

270  CALL  |NTRP0L<RR«2Z*URR« XXK,YVY,UU«|PLAGI 

EXTRUDE 

1729 

IRLAGal 

EXTRUDE 

1726 

CALL  |NTRP0L<PP«ZZ«U?Z«XXX«YYY,VV«1RLAGI 

EXTRUDE 

1727 

c 

EXTRUDE 

1720 

c 

EXTRUDE 

1729 

c 

determine  The  NEXT  LOCATION  DR  POINT  ON  RLOV  LINE 

EXTRUDE 

1730 

c 

444444 

EXTRUDE 

1731 

c 

EXTRUDE 

ITS* 

XR|KtNC0UNT4| l4XXX4UU«STEP 

EXTRUDE 

1733 

YYf N«NC0UNT4| I aYYY4VV«STCP 

EXTRUDE 

1734 

LOLDaL 

EXTRUDE 

1739 

MOLCaM 

EXTRUDE 

1736 

IR'N.NR.NTtMESt  60  TD  310 

EXTRUDE 

1T3T 

XKlsRBX IT4YYCN«NCrUNT4l lATANtH 

extrude 

1738 

tR<  XXX* LT»  XX(Nt8C0UNT4| f I K X IN * NCOUNT 41 taXXX 

EXTRUDE 

1734 

IR|XKCN«NC0UNT411*GT*1*I  KXCNtNCOUNT  |la|* 

EXTRUDE 

1740 

GC  TO  310 

113 


CKTouoe 

EMTQUOC 

exTQUOE 

e XT9UDE 
EXTQUOE 
EXTPUOE 
EXTPUOe 
ExTQunE 
fXTQUOE 
EXTPUOE 
EXTRUDE 
E X TR  UOE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 

extrude 

EXTRUDE 
EXTRUDE 
EXTRUDE 
extrude 
EXTRUDE 
EXTRUDE 
E XTRUDE 
EXTRUDE 
extrude 
EXTRUDE 
EXTRUDE 
EXTRUDE 
EXTRUDE 

extrude 

EXTR UOE 
EXTRUDE 
E XTRUDE 
EXTRUDE 
EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

•XTRUDE 

extrude 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

extrude 

EXTRUDE 


THE  RQINT  IS  MOVING  M|TH  EXIT  VELOCITY 


300  XX(li.,NC0UNT41)«XXk 

YV<N,NCOUNT^l l■VVY4VE  X I TESTER 
310  1F( NC0UNT«LT«NMAX*| I 6C  TO  210 


DETERMINE  THE  VALUES  OF  EFFECTIVE  STRAIN  AT  EACH  ROINT  FOR 
LINES 


FLO« 


KCCUNT«0 

EElN»n«0« 

MOLOaLOL 0«1 
AOO  NCQUNT«NCQUNT4t 
VVV«VV(N«NC0UNT) 
XXX»XX<K»NC0UKT | 

|Ff vvv.Lfi* YEXIT)  GO  TO  400 


DETERMINE  THE  LOCATION  OF  FOUR  SURROUNDING  CLEMENT  CENTERS* 


NJMlt«NJMl*i 
DO  AtO  X«M0LD«NJM1| 

AtO  !•! VVV«LE*BD(tl*ANO«VVV.GE*0R<l4|»l  MaX 
DO  420  l«l f NX 

4 20  XARD<  I ) aAA(  I AA(  I fMl-AAf  I tM^l  ll•IOeCM|•VYV  ) / ( OBI  M |-M  |M4  f | ) 

DO  A30  lal.NtMl 

430  XF(xxX*GE*XAREf t»*AN0*XXX*LE*XARB(|4l))  L*  I 
1F<nC0UNT«E0*1)  go  to  480 
IFf m«EO*mOLO*ANC*L*EO*LOLD)  go  to  460 


INTERRCLATC  the  values  of  STRAtN*RATE$* 


450  RR< 1»  aAAlLfM) 

PP(2)aAA(L«M4l> 

RP  f 31 aAA(L4l tMTl  I 
RR<4}aAA(L4l«M> 

Z7t 1 )aZZ(4laPBIMl 
ZZf  2) a?Z<3l«Pe(M4t I 
URR^l »aSrR«L«M> 

UPRf2)a$rRCLfM4l  I 
UPR  < 3 1 a$TR  <L4t«M«>it 
UPR<41aSTRfL41tM| 

GO  TO  4T0 

460  TFLAGal 

470  CALL  INTRPOHRR.ZZ.URK*XXX,YVY*esTEP,IFLAGI 

4 4*444 Y4«a«««**44«*#4****4»*«4**a4**44444«***«444444«44 4*44*4 4*4*44 444 

AOO  THE  INTERPOLATED  STRAIN-RATE  INCREMENTALLY  TO  THE  VALUE  OF 
EFFECTIVE  STRAIN  AT  PREVIOUS  LOCATION  TO  OETEPMINE  THE  EFFECTIVE 
STRAIN  AT  NEW  LOCATION* 


EE  < N* NC0UNT4 1 1 aEE  < N« NCCUNT  > ♦CSTEP4STEP 

LOLO«L 

MOLDaM 

60  TO  400 

ABO  CCf N*NC0UNT41 )aCECN«NCOUNT) 

440  1F|NC0UNT,LT*NMAX-1»  GC  TC  400 
500  CONTINUE 


interpolate  the  effect IVE  strain  OISTRtBUTtON  FOR  ELEMENTS  FROM  VALUE 
OF  EFFECTIVE  STRAINS  ALONG  FLOY  LINES* 


00  700  Na| ,NEL 

|FCF<NI*GT«YY|1,1I1  GC  TO  580 
DO  550  Ja2,NTIMCS 
DO  520  |a]yNMAX 
XF(8(NI*L^*YY(U«tn  GO  TO  520 
L*| 

GO  TO  S30 
520  CONTINUE 

530  !F(  xx<  J«L)  •LT.ACNM  60  TO  550 

K«J 

GC  TO  8 * 0 
550  C0NTINU5 


exTRuoe 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUOE 

EXTRUOE 

EXTRUDE 

EXTRUDE 

EXTRUDE 

EXTRUOE 

EXTRUOE 

EXTRUDE 

EXTRUOE 

EXTRUOE 

EXTRUOE 

EXTRUDE 

EXTRUOE 

EXTRUOE 


060  iruasro 

RRniRXX(K-l.Li 
RR(2IRXX<K-I,L-| I 
RR<  3taXXIK,L-ll 
RR(6)aXX(K,LI 
ZZ(lt«VV(K.| ,U) 

ZZI ZlavVIX-t,L-l| 

ZZ(3taYV(x,L-|t 

ZZ<4|aVV(K,|.l 

URR(ltaEElK-t«Lt 

URR|2|aEEt«-l ,L-t I 

URRf3laEE<K,L-ll 

URR(«|aEEIK,l.l 

CALL  INTRROL  <RR • 7Z .URR . A( Ml t BfNI .TERSIStNI • lELACI 
60  TO  700 
650  TERSfStMlaO. 

700  CONTINUE 

RETURN 

END 
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